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SYMOPS i S 


Jagj iwan Lai Hod i 
M.Tech. (Hech.) 

Indian Institute of Technology, Kanpur 
September 197^ 

"DETERMINATION OF OPTIMAL ROUTES AND VEHICLE DISPATCHING 

STRATEGIES FOR DISTRIBUTION 
NETWORK" 

In the present work, an heuristic approach has been 
developed for designing the optimal routes through customers 
in a distribution network, for selecting the optimal set of 
vehicles and for their allocation to the routes designed so 
as to satisfy customers' demands at the minimum distribution 
cost. The proposed model can handle the allocation of equal 
as well as multiple capacity vehicles. This model also takes 
into account the time constraint - that has not been consi- 
dered heretofore - for ensuring the distribution of commodity 
to customers within prescribed time. 

Method of solution starts with constructing initial routes 
by Multiple Travelling Salesman Algorithm. Initial routes 
are adjusted for the capacity and time constraint, keeping 
the distribution cost at the minimum possible level. Finally, 
the refining heuristic are used to improve the routes in 
order to reduce the distribution cost. A computer package 
has been developed for the proposed methodology. 

Nine test problems have been solved using the proposed 
methodology. Same set of problems have been solved using 
Clarke and Wright's method. A comparison of the results 



indicate the superiority of the proposed methodology over 
Clarke & Wright's method i n terms of total distribution 
cost. Further, the proposed model has been tested on a 
case study based on the distribution system of Karspur 
Sahakari Milk Board. It is observed that if the routes and 
dispatching strategies obtained by the proposed methodology 
are used, it results in a saving over the currently used 
policies of Kanpur Sahakari Milk Board. Savings have been 
observed to vary from 2Q% to 3,2% depending upon the time 
period under consideration. 



CHAPTER I 


INTRODUCTION 


The primary objective of any business organization is to make 
money. Managers are always on the look out for finding ways and means 
to reduce the costs of production and distribbtion so that product 
reaches the customer at the lowest possible price. If a business is to 
succeed it must exert no less effort towards the efficient direction 
and control of distribution costs than towards the production costs. 
Therefore, most business executives are trying to find ways and means 
to reduce distribution costs. In many concerns distribution cost 
constitutes the major cost and in most industries it is a factor of 
major proportion. 


An efficient and economic distribution system is particularly 
very important for those industries in which the distribution expense 
represents a major share of the price paid by the ultimate consumer. 

When there is a buyer's market in any industry, the low-cost distributor 
has an important edge on competition. Inefficient distribution system 
leads to reductions in profit and/or dissatisfaction of customers. 

Distribution embraces all the activities in an enterprise that 
are required to move finished goods frcm the ends of production lines to 
the points of ultimate sake or use. In an industry which uses vehicles for 
the distribution of Its goods., an. attempt for getting an efficient 
distribution system creats the vehicle-dispatching problem i.e., the 
problem of designing optimal routes through customers in distribution 
network., finding an optimal size of fleet and its allocation to the routes 
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designed, so as to satisfy customers' demands at the minimum delivery 
cost. Present work attempts to solve such a problem. 

literature on vehicle dispatching problems can be traced back 
to I960 when Dantzig & Panzer (8) attempted the problem as a linear 
program . Later on, the problem was also formulated in Integer program- 
ming. terms (1) but large number of variables and constraints precluded the 
solution of such models . This 'ltd to the development of heuristic 
approaches which could successfully handle real life distribution problems 
which happened to be large in size. It is to be pointed out that most 
of the currently available approaches give near optimal solution 
the problem. Therefore, better algorithms need to be developed for 
the efficient generations of optimal solutions. Further , a thorough survey 
of the literature indicated that most of the research workers have concentrated 
on models with assumptions of constant capacity vehicles and single product 
delivery. No extensive work has been reported for the problems in- 
corporating traffic constraints in distribution network and time cons- 
traint ensuring the delivery of commodities to customer points in pres- 
cribed time. On the other hand, if the model is to be of any practical 
significance., it should also incorporate the cases of multiple capacity 
vehicles, multi-product distribution, traffic and time constraint etc. 

The literature review reveals that, but for a brief mention, no serious 
attempts have been made by the researchers to incorporate the above listed 
features into their models. In the present work, seme of these features 
have been incorporated, An approach, which is heuristic in character, has 
been developed for solving vehicle dispatching problems with vehicles 
having same capacities as well as different capacities and incorporating 
the constraints for (1) customers' demand satisfaction (2) maximum 
load on a route depending upon the capacity of vehicle on that route and 



3 


(3) route and vehicle time for distribution of conmodity to customers 
within prescribed time. The proposed methodology has been validated 
using a case study on the distribution system of Kanpur Sahakari 
Milk Board. Further 3 the proposed method has been compared with that of 
Clarke & Mright 's method (6) for nine test problems in which the time 
constraint has been relaxed. The practical application of the proposed 
model extends to all distribution systems where a fleet of delivery 
vehicles deliver a commodity from a supplying source to a set of 
customers with known location. 



C3RfiPTER-II 


LITERATURE SURVEY 


2.: 


Critical Review: 


The problem of producing routes for vehicles servicing a 
number of customers from a single depot has aroused much interest among th„ 
researchers. In general s the objective is to produce a set of routes which 
minimize the total cost of delivery subject to the restrictions due to 
the maximum load carrying; capacity of the vehicles and the maximum 
allowable elapsed time (or distance) for any route. The cost factors 
generally considered are the number of trucks and the total distance 
travelled to satisfy the customer's demands. The approaches that are 
cited in literature for solving such problems are discussed in the paragraphs 
that ensue. 

Balinski an d Quandt (1) formulated the above problem as an 
Integer Program. They listed various possible feasible combinations 
depending upon the number of customers, number of permissible geograph- 
ical routes and number of customers delivered together by a given 
vehicle and computed their costs. Problem was then to find a set of 
activities that satisfy all orders and minimizes total cost. They 
defined an activity to be a single feasible combination of customers’ 
demands. This model is not useful for large size practical problems 
because of high number of variables involved. 

Dantzig and Pamser (8) proposed a linear programming approach 
for the truck dispatching problems. The objective ms to minimize the 
interpair distances between delivery points. Their method of solution 
starts from the. basic idea to synthesize the solution in a number of 
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■Stages of aggregation' 1 in which suboptimizations are carried out on 
pairs of delivery points or groups. This method tends to lay more 
emphasis on filling the delivery trucks to near capacity than on 
minimizing total distance . Therefore s it does not give optimal solution 
and sometimes results might be far from optimal. 

Clarke and Wright (6) developed an heuristic approach to 
incorporate the practical constraints and formed the basis for most 
methods developed to date. This approach is called "Savings approach 1 ' 
because in this method the savings for each possible link are calculated 
and arranged in descending order and then each link is examined in turn 
by going down the list. If the link under consideration satisfies the 
constraints of capacity and distance, it is added, otherwise another 
link is examined, until no more links can be added and we are then 
left with the final tours. Insoite of its simplicity fdr use, this 
method does not guarantee optimality for vehicle routes. Nevertheless, 
the results obtained by this method become progressively worse as the 
constraints are made more stringent. Gaskell (15) has shown that 
the results obtained by Clarke & Wrights' method are often far frcm 
optimal results. This approach has been successfully used, for some 
cases, by Norman (28) and by National Computing Centre (29). 

Christofides and Eilon (4) generalized the simple methods of 
solution of the "Travelling Salesman Problem" for solving the truck 
dispatching problems. They used two techniques namely Branch and Bound 
technique and ‘r-optimal tour" technique. Christofides and Eilon observed 
that the Branch and Bound technique was not practicable for large size ve- 
hicle dispatching problems because it required excessive computer time 
and memory. Further, they found that the computational efficiency of 
the Branch and Bound algorithm when applied to the vehicle scheduling 
problems was substantially reduced compared with its efficiency in 
solving -an equivalent Travelling Salesman Problem. Therefore, for large 
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size problems, they developed ! r-optimal tour’ method. They defined 
’r- -optimal tour* as the tour which cannot be reduced in length by 
replacing r of its links by r other links. Their method starts 
with .an arbitrary random route structure and each route is examined 
to see whether mileage can be reduced by re-arranging the customers 
of that route. The method is then repeated using different initial 
solutions. In general, the probability that the r-optimal tour is minimal 7 
increases with the value of r. But on other hand s the amount of 
computations required increases very rapidly with value of r. They 
claim that 3- optimal tour produced, very good results. To substantiate 
this claim.. Christofides and Eilon solved ten problems and compared the 
results with the results obtained by Clarke and Wrights approach. The 
results were found to be quite comparable with Clarke & Wright's 
approach , 

The approaches discussed above s are very restrictive in 
nature because they cannot handle situations when a customer may be 
served by one of the several depots. Recently., Wren and Holliday ( 3*0 
presented an algorithm which allows routes frcm several depots to be 
constructed simultaneously subject to restrictions on number of vehicles 
at individual depots and constraints on load and distance. Their 
approach is first to construct initial feasible routes and then to 
apply a number of different refining heuristics to the initial routes. 

The quality of final routes depaids very much upon the initial routes 
constructed. For constructing .initial routes, the customers are 
presented to the algorithm in a predetermined order which affects the 
quality of initial routes. Refining heuristics further improve the 
initial routes (improvement being measured in terms of reduction of 
distance or number of vehicles) . Theoretically, the programme might be 
allowed to continue until none of heuristicscould produce an Improvement . 
However, they found that after a spectacular initial leap, the rate of 
imprcvanent often slowed down considerably. Therefore, they suggested 
that the program should be terminated after a specified amount of computer 
time. 
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2 . 2 Summary 

The review of literature has revealed a number of 
approaches developed for tackling dispatching problem. The 
limitations and the relative performances of different models 
in terms of results obtained and computational efficiencies 
have also been discussed. However, but for a brief mention, 
none of these approaches handle the problem of route selection 
and vehicle dispatching when the delivery vehicles are of 
different capacities, multiple products &r« to be delivered 
to the customers in a pre-specif ied duration of time. 



CHAPTER-111 


MATHEMATICAL TORMEATION 


3 . 1 Problem Description 

Broadly speaking- s the problem can be viewed as designing 
the optimal delivery routes and the determination of the optimal set 
of vehicles to be plied on these routes out of the available vehicles 
of known capacities . Delivery vehicles distribute the commodity from 
the depot to a set of customers each with a known location .and antici- 
pated demand. 

Before the mathematical formulation is presented s author feels 
it in order to discuss briefly the objective function 3 the various 
constraints and. the as sumptions made for the development of the model. 

3.1.1 Objective Function 

The objective is to minimize total cost of distribution to satisfy 
a known demand. The various distribution costs are : 

1) Cost of Transportation : This is the cost of fuel 
consumed by vehicles in transporting the conmodity to 
customers. Obviously the transportation cost for 

a delivery vehicle is approximately a linear function 
of distance travelled by that vehicle. 

2) L abour Maintenance and Depreciation Costs : Labour cost 
includes the wages for drivers 3 helpers and route clerks involved 
in distribution operations . This cost is a linear function 

of number of vehicles in operations because each vehicle 
is accompanied by driver , route clerk and helpers. Maintenance- 
cost is - the cost for maintaining the vehicles for keeping 
than ready for the operations . It depends upon the condition 
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of vehicle, frequency of accidents etc. Depreciation cost 
refers to the fall in value of the vehicle with respect: to 
time. 

Constraints: ■ •. - 

1) Demand Constraint : Routes designed should be such that 
anticipated demand at each customer point is fully 
satisfied. In case the management decides to fix up 
some quota to he delivered to the customers , demand is 
considered as to be equal to that quota. 

2) Capacity Constrai nt; Vehicle routes should be designed- 
such that the total amount to be delivered by any vehicle, 
on the route allocated to it, should not be more than its 
capacity . 

3) Tune Constrai nt • This constraint is of much importance 
when commodity to be delivered to customer is perishable 
one. This constraint is incorporated to restrict the 
delivery to each customer within permissible time limits. 

i) Route time constraint : The upper limit for, the time 
elapsed on a route is kept to make sure of deliveries 
to all customers during the specified time. Route time 
includes total travelled time on the route and the total 
stop - off- time at customers ; points on that route . This 
constraint restricts the length of the route and the 
number of customers to be served on a route. 

ii) Vehicle time constraint : This constraint comes into 

picture when a vehicle is assigned to more than one route 
In this case,, vehicle has to return to depot, after 
delivering one route, to get loaded for serving next 
route assigned to it . So loading and unloading of vehicl 
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at depot 'is also involved. The total time elapsed 
from the moment the vehicle is loaded ..for serving first 
route assigned to it, to the time this vehicle is 
unloaded at depot after serving the last route 
assigned to it, is the vehicle time because for 
that much time vehicle was in operation continuously. 
This constraint restricts the number of routes 
which can be served by the same vehicle. Significance 
of this constraint is when delivery vehicles are being 
operated on a shift basis. 


At this stage, it is opportune to explain sane of the terms used 

above . 

a) Capacity of Vehicle : Capacity of any vehicle is the 
maximum permissible load (in terms of weight or volume) 
which can be carried safely by that vehicle. Capacity 
of vehicle is also affected by the condition of vehicle . 

b) Stop-off --time at each customer point : Stop-off- 
time at a customer point includes 

i) time for unloading the items to be 
delivered to that customer. 

ii) time for loading the items returned from 
that customer. 

ill) time for preparing bill for money 
transaction etc. 

In General s stop-off time at each customer point is 
different, depending upon- the amount of loading and 
unloading involved. 
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c) S top- r ff-tlme at depot : It includes the time for loading 
and unloading of the vehicle at depot. 

d) Maximum remissible route time ’ Route time is the elapsed 
time for which vehicle is on that route, It includes 

i) travelling time on the route 

ii.) stop off -time for all the customers on the route. 

r-feximum limit for route time is put for making the deliveries 
tc all customers in a specified time. 

e) Maximum permissible vehicle time : It is the time for which 
the vehicles are available for the distribution. It 
includes 

i) loading and. unloading time for the vehicle at depot, 

ii) vehicle idle time. 

iii) route time for the route(s) assigned to that vehicle. 

f ) Average speed of vehicles Speed, of vehicle depends upon 
its size aid condition. In addition to it, average speed 
of the vehicle depends on the traffic intensity of the route 
assigned to the vehicle. 

3.1.3 Basic Assumpt io n s: 

1) All delivery vehicles would return to the depot after 
distributing the commodity to customers. This is usual 
practice of most of the companies engaged in distribution. 

2) Only one route passes through a customer and full delivery 
to that customer is made by the vehicle assigned to that 
route passing through that customer. This assumption is 
valid for only those customers whose demand is less or equal 
to the capacity of vehicle. In fact, most of the 
distributing agencies follow this strategy. 
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3) All vehicles are in sane condition. This assumption implies 
that maximum allowable load to be carried by vehicles 

of the same capacity is same. 

4) Average permissible speed on the routes is same for a 
vehicle throughout the distribution network. 

5) Supply from source i.e. depot is sufficient to fulfill 
demands of all the customers. 

6) labour and maintenance cost is same for every vehicle and 
therefore, the total labour and maintenance cost is considered as 
a linear function of the number of vehicles. 

7) There is no vehicle idle time at the depot due to delay 
in loading and unloading of vehicle i.e., vehicle is ready 
for loading or unloading as soon as it reaches the depot. 
Similary assumption has been made that no vehicle idle time 
occurs on the route allocated to it. 

~ ' Assumptions (33- and (4) follow that average speed of vehicles 
is no more a function of condition of vehicle or traffic 
intensity on routes. It is only a function of type of vehicle 
and its size. 

A. mathematical formulation of the above stated problem has been 
•presented in the following section. 

3.2 r fathematical Formulation 

The notations and terminology used for the development of 
mathematical model is presented below: 

3.2.1 Notations and terminology 

NC - Number of cus toners to be served. 

N ~ Total number of nodes in distribution network. 
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i 

J 

G 


S 

Xg 

NV 

d ij 

r 


y ij = 



% 

6 r = 

g 




Index for the preeeeding node 1 = 1,2, ...» N 
Index for the suceeding node, j = 1,2, ..., N 
Number of routes 

Index for the routes, g * 1,2, ..., G 
gth route 

Number of vehicles employed for the distribution 
Shortest distance between ith & jth node 
Index for the vehicles, r = 1,2, . .., NV. 

Capacity of rth vehicle 
Demand at jth node, j = 2,3, ..., N 
D]_=0 for j = 1, i.e., for the depot 
1 if ith node is connected to jth node in network 
0 otherwise 

Set of delivery point forming the routeAg 

i.e., jgl x f' such that j is on routed; for j^l. 
g 

Number of delivery points forming the routeAg, 

1 if rth verhiele is assigned to routeAg 
0 otherwise 

Stop off time at jth customer j = 2,3, ...» N 
Loading and unloading time for rth vehicle at depot. 
Cost of transportation/unit distance for rth vehicle 
Maximum permissible route time. 

Maximum permissible vehicle time 
Average speed of rth vehicle 
Labour and Maintenance Costs per vehicle 
Depreciation cost of rth vehicle. 
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3.2.2 Objective Function: Objective function is to minimize 
distribution costs. Total cost of distribution is determined 
as follows 


1) Cost of transportation 

Total distance travelled by vehicle on gth route* 


’ C j e I- dl 3 yi 3 + ? ,„r d iJ y i3 + Lt 


i,JeI 


Ag 


j el 


Ag 


( 1 ) 

Total cost of transportation for gth route = 

NV r , 

TC > = E V^W 

8 r=l g s r 

Therefore total cost of transportation for all 
the routes 

G NY 

TC X = E E 6^(TD k .K ) (2) 

g=l r=l g 8 T 

2) Labour and Maintenance Cost : Since the labcbir 
and maintenance cost is assumed to be the same 
for all the vek&lles, the total labour and 
maintenance cost can be expressed as follows: 

TC g - NV * P (3) 

3) Depreciation Cost : Let DC r be the depreciation 
cost for rth vehicle. Then the total depreciation 
cost is given by 

NV / 

TC- = E DC 
3 r=l r 




( 4 ) 
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The total cost of distribution is the sum of transportation 
cost, labour maintenance and depreciation cost. Therefore, 

Total cost = TC^ + TCg + TC^ 

c* MV N. V 

= 2: S (~y Dc: K^. "*) -f N V P ■+* ^ 

3 ' ^ 

where TDg is distance travelled by vehicle on gth route. 


3.2.3 Constraints . 

1) Except for depot, each one of other nodes is connected to 
one and only one other node. Mathematically this can be 
expressed as t 

N 

i) ^ yij - l for j = 2,3, jN (6) 

{ - % 

ii) 4. yi|'= ^ for i = 2,3, ,N (7) 

2) Since each route is to be assigned to only one vehicle, 
we obtain following constraint: 

c ^ , 

22- O c, ~ i for g = 1,2, ...... ,G (8) 

J 

However, a vehicle can be assigned to many routes. Therefore, 

<5- 

X *2 , 

'■y- \ <) ~~ 0 ^ -Where K is an integer value ... ( 

S 

3) Demand Constraint : The restriction that one route should pass 
through a. customer and full requirement of the customer is 
delivered at a stretch, implies that the capacity of the 
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vehicle assigned to the route 5 on which the considered 
customer is, should be equal to or more than demand of 
that customer, 

Kty 

i.e, S r C ^ S ^ hT D ) for j = 2 S 3, • • • • jN (10) 

. P-j - ■ i ( \ v 

1 J , 

where j(- 1 

Ao 

Constraints (3) & (4) imply that if rth vehicle is assigned 
to the route., say gth 3 on which the customer ,j is 3 then 

C > D. (11) 

■T J 

4) Capacity Constraint: If rth vehicle is assigned to gth 

route then capacity of rth vehicle should be equal to or more 
than the sum of demands of all customers on the gth route. 



1) Koute time constraint : This constraint restricts the length 
of the route and number of customer to be included in that 
route, such that delivery to each customer is made within 
specified time . Route time of any route is the sum of the total 
time . travelled by the vehicle on the route and the total 
stop-off -time at all the customers on that route. Let 
be the maximum permissible route time. Then this 
constraint implies that route time of any route should 
not exceed T^. 
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Total distance travelled by vehicle on gth route , 

'TOg-C ■2a l jy 1 j + 2 z J-ji-yjil ••• cuo 
LJtI A 3 i)i£l X3 ^>3 J 

If rth vehicle is assigned to gth route * i.e. s rv =1 

” O a 

then route time constraint is 


(TD /v ) + ■< t . 

g r <-• sj 


<r 


... (15) 


J (? 1 


AS 


In general,, route time constraint is 


-+ - z > s J 

2) Vehicle time constraint : This constraint implies that a 
delivery vehicle should make all the distribution 
operations assigned to it within a specified time i.e. 9 
maximum permissible vehicle time, so that they are 
available for the next operation. Vehicle time includes 
the route time, for the route allocated to it and loading 5 

unloading time of vehicle at depot. Mathematically, this 
constraint can be expressed as. 




<C T 


g=l,2. . . .G 


(16) 



CF f 

-g: 

f) 


7-,°3 L^ t T)/' 0+-_ i sj + i t 

' *3 


for r=l,2. . .NV 


Where TD., is total distance travelled by vehicle. 
6 ) Node-connecti^ l ty Constraint : 


... (17) 


y. = 1 if ith node is connected bo jth node , 


for i = 1,2, N: j = 1,2,, 


..N 


0 otherwise 


(18) 



Equations (5) through (18) represent the mathematical 
f emulation of the route selection and vehicle assignment problem 
under consideration- A solution methodology for the problem is 
presented in the next Chapter. 



CHAPTER 4 


METHODOLOGY 

In this chapter, the methodology for solving the 
problem .formulated in Chapter 3 S is presented. The approach 
is heuristic in character and gives near optimal solution. 

Methodology for the solution of the problem broadly 
consists of following phases: 

(1) Determination of the shortest distances 
between the nodes in distribution network 
using shortest distance model. 

(2) Selection of the optional set of vehicles 
out of available vehicles. 

(3) Design of optimal routes and the optimal 
allocation of vehicles to the various 
routes by Route design model. 

4 . 1 Shortest Distance Model 

Shortest path algorithmic) has been used to find the 
shortest distances between the delivery points, in the distri- 
bution network. The procedure is to examine simultaneously 
all the routes out of the starting point i and into the ter- 
minal point j as far as their adjacent connected points and 
to extend further the routes which have so far covered the 
least distance. This process is repeated step by step until 
a route out' of i has’- a delivery point on it which has already 
occurred on a route Into j or vice-versa. All such complete 
routes between I and j are examined and shortest possible 
route between i and j is found. Similarly shortest path 
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distance paths for each pair of nodes are determined. Computer 
program for the above algorithm has been written in FORTRAN IV 
language . 

4 . 2 Routes Design Model 

Looking back into the mathematical formulation presen- 
ted in Chapter S» it is observed that objective function, i.e., 
total cost of delivery depends upon (i) total distance travelled 
by vehicles, i.e., the sum of distances of all the routes for 
distribution, and (ii) the number and types of vehicles used 
for distribution. For a given set of vehicles, the total deli- 
very cost depends upon the total distance travelled by vehicles 
in the set. So, for a particular set of vehicles, the problem 
consists of designing the routes for the vehicles for minimum 
total distance and the allocation of the vehicles to the designed 
routes for minimum delivery cost. As the set of vehicles is 
chosen from the available vehicles for distribution, there will 
be few feasible sets of vehicles to be considered. Problem is 
solved for each set of vehicles and the one which yields best 
results is chosen. Main steps involved are as follows: 

(I) Choose a particular set of vehicles to be used for 

distribution . 

(ii) For the above chosen set of vehicles, design the 

rpyfe.ts such that the total distance is minimum. 

Designing of routes is carried out in three stages: 

(a) Constructing initial routes 

(b) Adjusting initial routes for satisfying the 
constraints . 

(c) Improving the adjusted routes by refining 
heuristics , 

(iii) Total cost of delivery is calculated. 

(iv) •''PShpeat the Steps (ii) and Step (iii) for various 

feasible set of vehicles to be used for distribution • 

■ ... ■ - .... : ... 

and select the best solution yielding the minimum 
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Next, we have to find feasible sets of vehicles which 
can be used for distribution, 

4.2.1 Determination of feasible ’sets* of vehicles which 

can he used for distribution 

Initially, It is assumed that each vehicle is allocated 
to one route only. Later on, in the refining heuristics, this 
assumption is relaxed and if possible and economical, one vehicle 
Is allocated to more than one routes. With this assumption, 
the number of vehicles will be equal to the number of trips which 
naturally depend upon the capacity of the vehicles and the total 
delivery at all the delivery points. 

In case the vehicle set comprises of vehicles of same 
type and capacity, then the problem is to determine optimal 
number of vehicles to be used for distribution. For obtaining 
the lower bound on the number of vehicles to be used, following 
formula can be used: 

n 

Number of vehicles = N = £ D,-/C 

j=2 J 

where Dj -»■ Load requirement at jth delivery point 
C -* Capacity of the vehicle 

at j =1, Demand is zero because it is supplying source 

If N is not a whole number, it is rounded off to next 
higher whole number. The problem of designing the routes may be 
solved for several values of N and the best solution is chosen 
It has been experienced by solving some test problems that usually 
no more than throe values of N need be considered (starting with 
lowest fsdSsible value of N and then increasing it in steps of 
one) for the final solution to be obtained. 

Above formula holds good only in the case of same 
capacity vehicles. For different capacity vehicles, the cost 

of delivery, in addition to number of vehicles, depends upon 
the type of vehicles. 'Therefore, in such cases, the feasible 
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sets of vehicles for distribution are determined and problem 
may be solved for each of these sets to achieve the best sclU' 
ticn. Vehicle-set-selection can be done by one of the follow- 
ing strategies; 

(i) VehMle-set -selection policy 1: Vehicles are 
picked from the available set of vehicles, 
for assignment, in descending order of their 
capacities . 

(ii) Vehicle-set-selection policy 2: Vehicles are 
picked from the available set of vehicles, 
for assignment, in ascending order of their 
capacities . 

(iii) Vehicle-set-solution policy 3: Vehicles are 
picked randomly from the available set of 
vehicles. 

The minimum number of vehicles to be employed is r, such 

that 

r-1 n 

I C, < E D - 

i=l j =2 J 

and 

r n 

S C ± > E D 1 . 

i=l ~ .1=2 

where Dj load requirement at jth customer 
-*■ capacity of ith vehicle 

Mathematically the various set selection policies can 
be represented as follows; 

Policy I ; C-j. > C i+1 > C i+2 

Policy II : Cj_ < C^ + - L < C^ +2 

Policy III : Cj_, selected randomly. 
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For each policy, problem is solved for r+s vehicles, where 
s = 0,1,2, . .., S; such that S+r = total number of available 
vehicles. 

4.2.2 Designing the Routes 

Having determined the sets of vehicles to be used, 
the problem of designing the routes is solved for each of 
these sets as follows: 

(i) Construct the unconstrained initial routes, 
for minimizing the total distance travelled 
on the routes by vehicles. 

(ii) Adjust the Initial routes for satisfying 
constraints, increment in total distance, 
due to re-arrangement of customers on routes , 
is kept minimum. 

(iii) Attempt to accomodate ommitted customers. 

(iV) Improve the routes obtained by refining 
heuristics . 

4.2.3 Construction of Initial Routes 

Initial routes, one for each vehicle, are constructed 
to cover all the delivery points in the distribution network 
such that the sum of total distance of all the routes is mini- 
mum. IT hile constructing initial routes, constraints are not 
imposed. Initial routes are developed by using Multiple 
travelling salesman algorithm. Multiple travelling salesman 
problem is the extension of single salesman problem. 

In general, the objective of single travelling salesman proble 
is defined as to find a minimum cost (or distance) route for 
a salesman to visit each customer once and only once, starting 
from depot tr coming back to depot. However, if the real 
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depot is eliminated and replaced by NV artificial depots 
(where NV is the number of vehicles used for distribution), 
all located in same position, then problem of constructing 
initial routes is similar to multiple travelling salesman 
problem and reads' Given NV vehicles and NC customers, find 
NV routes such that every customer is visited exactly once 
by exactly one vehicle, so that the total distance travelled 
by the vehicles is minimum. ! An efficient algorithm for 
solving such problem is M-salesman travelling salesman algo- 
rithrrr . Initial routes have been obtained by this algo- 
rithm as follows. 

As stated above, real depot is replaced by NV artlfic. 
depots. Thus, 

NN = NC + NV = (N-l) + (NV) = N + NV - 1 
where N is total number of nodes in original network 

NN is total, number of nodes after replacement. 

NV is number of vehicles. 

The new distances [C^] are obtained from the original distances 
[dj_j] by augmenting the matrix [d^j ] with (NV-1) new rows and 
columns, where each new row and column is a duplicate of the 
first row and column of the matrix [dj_j] (because first row $-.■ 
column correspond to depot). All other new elements of the 
augmented matrix are set to infinity. 

Using the notations used in Chapter 3 S the problem 
for constructing initial routes can be structed similar to 
Miller-Tucker-Zemlin (27) *s- formulation. Thus 

NN NN 

min Z = I E C 11 y., (1) 

i=l j=l 0 

where Z is the total distance of all the routes. 
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Subject to 


NN 







I 

i=l 

y ij 

= 1, 

3 = 

a e a 



NN 







1 

j=l 


= 1, 

i = 

1 2 

NN 

(2) 

y_> • 

= 0 

or 1 for 

all 

i.sj 

e <s 

(3) 

x. - 

i 

X j 

+ [y~f-]y. , 
L N\r j « / X; 

! i £ 

+ m ■ 

-2 

.... (4) 


for all i ^ j and i,j ^ 1 0 y l Q = (1,2, . . . ,NV) 

The variables x-j_ and Xj are arbitrary real numbers which 
satisfy the constraint (4). The braces [.] denote the 
largest integer not to exceed the value within the braces. 

In above formulation. Equations (1), (2), (3) represent 
assignment problem formulation. The constraints (4) are the 
loop constraints which block the formation of infeasible 
subtours. Subtour is feasible only and only if it contains 
at least one artificial depot. 

Stated in broad terms, the algorithm works as follows: 

An assignment problem with constraints (2) and (3) is solved. 
If the solution satisfies constraints (4), solution obtained, 
i.e., the initial routes constructed are optimal. If one or 
more of the constraints (4) is violated, a subset of the 
constraints (4) is implicitly introduced by setting certain 
elements of the distance matrix to infinity. The new problem 
is another assignment problem, which is solved and so on. 
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Hence, the algorithms for generating initial routes, solves 
a sequence of assignment problems, which ends x^hen the 
solution to one of the problems satisfies* all the constraints 
(2), (3) and (4). The sequence of problems is defined in 
such a- way as to assure that the terminal solution is optimal. 

Briefly, the steps of the algorithm are as follows: 

Step 1: Solve the- associated assignment problem. If all 
subtcurs are feasible - terminate - the solution 
is optimal. If not, set Zgg equal to the value of 
assignment problem. Construct a feasible solution 
by initial heuristics, set Z-qb (upper bound) equal 
to the value of the feasible solution and go to 
Step 2 . 

Step 2: Locate .the shortest infeasible subtour in that 
assignment problem solution which defines Z^g. 
Generate and solve k subproblems, one for each arc 
in the infeasible subtour. If the value of a sub- 
problem Z >_ Zub, delete it from the set of sub- 
problems, S, to be considered further. If Z < Z UB 
and the subproblem solution is feasible, redefine 
Zyg, delte the subproblem, continue. 

Step 3: If set S is empty, terminate, otherwise scan the 
set S for that problem such that Z = min Z± and 

' i es 

set Z LS = Z. If Z LB = Zjjq - terminate - otherwise 

go to Step 2 . 

4.2.4 

In the previous section, unconstrained initial routes 
irere obtained for the minimum total distance. These initial 
routes are now adjusted, to satisfy the constraints, by re- 
arranging the customers on the routes. Re-arrangement of the 
customers on the routes involves the increment in the total 


Adjusting the initial routes for satisfying 
constraints 



27 


distance of initial routes. Proposed algorithm s for adjust- ' 
ment of initial routes, attempts to re-arrange the customers 
routes for adjustment, such that the increment in total 
distance of initial routes is minimized. 

Algorithm starts with allocating the vehicles to the 
initial routes by following policies: 

(a) Vehicle allocation policy 1: Allocate vehicles 
to routes in random order. 

(b) Vehicle allocation policy 2: Allocate vehicles 
to routes in descending order of capacities, i.e., 
maximum capacity vehicle is allocated to first 
route and minimum capacity .vehicle - ■*» ■ last route. 

(c) Vehicle allocation Policy 3: Allocate vehicles 
to routes in ascending order of capacities. 

For each policy, allocation of vehicles to routes causes some 
routes to be overloaded* and some to be underloaded*. An 
attempt is made to accomodate customers from the overloaded 
routeS'-on- underloaded routes for minimal increment in distance 

-three. 

due. to reshuffling of customers. All the ^jtr-ee policies listed 
above, are tested and the policy giving the best solution is 

chosen. For each of the policies, the routes are adjusted 

*(1) Overloaded Route : If the "total demand of customers on a 
route is more than capacity of vehicle allocated to it, 
is more , than route is called an overloaded route. 

(2) U nd e r 1 o a d e cl Route : If total demand of customers on a 
route is less than capacity of vehicle allocated, then 
route is called an underloaded route. 
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as follows: 

Step ^l : Firstly the proposed algorithm lists out all over- 

* 

loaded and/or overtimed routes* and adjusts one 
such route at a time. An immediate problem arises 
as to which overloaded route should be considered, 

because the order of presentation of overloaded 
routes might affect the solution. In order to 

observe the affect of order of presentation of 
overloaded routes for adjustment to algorithm, 
the following three route policies have been consi- 
dered: 

Route Policy 1: Overloaded route is selected randomly 

out of the overloaded routes. 

Route Policy 2; Most overloaded route is selected first 
Route Policy 3 : Least overloaded route is selected 

first . 

Overloaded route, obtained by a particular route 
policy, is adjusted as follows. 

Step 2: Suppose the overloaded route to be considered for 

adjustment is OVR^. Each customer on route OVR^ is 

considered to be discarded from route OVR^ and to 
be accomodated in some of other routes. Suppose 

customer A-, is considered to be discarded from route 

_L 

OVRj . This customer can be accomodated in any of 
other unadjusted routes and also to already adjusted 
routes if inclusion of it does not violate the capacity 
and time constraint of that adjusted route. Now a 
search is made to find the pair of linked nodes B]_ 
and C-^ on any of these permissible routes for 

accomodating A^, such that new connections of A-j_ 
with and involve the minimum distance Increment 
MD^, MD]_, B^ and Cj corresponding to discarded custom 

aC*. — — . « , . • - ------ -- " 

* Overtimed route is the route for which route time is more tfrayi 
specified maximum route time . 
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Step 3: 


are noted and customer is restored to 
route OVR ± . Then another customer A 2 is consi- 
dered. Similarly for all the customers Aj on 
route OVR , search is made to find minimum 
distance incdrement MD. and corresponding new 

J 

connections for A^ , i.e., the pair of linked 

nodes 3 . arid C4 . 

J ' d 

In this step 3 a combination of customers on 
route OVRj_ is sorted such that (i) removal of 
these customers adjusts the overloaded route 
under consideration, (ii) sum of minimum 
distance increments of discarded customers is 
minimized. This problem for sorting the best 
combination of customers, has been tackled by 
zero-one programming. Problem can be expressed 
in mathematical form as follows: 

Let J-.set of customers on overloaded route OVR-- 

MDj -minimum distance increment involved with 
removal of jth customer, where j e J. 

OVC(i) - the overloaded quantity on route OVR^. 

QVT(i) - the route time exceeded for route OVR-j 
D,- - demand of jth customer, where j e J 
t j - route time saving due to removal of jth 
customer from route OVRj_. 

X -4 - 0,1 variables. 

J > 

X • = 1 if jth customer is discarded 
J 

= 0 otherwise . 

The combination of customers is to be chosen 
from the set J so as to. 
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minimize E MD. X. (1) 

jeJ 3 3 

subject to I t- X. > OVT ( i ) 
j eJ. 3 3 

Z D, E, > OVC(i) . . (2) 

jeJ J ° 

X. = 0 or 1 ...... (3) 

Zero-one programming solves the above problem and gives the 
optimal combination of customers to be discarded from route 

0VRj_ in order to adjust it. The discarded customers are 
connected to bh'eir new edHheetionS Bj and Cj as found in 
Step 2. If any customer j on route 0VRj_ can't be included 
in any other route 3 MD-? is set to infinity so that this does 
not come into the solution giving the combination of custome? 
to be discarded. Such customers are kept aside and listed 
as ommitted customers. 

Thus Step 3 adjusts the overloaded route OVR^ . 

This route is now called adjusted route for capacity and tinr 
constraint and in further adjustments, customer can be 
accomodated to this route if and only if that customer does 
not violate the capacity and time constraints of this route. 

Repeat the Step (1) through Step ( 3 ) till all the 
routes are adjusted for capacity and time constraints. 
Ommitted customers are also listed to be considered in next 
section. 
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4.2.5 Consideration of Ommitted Customers 

Next, an attempt is made to accomodate the ommitted 

customers into some of the routes. Each ommitted customer 
is taken in turn and every route is inspected in order to 

determine whether by removing another customer, the ommitted 

one could be accepted; the customer removed .is then fitted, 

if possible, into another route. If there is no customer 

left unassigned, we go ahead for improving the solution by 

refining heuristics. However, if some ommitted customers 

remain, we go back .to Section 4.2.2 for considering another 

set of v ehicles to be used for distribution. 

4.2.6 Refining Heuristics 

Routes have been improved further for reducing total 
delivery costs by the following refining heuristics. 

(1) In first heuristic, each customer is considered in 
turn and the possibilities of moving him to all other positions 
in the same or another route are considered. He is moved to 
the position found for which the total distance is reduced and 
no constraints are violated. 

(2) Second heuristic also attempts to reduce total distance 
hi customer from one route is placed in another route, if a 
customer from the second route could be inserted elsewhere in 
the system to reduce total distance. This procedure was 
developed because it was felt that the first heuristic might 
fail to place a customer in a route where this was desirable 

because of violation of capacity constraint. In such cases. . 
the room could have been made in the route by removing another 
customer from it by this heuristic. 
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(3) Lastly, a refining heuristic is applied for reducing 

the number of vehicles In use and hence reducing the total 
delivery cost. If some vehicles return to depot (after dis- 
tribution) after only a few hours, it might be desirable to 
give a second trip to one or more vehicles if second trip can 
be completed within maximum allowed vehicle time. This pro- 
cedure examines every pair of routes and searches for the 

pairs whose combined duration, after allowing for time to 
reload the vehicle, is below the maximum permitted time. 

This process does not reduce the number of distinct routes, 
therefore no savings in terms of total distance involved occur 

But If two routes are combined and allocated to cuaevehicle, 
if possible, a vehicle is saved which in turn reduces the 
total delivery cost, 

The computer programmes have been developed for 
all the algorithms discussed above. The computer package 
has been written in FORTRAN IV language and has been run on 
IBM 7044. 



CHAPTER - V 


COMPARISON OF PROPOSED MODEL WITH CLARKE & WRIGHT'S METHOD 


In this chapter, the proposed model has been compared with 
the Clarke & Wright's (6) method for solving the vehicle dispatching 
problem. Mine test problems have been structured and solved by 
both the methods and the results obtained have been compared. Since 
Clarke & Wright s method, does not take care of time constraint, the 
time constraint has not been considered in the test problems. The 


various costs of distribution arc assumed to be the same for all 


the test nroblams. The various costs of distribution are tabulated in 
Table 5-1. The data on shortest distance between customers , customers' 
demand and the capacities of the vehicles for the various structured 
problems are riven in tables 5-2 through 5-10. The shortest distance 
matrix has been generated by generating psuedo - random numbers between 


0 and 1. The psuedo ~ random numbers have been generated * : 

by using the function KNDY5'“ . The shortest distance between the nodes 

1 & j is determined using the relationship (d„ ) = (A * ENDY5) where A 
is some arbitrary chosen real number. For each problem a different 
value of A has been selected. Input data and results obtained from 
the Clarke & Wright's method and proposed model, for each test problem. 


are given in Tables 5-2 through 5*10. Results indicate that in a.ll the 
test problems proposed, model has performed better in terms of total cost 
of delivery. A summary of the results obtained by the two methods 
is given in Table 5.11. 


The results have been summarized in following paragraphs. 

The results indicate that 

i) For each problem the total distance travelled by vehicles 


« Note: ENDY5 is a built in function in JJW 7044. 
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to satisfy the demands of the set of customers 
considered in problem, is less in the case of 
the routes obtained by proposed method, 

ii) The number of vehicles required for distribution by 
proposed method has been observed to be either equal 
to or less than found by the Clarke & Wright's method. 

In test problem numbers 2,3 and 9. the saving of 
a vehicle has been observed with proposed method. 

iii) The cost of transportation is less with use of proposed 
method in all the test problems. Percentage saving in 
the cost of transportation with proposed method over 
Clarke & Wright’s method varies from 6 %, (for test problem 1) 
to 32 $ (for test problem 5). 

iv) Total cost of delivery has been observed to be less In 
case of distribution by proposed method than that in 
case of the Clarke & Fright’s method in all the test 
problems. Percentage saving In the total cost of 
delivery varies from 2 % (for test problem 1) to 23 % 

(for test problem 2). 

v) Computational efficiency of the two methods has been 
compared in terms of execution time and caopgiter 
memory' heeded. Table 5.11 shows the execution time 
ler both the methods. Execution time for the proposed 
method has been observed, to be higher than that for 
the Clarke & Wright’s method for all the test problems. 
Computer memory requirement is also more for the oreposed 
method than that required for Clarke & Wright’s method. 

To conclude, the proposed method is not computationally as 
efficient as Clarke & Wright’s method- but this is more than made up by 
the fact that the proposed method results in considerable higher 
savings in cost to the tune of 10 to 20$ and hence, is justifiably 
preferable . 



Table 5.1 


1) Operations Characteristics of vehicles 
1) Average speed = 20 fin. /Hour 

11) Consumption of fuel (netrol) 0 . 25 litre/km. 
ill) Fuel cost = Rs. 3.27/litre 

2) Cost of transportation ~ Rs. 0.8l/fin. 

3) Maintenance cost per vehicle- = Rs. 400 per month 

4. Depreciation cost : 

Initial cost of vehicle = Rs. 80,000 

Depreciation rate (straight line method) = 10$ per year. 

Depreciation cost per vehicle = Rs. 666.6 7 per month. 

5) Labour cost : 

Salary to driver = Rs . 200/ ■ per month 
Salary to helper = Rs. 175/- fher month 
Salary to Route clerk = Rs. 250/-- per month. 

Total Labour cost per vehicle = Rs. 800/ - per month. 
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T A B L t; 5 . 2 


TbST PROBLEM 1 


iNPUT INFORMATIONS 


NO* OF CUSTOMERS = 10 

DISTANCE MATRIX GENERATED BY USING RNDY5 
8(I,J) » I FIX (A. RNDY5( Y 5) ) 

A = 17.5 

DEMANDS FOR CUSTOMERS 
10 9 15 7 6 12 11 10 9 6 

CAPACITY OF VEHICLES = 50 


RESULTS 


CLARKE 

AND WRIGHT METHOD 



ROUTE 

UNUSED CAPACITY 

ROUTE PATH 


1 

C 

16 2 3 

5 11 

2 

5 

1 b 4 10 

9 1 

OPTIMAL 

NO. OF VEHICLES = 2 




TOTAL DISTANC£= 67 

COST OF TRANSPORTATION = 164-1.50 

MAINTENANCE, DEPRECIATION AND LABOUR COST= 3733.32 
TOTAL COST = 5374.82 

PROPOSED METHOD 


ROUTE UNUSED CAPACITY ROUTE PATH 

1 4 1.9 7 10 6 1 

2 1 I A 2 5 11 8 1 

BEST ROUTE POLICY IS 1 

OPTIMAL NO. OF VEHICLES = 2 

TOTAL DISTANCE 63 

COST OF TRANSPORTATION = 1543.50 

MAINTENANCE, DEPRECIATION AND LABOUR COST= 3733.32 
TOTAL CQST= 5276.82 
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T A B L t 5. 3 
TEST PROBLEM 2 


NPUT INFORMATIONS 


NO. OF CUSTOMERS = 12 

DISTANCE MATRIX GENERATED BY USING RNDYS 
B ( I * J ) = I FI X( A .RND Y5 ( Y 5) ) 

= 10.0 

DEMANDS FOR CUSTOMERS 

23 25 24 18 25 2D 24 22 18 19 20 17 26 

CAPACITY OF VEHICLES = 100 


RESULTS 


CLARKE 

AND WRIGHT METHOD 





ROUTE 

UNUSED CAPACITY 

ROUTE PATH 




1 

76 

1 4 1 




2 

15 

16 9 

5 

7 

1 

3 

• 13' 

1 N b 11 

14 

10 

1 

. 4J. 

15 

1 12 2 

13 

3 

1 


OPTIMAL NO. OF VEHICLES = 4 

TOTAL DIST ANC£= 53 


COST OF TRANSPORTATION = 1298.50 

MAINTENANCE, DEPRECIATION AND LABOUR CQST= 7466.64 
TUTAL COST = 8765.14 

PROPOSED METHOD 


ROUTE 

1 

2 

3 


UNUSED CAPACITY 

ROUTE PATH 




12 

1 12 8 

5 

14 

1 

7 

13 2 

7 

6 

1 

0 

1 11 13 

4 

9 

10 1 


BEST ROUTE POLICY IS 1 

OPTIMAL NO. OF VEHICLES « 3 

TOTAL DI ST ANCE= 48 


COST OF TRANSPORTATION = 1176.00 

MAINTENANCE, DEPRECIATION AND LABOUR COST= 5599.98 
TOTAL COST = 6775.98 
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T A 8 L E 5. 4 


TEST PROBLEM 3 


INPUT INFORMATIONS 


NO. OF CUSTOMERS = 15 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B( I,J) = IFlXt A.RNDY5(Y5) 3 

A = 22.5 

DEMANDS FOR CUSTOMERS 

15 12 11 12 8 18 17 11 12 10 9 18 17 8 10 

CAPACITY OF VEHICLES = 50 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE 

1 

2 

3 

4 

5 


UNUSED CAPACITY 
6 
5 
9 
40 
2 


ROUTE PATH 

1 ' 6 7 13 1 

1 8 14 4 1 

1 9 16 15 10 1 

1 11 1 

1 12 5 3 2 1 



COST OF TRANSPORTATION = 3185.00 

MAINTENANCE, DEPRECIATION AND LABOUR COST 
TOTAL COST = 12518.30 


PROPOSED METHOD 


UNUSED CAPACITY 


ROUTE 


, r ni i i 

- 14 13 


BUST ROUTE POLICY IS 3 

* i" , I? ' ' L j 

OPTIMAL NO. OF VEHICLES = ' 

TOTAL DISTANCE 108 


‘vV , 




: , 


■ 


? / f i 


OPTIMAL NO. OF VEHICLES = 5 

TOTAL DISTANCE* 130 
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T A B L E 5. 5 


TEST PROBLEM 4 



INPUT INFORMATIONS 


NO. OF CUSTOMERS = 18 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B(I» J) = I F I X ( A o RND Y5 ( Y 5 ) ) » 

A = 25.0 

DEMANDS FOR CUSTOMERS 

£ 0 60 70 40 35 3 2 48 23 32 41 35 58 35 31 37 42 

44 35 

CAPACITY OF VEHICLES = 175 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE 

1 

2 

3 

4 

5 


# 

UNUSED CAPACITY ROUTE PATH 


2 

1 

10 

4 

8 

9 

1 

1 34 

1 

1 

1 




6 

1 

13 

7 

18 

6 

1 

2 

1 

13 

17 

2 

12 

14 1 

3 

1 

16 

3 

5‘ 

19 

1 7 



OPTIMAL NO. OF VEHICLES 
TOTAL DI STANCE= 170 


COST OF TRANSPORTATION = 4165.00 

MAINTENANCE, DEPRECIATION AND LABOUR COST 
TOTAL COST = 13498.30 


PROPOSED METHOD 


ROUTE PATH 

1 I: 10 15 17 1 

1 8 19 1 

I 9 2 18 6 11 1 

1 12 7 5 3 1 

1 16 14 1 


UNUSED CAPACITY 

12 

22 


ROUTE 

1 


BEST ROUTE POLICY IS 2 

■ ■ 

OPTIMAL NO. OF VEHICLES 
TOTAL DISTANCE 152 


'■\ , f»T* . w / 

__ _ 


•V,;, 



TEST PROBLEM 5 


INPUT INFORMATIONS 


wO. OF CUSTOMERS = 20 

DISTANCE MATRIX GENERATED BY USING RNDY5 

= I FIX t A.RND Y5( Y 5) ) * 

' A = 27.5 

DEMANDS FOR CUSTOMERS 

25 18 16 23 30 12 20 18 20 15 20 15 17 23 

20 15 16 16 

CAPACITY OF VEHICLES = 100 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE UNUSED CAPACITY 

1 2 

2 2 

3 12 

i 4 A 


ROUTE PATH ** 

1 4 16 20 11 9 

1 12 21 2 17 1* 

1 13 8 5 6 1 

1 IS 15 3 10 19 


OPTIMAL NO. OF VEHICLES = 4 

TOTAL DISTANCE* 176 

COST OF TRANSPORTATION = 4312,00 

MAINTENANCE, DEPRECIATION AND LABOUR COST 
TOTAL COST* 11773.64 

PROPOSED METHOD 


7 

1 


= 7466.64 


ROUTE 

1 

2 

4| 

3 

4 


UNUSED CAPACITY 
2 

2 ■ 

it 

i 


: 

. . : ’ ' 'l i i ■■ ’ 

BEST ROUTE POLICY IS 2 




OPTIMAL NO. OF VEHICLES * 
TOTAL DISTANCE* 119 

I 

COST OF TRANSPORTATION = 
MAINTENANCE, DEPRECIATION AND 


ROUTE PATH 

1 1 11 9 17 6 1 

1 2 3 2 15 4 1 

1 1 10 21 19 12 1 

1 ? 16 8 5 18 1 

t 4 ,, A- $ 1 > 

: 1 : : f : . 

: | *. | « : T- 


1UTAL COST- K . 3 2.U 

...... ■ 


- . Vi :■ : y , f, '< ' ' i.T;. t ? ■ 

■ ■ ■■■• : 

29l ClnU COST. 74S6.6, 


1 * i „ i i "* 7 1* * F L P' . / T ? A ) ' 

’#> fKi r‘ ' J 
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TABLES,? 


TEST PROBLEM 6 

INPUT INFORMATIONS 


NO. OF CUSTOMERS = 22 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B( I, J) = I FIX ( A.RND Y5I Y5J I 

A = 15.0 

DEMANDS FOR CUSTOMERS 

13 14 16 14 8 14 12 15 14 10 14 8 14 12 12 15 

9 10 15 11 14 15 

CAPACITY OF VEHICLES = 50 

RESULTS 


CLARKE AND BRIGHT METHOD 


ROUTE 

1 

2 

3 

4 

5 

6 


UNUSED CAPACITY 

0 

ROUTE PATH 
1 2 11 

17 

15 

1 

7 

1 3 22 

9 

1 


2 

1 5 10 

13 

8 

1 

3 

1 12 19 

6 

23 

1 

1 

1 16 7 

18 

14 

1 

8 

1 20 4 

21 

1 



OPTIMAL NO. OF VEHICLES * 6 

TOTAL DI ST ANCE= 124 


COST OF TRANSPORTATION *» 3038.00 

MAINTENANCE, DEPRECIATION AND LABOUR COST= 11199.96 



TOTAL COST = 14237.96 


PROPOSED METHOD 


ROUTE PATH 

1 7 10 4 1 

1 13 8 14 5 

1 21 20 19 12 
1 ,3 11 11 6 
1 16 18 15 9 

1 2 23 22 1 


UNUSED CAPACITY 


ROUTE 


BEST ROUTE POLICY IS 




; ' 




IATION AND 
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TAB LfrS. 8 


TEST PROBLEM 7 


INPUT INFORMATIONS 


NO* OF CUSTOMERS = 25 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B(lfJ) = I F IXl A.RNDY5 { Y 5 ) ) 

A = 20 «0 

DEMANDS FDR .CUSTOMERS 

30 60 7 0 40 35 32 48 23 32 41 3 5 58 35 31 37 

4% 35 25 28 28 32 23 31 32 

CAPACITY OF VEHICLES » 200 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE 

1 

UNUSED CAPACITY 

19 

ROUTE PATH 
1 6 24 

7 

21 

19 

9 

2 

13 

1 8 16 

10 

4 

1 


■ 3 

18 

1 12 3 

20 

23 

2 

1 

4 

17 

1 14 17 

15 

25 

18 

1 

5 

1 

1 26 11 

5 

13 

22 

1 


OPTIMAL NO* OF VEHICLES = 5 

TOTAL DIST ANCE= 162 



PROPOSED METHOD 


ROUTE PATH 


ROUTE 


UNUSED CAPACITY 

41 


BEST ROUTE POLICY IS 2 

• , 

OPTIMAL NO. OF VEHICLES = 5 

TOTAL DI ST ANC£= 126 




: ' 


'Mtl.v; > ! 


COST OF TRANSPORTATION = 3969.00 

MAINTENANCE, DEPRECIATION AND LABOLR COST= 
TOTAL COST = 13302.30 
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TABLES* 9 


TEST PROBLEM 3 


INPUT INFORMATIONS 


NO* OF CUSTOMERS = 28 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B{ I, J) * I FIX I A.RNDY5C Y 5) I 

A = 20.0 

DEMANDS FOR CUSTOMERS 

13 14 16 14 8 14 12 15 14 10 14 8 14 12 12 15 

9 10 15 11 14 15 15 12 11 12 13 18 

CAPACITY OF VEHICLES = 55 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE 

UNUSED CAPACITY 

ROUTE PATH 




1 

14 

1 4 11 

24 

1 


2 

1 

1 5 29 

10 

13 

1 

3 

0 

1 12 3 

22 

2 

1 

4 

5 

1 16 26 

23 

15 

1 

5 

2 

1 18 20 

9 

7 

1 

6 

2 

1 25 21 

8 

19 

6 

7 

1 

1 27 1.4 

17 

28 

1 


t 



ROUTE PATH 

1 13 14 13 11 
1 21 24 IT 25 
1 V 9 20 6 

1 4 16 26 2 

1 22 23 15 27 
1 1 29 12 1 


UNUSED CAPACITY 


ROUTE 


BEST ROUTE 


TRANSPORTATION = 3479.00 

NCE, DEPRECIATION AND LABOUR COST*- 13066*62- 
ST* 16545.62 


■ 4 : 


OPTIMAL NO. OF VEHICLES 
TOTAL DISTANCE* 2C5 


COST OF TRANSPORTATION = 5022*50 

MAINTENANCE, DEPRECIATION AND LABOUR COST* 13066.62 
TOTAL COST* 18089.12 

PROPOSED METHOD 
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TABLE 5,10 


TEST PROBLEM 9 


INPUT INFORMATIONS 


NO. OF CUSTOMERS = 30 

DISTANCE MATRIX GENERATED BY USING RNDY5 
B( I, J) = I FIX( A.RNDY5( Y 5) ) 

A = 23.0 

DEMANDS FOR CUSTOMERS 

20 24 25 20 14 19 14 20 22 20 18 15 17 23 21 17 

16 19 15 18 20 19 24 16 18 17 15 25 16 19 

CAPACITY OF VEHICLES = 75 

RESULTS 


CLARKE AND WRIGHT METHOD 


ROUTE 

1 

UNUSED CAPACITY 

12 

ROUTE PATH 

1 4 21 2 

1 


2 

2 

1 

12 

9 5 

13 

1 

3 

7 

1 

14 

8 22 

27 

1 

4 

10 

1 

16 

3 11 

1 


5 

12 

1 

18 

17 28 

20 

1 

6 

4 

1 

19 

6 10 

25 

1 

7 .4 

3 

1 

23 

30 26 

7 

1 

8 

9 

1 

24 

15 31 

1 




OPTIMAL NO. OF VEHICLES 
TOTAL DISTANCE- 229 


COST OF TRANSPORTATION = 5610.50 

MAINTENANCE. DEPRECIATION AND LABOUR COST= 16799.94 
TOTAL COST = 22410.44 


PROPOSED METHOD 


ROUTE 

1 


UNUSED CAPACITY 


ROUTE PATH 

1 2 20 7 12 

1 19 22 28 17 




litl 


' 




.TAT I ON » 
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COMPARISON OF RESULTS - SUMMARY 
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CHAPTER VI 


VEHICLE DISPATCHING FOR KANPUR 3AKAKARI MILK BOARD - A CASE STUDY 


6 . 1 The proposed, model has been used for the development of distribution 
routes and vehicle dispatching strategies for Kanpur Sahakari Milk Board 
which is enraged in distribution of milk products in Kanpur City. The company 
has a wide distribution network and a fleet of delivery trucks. A brief 
description of the operations of hempur Rabakari Milk Board (KSMP). is 
given in the following paragraphs * 

6.1.1 Procurement Raw milk is collected at various collection centres 
leaving milk potentiality . At some of these collection centres where the 
milk potentiality Is high, company has got its chilling centres for 
preserving the collected mile,. Presently three company owned tankers and 
six hired, private vehicles transport this collected milk to the Dairy for 
processing . 

6.1.2 Pr ocessing : Procured raw milk and imported milk powder is 
processed to yield pasteurized mile products namely Milk, butter, ghee and. 
cream etc. Processing is done for preparing termed, double termed and 
skimmed milk etc . having different fat percentages . Processed milk Is 
bottled in half-litre capacity bottles and kept in cold storage, ready for 
distribution . 

6.1.3 Distr ib ution : The Company has a wide distribution network and a fleet; 
of delivery trucks. The milk products are distributed to the general public ! 
through the agents at various delivery centres. In other words. Distribution j 
Chqnnel followed is Producer to retailer to consumer. Company is distributing i 
the products twice a day (morning & evening) to about 250 agents in Kanpur 

'.A .A. A e- . ,;A e/APe, ' • Avv’A v ,P : ; lA'Ap, .A,; AA ; '■ /AAA^AAAA/ee/ p A’ : e A A '-ALA ; A; A; .'.A; : py, ■ A ; 

city who sell them to public. Therefore, so far the company's distribution 
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system is concerned, agents are the customers. Main orcducts to be delivered 
are milk and butter. For Ghee, there are only two delivery centres, one 
at the dairy itself and other at Super Market, Kanpur. There are just 
two or three deliveries rer month to Super Market , therefore the 
distribution of Ghee has not been scheduled in the regular dis tribution 
system. Presently, the company is engaging Nine delivery trucks for 
distribution, seven out of them are bigger (Capacity of 4000 bottles) and 
two smaller (capacity of 1600 bottles) . All the nine trucks are allocated 
to 9 routes in the morning shift and six trucks (5 bigger and 1 s mall er) 
trucks for 6 routes in evening shift, demand in evening being lesser. 

Delivery trucks are loaded with crates (each containing 2^ bottles) 
at the loading centre. The loaded truck is accompanied by two helpers, 
one route clerk and obviously one driver for the distribution of products. 
Helpers load and unload the crates at the delivery points and depot. 

Route clerk is responsible for the proper distribution and money transaction 
etc. Full delivery to each customer is being made simultaneously. At 
each customer point , the crates of returned empty bottles are loaded 
in truck. After delivering to every customer on the routes, trucks 
return to the dairy ' (depot) and all the crates of empty bottles in truck 
are unloaded. This completes the distribution cycle. 

For conducting the case study, various data need to be collected 
as an input information to the model. Following section discusses about 
the data collection for the use of present study. 


■■6.2 Data Collection ; The strategy used for the collection of data 

has been discussed in this section. Mainly five set of data are req uir ed 
viz.: 
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Before going ahead, it is opportune to disci^ss seme procedural 
simplifications made, for making the study convenient, as follows: 

^•2.1 Aggregation of delivery points 

The number of delivery points at which milk products are 
distributed in Kanpur City is around two hundred and fifty. Because of 
computational difficulties and limited computer memory, it is not possible 
to consider them individually. Therefore, delivery points have been grouped. 

A group of customers are supposed to be supplied from a single delivery 
point. Firstly the delivery points in a locality have been grouped 
and each group has been represented by a delivery point in that locality. 

These delivery points have been plotted on the map of Kanpur City. These 
delivery points turned out to be 88 in number. Since the proposed model 
can handle- about 30 customers because of computer storage limitations on 
IBM/7044, further grouping of these delivery points was needed. The 
plotted delivery points forming the clusters have been grouped first. The 
delivery points which do not form the clusters have been grouped, based on 
proximity criterion. The delivery points within a diameter of 1.5 km 
have been considered to be at a single delivery point . By this proximity 
criterion twenty three groups have been formed. But it was observed 
that the number of individual delivery points in some groups, which were 
in dcnsly populated area, was quite high as compared to that in thinly 
populated area. Therefore, second criterion was used for grouping. 

This criterion groups the indiviaual points in densly populated area within 
diameter of 1 km which yields the number of delivery points in each 
group more or less same 1 . e . , between 8 to 10 delivery points, in each 
group. Thus, .in nutshell, the criteria used for aggregation of 
delivery points can be spelled out as 'the delivery points within a diameter 
of 1.5 km in the densely populated area and delivery points within 1.0 km 
in a thinly populated area are considered to be at a single delivery 
point which has been named as zonal delivery point/ Above procedure of 
aggregation of delivery points yields Twenty seven zonal delivery points which 
have been considered as the nodes in the distribution network 



for the present case study. These nodes have been listed in table 6.1. 

The total requirement of individual delivery points group* named as 
zone is represented as the load at zonal delivery point of that zone. 

The total delivery to the delivery in any zone is supposed to be 
delivered at the zonal delivery point. 

6.2.2 Milk as the only product to be delivered : In the present 
distribution system, among the milk products which are distributed, milk 
has been considered as the main product for distribution. The ghee and 
butter are needed in smaller quantities and are being distributed at 
fewer delivery centres. ' Moreover, ; the capacity contribution of 
these products to the vehicles for distribution being very small, these 
products have not been considered for distribution. This assumption 
leave milk as the only product for distribution in our study. 

6.3 The strategies used, for the collection of data are outlined in 

the following paragraphs: 

6.3.1 Distances between the nodes in distribution network ’ As discussed 
in Chapter 3, the shortest distances between the nodes are required as 
an input data to the model. Shortest distances are obtained from 
direct linked distances as discussed in Chapter 4 , by shortest path 
algorithm. Therefore, data have been collected to prepare direct 
linked distance matrix showing the distances of direct path between 
every pair of nodes. For collecting these data, the zonal delivery 
points and milk depot are plotted oh the road map of the Kanpur city. 

The road distances between each pair of zonal delivery points, which are 
connected by road are determined. If there are many possible direct 
links between two zonal delivery points, the shortest direct link is 
considered. The unit of distance matrix is one fcrrth of a kilometer . 
Wherever there is no direct link between two delivery points distance 
between them has been assumed as infinity. Distance matrix has been sham 
in liable 6. 2. (a) . Shortest distance matrix,’ as computed by 
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NODE 

TABLE 6.1 

ZONE-NAME 

1 

Milk Board 

2 

Khapra Mohal 

3 

Generalganj 

4 

Ram Narain Bazar 

5 

Chowk 

6 

Birhana Road 

7 

Cantonment 

8 

Latouch Road 

9 

Chamanganj 

10 

Anwar ganj Station 

11 

Faithfulganj 

12 

Govindnagar (East ) 

13 

Govindnagar (West ) 

14 

Ordinance factory 

15 

Shastrinagar 

16 

Lajpat Nagar 

17 

L.L.R. Hospital 

18 

IIT 

19 

Aryanagar 

20 

Nawabganj 

21 

Azadnagar 

22 

Gwaltoli (South) 

23 

Gwaltoli (North) 
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NODE 

ZONE NAME 

24 

Fazalganj 

25 

Jawaharnagar 

26 

Kidwa inagar 

27 

Civil Lines 

28 

Parade 


1.1. T. KAN FUR 

rFNTR AL LIBRARY 


Acc. No. 
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TABLE 6.2 (M 

SHORT. ST DISTANCE MT. RIX 



1 

■ 2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 . 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

I 

0 

37 

32 

3 3 

31 

36 

36 

27 

22 

18 

32 

6 

6 

24 

22 

17 

25 

65 

26 

37 

44 

30 

35 

10 

18 

14 

30 

27 

2 

37 

0 

6 

7 

10 

6 

8 

1G 

20 

19 

8 

37 

37 

47 

41 

34 

38 

73 

25 

36 

43 

20 

21 

30 

22 

24 

16 

13 

3 ■ 

32 

8 

0 

4 


5 

16 

5 

12 

14 

12 

32 

32 

42 

35 

28 

25 

65 

17 

2 6 

35 

12 

13 

22 

14 

25 

12 

5 

4 

33 

7 

% 

0 

B 

4 

15 

8 

13 

17 

13 

35 

35 

45 

36 

29 

26 

66 

18 

29 

36 

13 

14 

23 

15 

28 

11 

6 

ft 

31 

10 

4 

3 

0 

7 

18 

? 

11 

16 

10 

33 

34 

44 

34 

27 

24 

64 

16 

27 

34 

11 

12 

21 

13 

27 

8 

4 

h 

36 

6 

\ K 

4 

7 

0 

14 

9 

16 

18 

9 

36 

36 

46 

39 

32 

29 

69 

21 

32 

39 

16 

17 

26 

18 

27 

12 

9 

i 

36' 

8 

16 

15 

18 

14 

0 

18 

28 

26 

12 

37 

42 

54 

48 

41 

40 

80 

33 

44 

51 

28 

29 

3 b 

30 

22 

24 

21 

* b 

27 

ID 

5 

8 

7 

9 

18 

0 

10 

9 

17 

27 

27 

37 

31 

24 

23 

63 

22 

33 

40 

17 

18 

23 

19 

20 

15 

10 

9 

22 

20 

12 

13 

11 

16 

28 

1C 

0 

10 

21 

22 

24 

38 

32 

25 

22 

62 

14 

25 

32 

14 

15 

18 

12 

25 

14 

7 

10 

18 

19 

14 

17 

16 

ia 

26 

9 

10 

0 

14 

18 

18 

28 

22 

15 

14 

54 

18 

26 

33 

24 

25 

16 

12 

20 

24 

17 

11 

32 

3 

12 

13 

10 

9 

12 

17 

21 

14 

0 

32 

32 

42 

36 

29 

28 

68 

26 

37 

44 

21 

22 

50 

23 

18 

18 

14 

12 

6 

37 

32 

35 

33 

36 

3? 

27 

22 

18 

32 

0 

10 

28 

18 

18 

26 

66 

32 

38 

45. 

36 

3? 

16 

24 

15 

36 

29 

13 

6 

37 

32 

35 

34 

36 

42 

27 

24 

18 

32 

10' 

0 

18 

16 

20 

28 

68 

32 

40 

47 

36 

39 

16 

24 

20 

36 

31 

14 

24 

47 

42 

45 

44 

46 

54 

37 

’38 

28 

42 

28 

18 

0 

19 

26 

32 

72 

40 

44 

51 

48 

53 

28 

36 

38 

48 

45 

15 

22 

41 

35 

36 

34 

39 

48 

31 

32 

22 

36 

18 

16 

19 

0 

7 

13 

53 

21 

25 

32 

31 

38 

14 

21 

33 

33 

30 

16 

17 

34 

28 

29 

27 

32 

41 

24 

25 

15 

29 

18 

20 

26 

7 

0 

6 

48 

16 

20 

2? 

26 

31 

7 

14 

31 

26 

23 

17 

■25 ' 

33 

25 

26 

24 

29 

40 

23 

22 

14 

28 

26 

28 

32 

13 

8 

0 

40 

, 8 

12 

19 

18 

25 

15 

11 

34 

23 

20 

18 

65 

73 

65 

66 

64 

69 

80 

6 3 

62 

54 

68 

66 

68 

72 

53 

48 

40 

0 

48 

52 

59 

58 ! 65 

55 

51 

74 

63 

60 

19 

'26 

25 

17 

18 

16 

21 

33 

22 

14 

18 

26 

32 

32 

40 

21 

16 

8 

48 

. . 0 

11 

18 

10 

17 

16 

8 

38 

18 

12 

20 

37 

35 

28 

29 

27 

32 

44 

’33 

25 

26 

37 

38 

40 

44 

25 

20 

12 

52 

11 

0 

* 7 

21 

28 

27 

19 

46 

29 

23 

21 


43 

35 

36 

34 

39 

51 

40 

32 

33 

44 

45 

47 

51 

32 

27 

19 

59 18 

I 

0 

28 

35 

34 

26 

53 

36 

30 

22 

30 

20 

12 

13 

11 

* 16 

28 

17 

14 

24 

21 

36 

36 

48 

ii 

26 

18 

58 

10 

21 

28 

0 

7 

20 

12 

< 37 

10 

7 

23 

35 

21 

13 

14 

12 

17 

29 

18 

15 

25 
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TABLE 6.2 (a) 

DISTANCE MATRIX „ .* 



1 

2 

3 

4 

5 

6 

? 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

2D 

21 

22 

23 

24 

25 

26 

27 

28 

1 

Q 

999 

999 

999 

999 

999 

999 

999 

22 

18 

999 

6 

6 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

10 

999 

14 

999 

999 

2 

999 

3 

8 

7 

10 

6 

8 

10 

999 

999 

8 

9 99 

999 

999 

9^9 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

24 

16 

13 

3 ' 

999 

3 

Q 

4 

4 

5 

999 

5 

999 

999 
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999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 

999 
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NOTE- D ( I» J ) MEANS SHORTEST DIRECT LINK BETWEEN I AND J NODE 

DU, J) =999 MEANS NO DIRECT LINK BETWEEN I AND J NODE 
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6*3*2 Demand at each delivery point : ' Demands at zonal delivery points 
are the input data to the model. Demand at any zonal delivery point is the 
summation of the demands of all delivery points grouped in that zone. 

Therefore the demand at each delivery point is found first. Demand at 
any delivery point is the quantity indented by the agent at that delivery 
point . Hiis quantity indented by the agent can be obtained frcm the 
indent book available with Dales Office of the Milk Board. Company is 
distributing milk twice- a- day (both in the morning and in the evening) . 
Therefore , the morning and evening demand data at each delivery point 
have been collected first. With these data, the average morning and 
evening demand at each delivery point is calculated for each month. Indent 
book was available from January, 1974 only. So, because of non-availability 
of enough demand data, no forecasting model could be tried for estimation 
of demands. However it was observed that morning and evening demand 
at each delivery point were not significantly varying for January, 

February and March months . Therefore , the average morning and evening 
demand for these three months were found. These average demands at 
delivery points have been called as Winter morning and evening demands. 
Similarly, the average morning and evening demands in months April, May and 
June haave been found- and have been called Summer morning and evening 
demands at delivery points. These average demands of all the delivery 
points in each zone are summed up to give average morning and evening demand 
(for sunnier and winter separately) at that zonal delivery point. Table 
6,3 indicates the average demands at zonal delivery points. 

6.3-3 Operating Characteristics of Vehicles 

As was mentioned in Chapter 3 S all vehicles are in same condition. 

So operating, characteristics depend upon the size of vehicle and not on the 
condition of the vehicle. Fcr each vehicle, the average travelling speed and 
maximum allowable load i.e. capacity is noted. These data have been tabulatec 
in Table 6.5. 



TABLE 6.3 
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DEMAND DATA (KANPUR SAHAKARI MILK BOARD ) 



SJMMLR 


WINTER 



M0RN|.ifG 


EVENING 


EVENING 


MORNING 


720 
5 50 
750 
7 70 
9 GO 
1100 
10 %0 
1150 
5.00 
740 
920 
670 
B20 
600 
7 60. 
1A-00 
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TABLE 6.4 


VEHICLE CAPACITIES 


Number of available vehicles - 9 

Vehicl e number Ca pacity 

(bottles) * 

1 4000 

2 4000 

3 1600 

4 4000 

5 1600 

6 4000 

7 4000 

8 4000 

9 4000 0 


*Each bottle of half litre capacity. 



TABLE 6.5 


OPERATING CHARACTERISTICS OF VEHICLES AM) DISTRIBUTION COSTS INVOLVED 


I) Average stop-off- -time at each individual agent = 3-5 minutes 

ii) Average stop-off-time at depot = 120 minutes 

iii) Maximum Permissible route time = 300 minutes 

iv) Maximum Permissible Vehicle time = 420 minutes 

v) Average fuel consumption of diesel for 

vehicles with capacity of 4000 bottles = 4 km. per litre 

vi) Average fuel consumption of petrol for for small - vehicle s = 5 km/litre 

vii) Average speed of bigger vehicles (capacity 4000 bottles) in 
morning = 20 k.m. : per hour 

viii) Average speed of bigger vehicles in evening 
= 15 km. per hour 

ix) Average speed of smaller vehicles (capacity 1600 bottles) in 
morning = 25 km. per hour 

x) Average speed of smaller vehicles in evening = 20 km. per hour 

xi) Maintenance cost per vehicle per shift = Rs. 167 per month 

xii) Depreciation cost* for bigger vehicle = Ps. 277 per month 

xiii) Depreciation cost* for smaller vehicle = Rs. 84 per month 

xiv) Salaries 

Drivers salary = Rs. 250/~ r.nu 

Route clerk's salary = Rs. 275/- p.m. . 

Helper's Salary = Rs. 210/- p.m. 

Total Labour cost = Rs. 045/*- p.m. 

xv) Cost of petrol = Rs. 3*27 per litre 
Ccst of diesel = Rs. 1.0c per litre 


* An nua l depreciation cost is 10$ of initial purchase cost. 



6.3-^ Distribution costs involved- Average consumption of fuel for 
each vehicle has been found. Knowing it, cost of fuel per kilometer for 
each vehicle has been found, labour Cost includes Salary of driver, 
helpers and route clerk. These data have been taken from the accounts 
office of the milk board. Maintenance and depreciation costs for each 
vehicle on monthly average basis have been taken for each vehicle from 
the transport section of the milk board. Table 6.5 shows the different 
costs involved in distribution. 

6.3.5 Stop off - time at milk de pot .and delivery points : Estimates of stop' 
off- time at the depot and the delivery points have been obtained from the 
dispatching section, Kanpur Sahakari Milk Board. St on-off -time at a 
zonal delivery point, which is required as input data to Model, has been 
taken as summation of average stop-off-time at all the delivery points 
in that zone. Thble 6.6 shows the stop-off -times at the zonal delivery 
points . 


* average stop- -off- time 
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CHAPTER-VII 


RESULTS & DISCUSSION 

In this chapter, the results obtained for the distribution system 
of the Kanpur Sahakari Milk Board have been presented and discussed. The 
distribution policies obtained using the model have been compared with the 
distribution policies presently adopted by the management of the Kanpur 
Sahakari Milk Board. A sensitivity analysis has been conducted for 
observing the sensitivity of results to the various parameters. Recomm- 
endations for the further work are given at the end of this Chapter. 

7.1 Results for the distribution system of the Kanpur Sahakari Milk 

Board . The proposed distribution routes and vehicle dispatching strategies, 
for the Kanpur Sahakari Milk Board presented in this section. Ihe best 
possible results obtained by the proposed model, for the morning and evening 
shifts of the sunnier and winter months, are represented in tables 7.1 and 
7.2, respectively .Table 7-1 indicates that the proposed vehicles, for the 
distribution in morning shift of summer months, are the vehicle numbers 
1,2, 3.-4. 6,7.8*. This table also indicates the route path for each 

route and the vehicle allocation to the routes. Untilized capacities of 
vehicles (UVC) are given in column 4. The time (RT) for which vehicles 
are on the routes and the time for which each veliicle is engaged for 
distribution is indicated in column 5 • The cost of transportation and the 
total cost, corresponding to the proposed routes and vehicle dispatching 
strategies for morning shift of surrmer months, are Rs. 1,239.77. per month 
and Rs. 10,768.77 per month respectively, as indicated by the table 7 -l.(a) 

The results for the morning and evening shifts of winter months 
are tabulated in table 1.2 which indicates the route path for the routes 
obtained and allocation of the selected vehicles to the routes. Costs 
of transportation for morning and evening shift are Rs. 1,407.10 per 
month and Rs. 1.168/- per month respectively. Total distribution costs 
are Rs. 9 547.10 per month and Rs. 7,725/- per month respectively, as 
indicated by table 7 .2 .(a) and 1.2 (b). 

#The available vehicles have be^n numbered as indicated in table 6.4 


TABLE 7.1(a) 

EOF SpyiEH-ffiRNrrNG DEMANDS 
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Route 

Vehicle Allocation 

Number Capacity 
(bottles) 

Unused 

capacity 

(UVC) 

Route 

time 

(RT) 

Mints. 

Vehicle 

Time 

(VT) 

Mints. 

Route Path 

1 

3 

1600 

10 

76 

196 

1-13-12-1 

2 

1 

4000 

270 

188 

308 

l_3„2-7-26-l 

3 

2 

4000 

100 

185 

305 

1- 27-23-22-28-1 

4 

4 

4000 

160 

255 

375 

1-18-21-20-19-1 

5 

6 

4000 

CO 

o 

18.1 

301 

1-25-9-10-2 4-1 

6 

7 

4000 

160 

216 

336 

1-11-6-4-5-3-1 

7 

8 

4000 

420 

185 

305 

1-14-15-17-16-1 


Cost of transportation = Rs. 1 ,,239 -77 per month 
Total cost of distribution = Rs, 10 2 768. 77 per month 
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TABLE 7.1(b) 


FOR SUMMERrEVEraftG DEMANDS 


Route 

Vehicle-< 

Number 

allocation 

Capacity 

(bottles) 

Unused 

capacity 

(IJVC) 

Route 

time 

(RT) 

Mints. 

Vehicle 
time (RT) 
Mints. 

Route path 

1 

3 

1600 

140 

69 

189 

1-25-24-1 

2 

1 

4000 

120 

259 

379 

l_19-,20-21-22- 

11-26-1 

3 

2 

4000 

75 

293 

413 

1- 16-17- -18-15- 
-14-13-1 

4 

4 

4000 

335 

233 

353 

1-7-2-6-4-3- 
8- 10-1 

5 

6 

4000 

80 

215 

335 

1-27-23“ 5-28- 
9-12-1 


■ ■ sssslfessiist isiiiiii ■ : si 11 iiiii 1 iiiiiiiiliiil is stf® 

Cost of transportation = Rs« 1 } 107.00 per month. 
Total cost of distribution = Rs. 7. 858 . 00 per month. 



TABLE 7.2(a) 


FOR WINTER-MORNING DEMANDS 


Route Vehicle-allocation 

Number Capacity 

(bottles) 

Unused 

Capacity 

(UVC) 

Route time 
. (RT) 

Mints , 

Vehicle 

Time 

(VT) 

Mints. 

Route Rath 

1 1 

4000 

200 

258 

378 

1-4-6-2-7-11- 

10-1 

2 2 

4000 

140 

228 

348 

1- 26-3-5-23-27- 

3 4 

4000 

250 

180 

300 

1-13-9-19-25-1 

4 6 

4000 

190 

283 

403 

1- 21 - 20 - 18 - 
17-16-1 

5 3 

lg00 

140 

105 

225 

1-28-22-1 

6 7 

4000 

150 

245 

365 

1 24-15-14-8 

12-1 


Cost of transportation = Rs. 1,407.10 per month 
Total cost of distribution* Rs. 9,547.10 per month 



TABLE 7.2(b) 


for IftJINTER- EVENING DEMANDS 



1 

3 

1600 

60 

207 

327 

1-18-20-21-16-1 

2 

5 

1600 

50 

94 

214 

1-17-24-1 

3 

1 

4000 

330 

263 

383 

1-9-25-19-15-14 

13-12-1 

4 

2 

4000 

O 

CO 

251 

371 

1-28-5-2-6- 4- 3-1 

10-1 

5 

4 

4000 

700 

245 

365 

1-26-7- 11-27-23 


22-1 


Cost cf transportation = Rs. 1,168.00 per month 

Total, c^st of distribution = Rs. 7 ,725.00 per month 


7.2 Comparison of the total distribution cost by proposed 

dispatching policies with that of the currently used policies of Kanpur 
Sahakari Milk Board: The proposed distribution routes and the vehicle 
dispatching strategies for the Kanpur Sahakari Milk Board have been 
discussed in Section 7.1 and are tabulated In Tables 7.1 and 7.2. .The 
current distribution costs for the existing distribution system of 
Kanpur Sahakari Milk Board s are given in table 7*3. Table 7 .4 
shows the comparison of the distribution cost by the proposed vehicle 
policies with that of the currently used policies of Kanpur Sahakari 
Milk Board. The savings in total distribution cost obtained by the use 
of proposed model, are indicated in last column of Table 7.4. F° r the 
monring shift of simmer months, the saving of 24$ has been observed. For 
the evening shift of sunmer months, the saving has been observed to the 
tune of 21$ . 


For the morning and evening shifts of winter months, the savings 
in total distribution costs have been observed more than that. for sunmer 
months . The reason for the higher savings in winter months can be 
explained as follows: The number of vehicles (obtained by proposed 
model) employed for satisfying customers ’ demands in winter months , are less 
than that employed in summer months, because of less demand in winter 
months. On the other hand, the number of vehicles enployed far the 
distribution, by the management, Kanpur Sahakari Milk Board, are same 
for both winter and summer months. For the morning shift of winter months, 
the saving obtained by proposed solution is 32.6$. For the evening shaft 
of winter months, the saving observed is 22.8$. In general, the gavihgs 
obtained by the proposed distribution policies over the currently used, 
policies by Kanpur Sahakari Milk Board are to the tune of 20 to 32$. 

7 p Fffoct of vari ous vehlcle - ^°t-.?;o1 ection policies and vehicl e 

allccatioTinlies.' In this section, the effect of various vehicle- set • 
selection policies and vehicle -allocation policies on the objective function 
has been discussed. To study the effect of various policies on 
the objective function, the problem has been solved with various policies 


TABLE 7.3 


DISTRIBUTION COSTS FOR EXISTING DISTRIBUTION SYSTEM OF KANPUR SAHAKARI MILK BOARD 
(Estimates given by Management * Kanpur Sahakari Milk Board) 

I) For morning shift 

Cost of transportation = Rs. 2,100.00 per month 

Maintenance , labour and depreciation cost = Rs. 12,113.00 per month 

Total cost = Rs'. 14,213.00 per month. 

ii) For evening shift 

Cost of transportation = Rs. 1,797.00 per month. 

Maintenance, labour and depreciation cost = Rs. 8,140.00 per month 
Total cost of distribution = Rs. 9,937.00 per month 


TABLE 7.4 

COMPARISON OF THE TOTAL DISTRIBUTION COST BY PROPOSED DISPATCHING 
POLICIES WITH THAT OF THE, CURRENTLY USED POLICIES OF KANPUR 
SAHAKARI MILK BOARD 


Total cost of distribution Savings % 

(Rs /month) 

Existing System Proposed System 


Summer-morning 

14213.00 

10768.77 

24.05? 

Summer- evening 

9937-70 

7858.00 

21 , 2 % 

Winter-morninm 

14213.00 

9547.10 

32 . 6 % 

Winter-evening 

9937.70 

m’ :/ ' :c 

7725.00 

. 

22 . 8 % 



for the different sets of demands . Different sets of demands have been 
obtained by variating the original demands at each customer from 

* 20$ to +20% in step of 5$ . For each demand set , the total costs of 
distribution for different vehicle set-selection-policies and vehicle- 
allocation policies ahve been indicated in Table 7,5. These results indicate 
that different combinations of vehicle-set-sel ection policies and vehicle- 
allocation policies yield different costs of distribution. 

’7*3<1 Ef fect of Vehicle' set-selection policies 

Vehielo-set-selection policy selects the set of vehicles for 
distribution to satisfy the customers 1 demands . In this section s the 
effect of different vehicle-set-selection policies on the cost of dis- 
tribution has been discussed. Fig, 1 shows the maintenance , depreciation 
and labour costs for different vehicle-set-selection policies with 
variation in sunmer-moming demands of each eustomer from - 20 % to 
+20 f in step of 5%. Fig, 1 shows that vehicle-set-selection policy 2, 
in general 3 gives higher cost, exceptional cases being only when demand 
at each node increases by 5% or reduces by 10$. Use of vehicle-fee- 
selection policy 2 results in higher cost because in this policy 3 the 
available vehicles are arranged in ascending order of their capacities and 
smaller capacity vehicles are selected first. Therefore for satisfying 
the customer demands s vehicle-set- selection plicy2 requires more number of vehicles 
than required by vehicle-set-selection policy 1 in which vehicles are 
•arranged in decending orcier of the capacities, then demand increases by 5$ s 
the reason for lower cost by using policy 2 than that by using policy 3-. is 
that policy 2 selects 2 Small end 6 big capacity vehicles whereas the 
policy 3 selects 1 small and 7 big capacity vehicles. The number of vehicles 
in both cases are same but higher depreciation cost s in case of big capacity 

• 3 to yield, higher cost. vpppi'.pg'p. bit 

In Fig, Z a graphical relationship between the total cost of 
distribution and the various demand sets considering, different vebicle- 
set-selection policies is presented. Fig. 2 shows that in all the cases, 
vehicle -set- selection policy 2 gives higher cost. When demand is varied by 



OP DISTRIBUTION (Rs , /MONTH) FOR DIFFERENT SETS OF DEMANDS, USING VARIOIB 
» SELECTION POLICIES AND VEHICLE ALLOCATION POLICIES 
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TABLE 7-6 
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Cost of maintenance;, depreciation and labour cost (Rs. /Month) , for 
different sets of demands, using various set policies. (Different 
sets of demands have beeb obtained by variating the summer-morning- 
demands in Step of + 5 %) 


Percentage 

Variation 

Total 

demands 

(bottles/ 

day) 

Cost of maintenance, depreciation & labour cost 
Vehicle- . Vehicle- Vehicle- • (Rs, /month) 

set- set- set- 

selection selection selection 

Policty 1 Policy 2 Policy 3 

-20 

* 19440 

8334 

9335 

8l40 

-15 • 

20655 

8334 

9335 

8l40 

-10 

21&70 

8334 

9335 

9529 

r 5 

23085 

8334 

10724 

9529 

0 

24300 

9723 

10724 

9529 

■5 

'25515 

9723 

10724 

10918 

10 

26730 

9723 

12113 

10918 

15 

27945 

10918 

12113 

10918 

20 

29100 

10918 

12113 

10918 



Maintenance, Labour & Depreciation costs (Rs./mont 


11000 


10000 


9000' 


8000 


V ehic l e - set - se lection policies 

Policy 1 • 

Policy 2 
Policy 3 o 
Policies 1 & 3 © 


1 


± 


-20 -15 -10 -5 0 5 10 15 20 

Percentage variation in summer -morning demand set 

Fig. 1 - Effect of various Vehicle -Set -Selection Policies on the 
Maintenance, Labour and Depreciation cost, for the different 


sets of Demands. 


Total Distribution Cost ( Rs./ month) 


13000 


12000 


11000 


10000 


9000 

»20 -15 -10 -5 0 5 10 15 20 

Percentage Variation in Summer-morning Demand Set 

Fig. 2 - Effect of Various Vehicle -Set-Selection Policies on the 
Total Distribution cost, for the different sets of Demands. 


Vehicle - Set -Selection Policies 

Policy 1 0 

Policy 2 □ 

Policy 3 ■ A 

Policies 1 & 3 <§> 
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o 
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O 


O 

A 


O 


o 
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5% s though the maintenance , labour and depreciation costs are 

observed to be lower In case of policy 2 than in case of oolicy 3 

(Fig. 1) _ ; . the total distribution cost is observed to be higher 

with use of policy 2 (as indicated by Pig. 2) This is due to the fact 

that the vehicle- set-selection policy 2 selects more number of small 

vehicles which use petrol (big vehicles use diesel) . Petrol being 

more expensive than diesel., the cost of transportation is more with 

use of vehicle- set- selection policy 2. Therefore, one can conclude 

that the vehicle-set-selecticn policy 2 for selecting the vehicles is 

uneconomical . However , as the results show 7 (Pig. 2) there is nothing 

much to choose between policies 1 & 3° On the other hand, policy 3 being that 

of random selection, it is recommended that policy 1 be used for 

inp lament at ion . 

7«3°2 Effect of Vehicle-allocation polic3.es 

Vehicle-allocation policy allocates the vehicles (selected by 

vehicle-set-selection policies) to the initial routes and then routes are 

adjusted, accordingly, for the capacity of vehicle allocated to it. 

Figs. 3 through 5 show the effect of various vehicle allocation policies 

on the total distribution cost with variation in demand for the particular 

vehicle, set-selection oolicy . When the vehicles selected by a vehicle-set--. 

selection policy, are of the same capacities, every vehicle -allocation policy 

gives the same results as indicated by Fig. 3* This is because the 

vehicle allocation to the initial routes is the same in every case. 

Howevef , the various policies yield different results when vehicles of 

different capacities are employed. No conclusion could be drawn about 

the superiority of one vehicle allocation policy over the other. Results 

show (Fig. 3) that there is nothing much to choose between vehicle™ 

c 

allocation policies 1, 2 and 3= Fig. - shows the relationship between 
cost of transportation and demand with use of various vehicle allocation 
policies . 



Total Distribution Cost (Rs 


Ve hide-allocation, policies 

Policy 1 

* 

Policy 2 

A 

Policy 3 

m 

Policies 1,2, 3. 

# 





Total distribution Cost (Rs. / month) 






tal Distribution Cost (Rs. / month ) 
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12000 


Vehicle -allocation policies 


Policy 2 a 
Policy 3 □ 

For Vehicle -Set-Selection 
Policy 3 


11000 


10000 


9000 


■20 -15 -10 -5 0 5 10 15 20 

Percentage Variation in Summer- morning Demand Set 


Fig. 5 -Relationship between the Demand and the Iota! Distribution 

Cost with use of different Vehicle -Allocation Policies and 
Considering Vehicle- Set -Sc lection Policy 2. 
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Percentage Variation in Summer-morning Demand Set 

Fig. 6 -Relationship between the Demand and the Cost of Trans- 
portation with use of various Vehicle-Allocation Policies 
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7.4 Sensitiv&y Analysis 

It consists of changing the values cf parameters either one 
at a time or simultaneously. Her e-in , changes have been made only one 
at a time. Among the various parameters , the parameter subjected to 
sensitive analysis, in present work* is 'demand! parameter for 
ob serving its effect on the objective function. This is due to the 
fact that demand varies frequently, depending upon the customers' 
requirement. Table 7-7 and Fig. 7 show the effect of variation In Winter- 
Morning: demand set cn the optimal cost of transportation and optimal 
total distribution cost obtained by the proposed method. It is 
observed that the total cost of distribution for demand variation of 
- 20 % to - 6 % does not change significantly s maximum being Rs. 8173/- per 
month and minimum being Rs. 79^3/-- per month. Similarly, no 
significant variation in total distribution cost for the present demand 
set, is observed till the demand at each customer increases by 12% 3 
maximum being Rs. 9547/" p.m. and minimum. being Rs. 949V-- p.m. . 

Because of practical difficulties involved in frequent change of the vehicle 
dispatching strategies with the change in customers' demand, it is 
recommended that upto 12 % increase in demand at each customer, same 
vehicle routes and vehicle dispatching strategies are followed, because 
cf no significant change in the optimal cost of distribution. Because 
of capacity constraints of vehicles, the routes for distribution in practice 
should be those obtained for maximum demand at each customer i.e., 
for the demands increased by 1255 at each customer. Table 7.8 shows the 
routes and vehicles allocated to them which can be used for satisfying 
the customers" demands .anywhere between present demands and demands increased 
by 12 % at each customer. 

Similarly, Fig. 7 shows that the total cost of distribution, for 
demand variation between -20? and -6?, does not change significantly, 
maximum being Rs. B173/- per month and minimum being Rs. 7983/- pen month. 
Table 7.9 shows the optimal routes and vehicle dispatching strategies 



TABIE 7.7 
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EFFECT OF VARIATION IN TOWER-KPIMG DEMAND SET ON 
OPTIMAL COST OF TRMSPORTCTTON -AND TOTAL COST OP 
DISTRIBUTION 


Percentage 

variation 

Cost of Transportation 
./&e©tto- 

Total cost of distribution 
' (Rs ./month) 

- 205 ? 

1422.70 

8173.7 0 

-15* 

1095.50 

3040.5° 

O 

H 

I 

1038.80 

7983.80 

- 7 % 

1113.75 

0058.75 

~ 6 % 

1103.60 

8048 . 60 

- 5 % 

1303.70 

0443.70 

05 ? 

1407.00 

9547.00 

5 % 

1162.35 

9496.35 

10 % 

1184.62 

9518.62 

12 % 

1150.20 

9494.20 

15 % 

1246.20 

10775.2Q 

20 % 

1368.9.0. 

IO897.90 
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Recommended routes and vehi cle-di spat ching strategies for the variation 
of demands at each customer between -20$ and - 6 % of present winter-morning 


demand set. 



Route 

Vehiclfe-allocat ion 
Number Capacity 

Route path 

1 

1 

4000 

1 n 5- 16-18-21-20-19-1 

2 

2 

4000 

1-12-0-17 24-1 

3 , 

4 

4000 

1-10-9-22-23-27-1 

4 

6 

4000 

1-26- 7-2-4-5- 28-1 

5 

7 

4000 

1-25-3-6-11-14-13 1 


Total cost of distribution = Rs « 

804-0 . 6 ©er month 


Table 7-9 

Recommended routes and vehicle- dispatching strategies fcr the variation of 
demands at each customer upto +12$ of present winter-morning demand set 


Route Vehicle-allocation Route path 

Number Capacity 


1 

*1 

JL 

4000 

19-3-4- 6- ll-l 

2 

2 

4000 

1-28-22-23-27-1 

3 

4 

4000 

1-10-8-2' '7- 26-1 

4 

6 

4000 

1-21-20 -18-17 -1 

5 

7 

4000 

1-25-5-19-16-1 

6 

r>, 

\ ) 

4000 

1-24-15-14' 13-12* 1 


Total, cost of distribution = Rs. 9494 » 2 per month 


Percentage Variation in Winter - morning Demand Set 

Effect of Variation in Winter-morning Derma nd on Optimal Cost of 
Transportation and Optimal Distribution Cost. 


TJ 

ih‘ 

t 


Optimal Total Distribution Cost (Rs ./ month) 


•xj 

03 

iD 

o 

o 

o 

O 

o 

o 

o 

O 


o 

o 

o 

o 
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EFFECT OF VARIATION IN .SUMMER, MORNING DEMAND SET ON 
OPTIMAL COST OF TRANSPORTATION AND TOTAL COST OF 
DISTRIBUTION 


Percentage 

Variation 

Cost of transportation 
(Rs. /month) 

Total cost of 
distribution (Rs. /month) 

- - 2 C % 

1122 

9262 

-1558 

uni 

9321 

-1055 

1172 

9506 

” 5 % 

1182 

9516 

0 % 

1239 

10768 

5 % 

1294 

11016 

10 % 

1320 

11043 

15 % 

1330 

12248 

20 % 

1360 

12472 
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Table 7.11 

Recommended routes and vehicle-dispatching strategies for the variation 
of demands at each customer upto +10$ of present summer-morning demand set 


Route 

Vehi 

.cle- -allocation 

Route path 


No. 

Capacity 


1 

1 

4000 

1-5-27-25-1 

2 

2 

4000 

1-2-7-11-24-1 

3 

4 

4000 

1-10-17- -18-16-1 

4 

6 

4000 

1-21-20-15- -14-13-1 

5 

7 

4000 

1-9-28-3-26-1 

6 

■ P: 

4000 

1-23-22-19-1 

7 

9 

4000 

1-4..-6-8-12-1 


Total cost of distribution = Rs 


11043.30 per month 




O CD CO r- cn 

(q^uouu/ s^) jsoo uo!)nqij|S|Q \ d\qi (Dujgdo £ 


Percentage Variation in Winter - morning Demand act 

iffect of Variation in Winter- morning Demon.! or Optimal Coot 
Transportation and Optimal Distribution Cost. 
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which are recommended for the variation in demands at each 
customer anywhere between - 6 % and - 20 %. Total cost of dis- 
tribution, for these recommended routes and policies, is 
Rs . 8o48/- per month, whereas the optimal cost of distribu- 
tion, for the demands reduced by 10$, is Rs . 7983/- per month. 
But this difference is net very significant. Therefore, to 
avoid frequent changes in vehicle dispatching strategies with 
change in demand, such vehicle dispatching strategies are 
recommended for the practical purposes. Similarly, the results 
have been obtained for the variation in Summer-Morning demand 
set as indicated by Tables 7.10, Table 7.11 and Pig. 8. 

7. 5 Scope for further work 

The work outlined in the previous pages attempts at 
developing methodology for solving single depot and single 
product delivery vehicle dispatching problems. For more 
real life practical problems, the cases of multi-depot and 
multi-product delivery should be included to develop a general 
model . 

In the problem, as formulated above, it was assumed 
that at every delivery point the demand must be satisfied by 
one delivery. As a further extension of the proposed model, 
this condition, may be relaxed by permitting several routes 
to pass through a delivery point and partial delivery at the 
delivery point by the vehicles. Proposed model does not take 
into consideration the traffic constraint. This constraint 
may be taken into account for developing a generalized model. 
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* # *•#*»■** *#*■*•* ■»#«*# *«*#***« «*•*•»•*■»*•¥****•*-»**•*«#•»*» 

* PROGRAM 1 FOR CONSTRUCTING INITIAL ROUTES * 

** ************ *******•**«•***■*•*****•«***«**** **■#***«■*« 

DI PIN SIGN KCQ(AC,4C>,MlN5gQ(20,4C) ,Z (SO) ,NGC( 20) ,L INSf Q 120,40) 
DIMENSION ID(40,4C),MMTP(40) ,LEG <20 ) * IES (20 } , LIN (50) , LF S( 20 ) 

DI MANSION SPTS(4C) ,SPPS (40) « DC (4 0*40 »,JC81( 40), JCB2( 40) 

DIMENSION KINSEQ <20,40 ) , KKEQ< 20 ) 

DIMENSION KDSD(;2,a?) 

INTEGER Z,ZZ, ZLE ,ZUB, DC, SPTS,SPPS 
COMMGN/SLM/MI NSFG, KK 
COMMCN/RnC/LI nsfq*kod,n 
COMMON/VSS/CD 
COMMON/NFEN/NV 
PFAC 60S'.., A 
6 0S-' FORMAT ( FF. I ) 

490 READ I ,NC 
i FORMAT (13) 

IF (NC. EQ.IOG ) GC TO 9000 

NCC=NC+1 

A=A+5. 0 

PRINT 6053, NC 

605? FORMAT (iHl ,// ICX »*NG. OF CUSTOMERS = *» 14 / ) 

PRINT 6054 

6054 FORMATl //2 OK, »D I STANCE MATRIX*//) 

PRINT 6055, (I I ) , I =1 ,NCC I 

6055 FORMAT <3X, 3214 ) 

DO 101 1=1, NC 
JJ=I+I 

RE AC 102, ! KDSDU ,11) ,II = JJ,NCC) 

101 CONTINUE 

102 FORM AT ( 201 4 ) 

CO 1C3 1=1 , NC C 
00 103 J=1 » NCC 

IF ( J.GT . I ) GO TO 103 
IF(J.EQ»I)GO TO 104 
KDSD (I , JI = KCS D( J, I ) 

GO TO 103 

104 KDSD ( I , J ) = 999 9 

103 CONTINUE 

DO 6020 1=1, NCC 

PRINT 6030, I, { KDSD ( I » J I » J=l, NCC I 

6020 CONTINUE 

6030 FORMAT (13, IX, 3214) 

6CC0 READ 2, NV • 

2 ■ FORMAT ( 13) 

IF (NV.EQ.5GC) GO TO 490 
PRINT 6007, NC ,NV 

6007 FORMATl 1HI, // /9X ,*N0. OF CUST MRS =* , I 3/9X * *N0. CF VEHICLE 3 *, 13/) 
NVV=NV+l 
N=NV+NC 

DO 6070 1 = 1 , NV 
DO 6070 J=1,N 
IF (J.GT.NV)C-O TO 6080 
KDD ( I * J ) =9999 
GO TO 6070 
6080 JJ = J-NV+! 

KDD II , J)=KDSD ( | , JJ ) 

6070 CONTINUE 

PRINT 5050 


PR0C0C3 

PR0CCC4 

PROCOCZj 

PR0GGC6 


PROGOI 1 
PROCOld 
PRO CGI Jj 
PROCC ^ 

pp.occM 

PROCC?/ 

prooot^ 

PROOOl* 
PROOC } 1 
PR 00 0 It * 
PROOOl/ 
PPOOO i 
PROOO ^ 
PROOO ? 
PROOO < 
pRocor 
PROCC 1 1 
PROGOI' 
PROCCI 1 

procov 

PROGOI 1 
PROOO' 1 ' 

prooo r 

PROOO? 

PROGOI: 

PROOO? 

PROOC? 

PROOO? 

PROOO? 

PROOO? 

PROOO-: 

PRO COT 

PROCC' 

PROCC > 

PROOO/ 

PROCC?: 

PROCCI 

PROOO? 

PR0C02 

PRO 002 

PRG002 

PROOO? 

PROCCI 

PROCCI 

PRO 002 

PROOC' 



5050 

FORMAT! //I 0Y* * INITIAL ROUTES BY MULT2PLF TRAVrt-L ING SALESMEN PRO* ) PR0QC?4 


DO t 090 I=NVV,N 

PR0C025 


II=I~NV+X 

PR0CG25 


DO €090 «J=1 * M 

PRO CO'" 6 


IFU.GT.NV1G0 TO 6095 

PROCC ? 7 


KCCU, J)=KDSD( II ,1) 

PPOCC’* 8 


GO TC 6090 

PROCC" 9 

6095 

JJ=J»NV+i 

PR0CC7C 


KDDCI,J)=KCSD(II,JJ) 

PROCC 1 

5090 

CONTINUE 

procc:: 2 


IF {NV.EQ.2IG0 TO 608^ 

PROCCF 7 


GO TO 6088 

PROCO -2 

€081 

PRINT 740, UI ),I=1,N) 

procc"’: 7 

140 - 

FORMAT (//15;<» *AU GMF.NT r D DISTANCE MATRIX*, //4X»*0( 13) ) 

DO 730 1=1., N 

PRINT 720, I* ( KDDCI,J), 

PROOO J A 

730 

continue 


120 

FORMAT (I ,2X,40l3) 

P ROC CO 8 

6088 

00 5 1*1, N 

CO 5 J=1 , N 
•CDU* J}*KDOU t J » 

PROCC’’*: 

5 

CONTINUE 

CALL ASGNMIMlNSEfif 2Z, KN,NEQ#N) 

IF (KN. EG *1000) GO TO 9008 

PRINT 710 


710 

FORMAT ( //15X, ^INITIAL ASSIGNMENT GIVES FOLLOWING SUBTOURS#,/) 



PRINT 6079 

PR 00 OF! 

6079 

FORMAT t/lX**SUBTGUR«, ?)! ,*PATH SECUENCF*,/) 

DO 1 1*1, KN 

KP S = NEQ ! I ) 

PRINT 8»I'» (MXNSEQf I,J),J*1,KPS> 

PR0004! 

8 

FORMAT (IX, I?, 5X ,41(13 ) ) 

PROG05< 

7 

CONTINUE 

PRINT 9,ZZ 


9 

FORMAT !/1.0X,# INITIAL LCfcER 80UN0=*,I5> 



CALL FEASt MlNSEQ ,NFQ, KN ,IES,KK) 

IF (KK.EQ.O)GO TO 500 

F-R0C09I 

C*****FINEING INITIAL UPPER BOUND 



KN Z=0 

PROCOSl 


DO 10 1=1, KN* 

KP5P=NFG!I ) 

DO 11 J=l, KPS P 

L1NSEQ! I, J )*MINSEQ<I, J) 


11 

CONTINUE 

LEQCI )=KPSP 


10 

CONTINUE 

LL =KK 

DO 12 1=1, KK 

LES(I)=IKS(I) 


12 

CONTINUE 

KMN=KN 

CALL UPRBNDIK ZUB »LEQ» LEG, LL, KMN ) 

ZUB=KZUB 

PRINT fclO, ZU8 

1 

810 

FORMAT! //I CX, *B€LOW TOURS BY INITIAL HE UR I STI CS* , 1 OX , «U * B . =* 1 5 / ) 

DO 811 1=1, KMN 

joii=Lea{i) 

PRINT 812, IL1NSFQU »J> ,0=1, JOIT) 

PPOCC’I 

1 

812 

FORMAT! 5X, 41! 13 ) ) ' .. 

PROCC ’-j 



£11 CONTINUE 

IF (ZUB» ?Q. ZZ) GO TO 600 

GO TO 650 

TOC PRIM i3 

13 FORMAT ( / X»*ABCV£ TOURS ARE FINAL AND MIN CIS IS INITIAL L.B.*> 

GO TO 6000 P1QCG-? 

6 GO PRINT IS 

16 FORMAT ( / /I0^, *F I NAL TOURS ARE AS OBTAINED BY INITIAL HEURISTICS'*} PROCCV- 
GO TO 6000 PR0CC6S 

£5C DO £63 I =X , N 
DO 663 J=I » N 
IC(I, J)=KDD(I ,J» 

€63 CONTI NL'- 

7CC CALL SHCRT ( IE5*NEG »HMTP , LMS, LIN! 

LNS=LMS»X 

MM=C 

DO I?0 I =1 1 LNS 
1 1 = I + ? 

SPTS(MM)=MMTPtI) 

SPPSl MM) = MMTPtII } 

120 CONTINUE I 

KONT=G PROOlCjj 

DO 130 JHUM=1 t MM 
DO 131. 1=1, N 
DO 131 J=1,N 
DD ( I , J )=ID ( I, J) 

131 CONTINUE 

KTEI,=SPTS( JHUM) 

KTE?=SPPS ( JHU M) 

DO ( KTEi»KTE2) =9999 

CALL ASCNMt MI NSE Q,ZZ* KN * NEQ, N! 

IF(«N.FQ.1CC0)GC TO 0003 

IF IZZ.GE.ZU8I GO TO 135 

KONT=KONT+1 

Z C KONT ) =ZZ 

JO El (KONT) =KT El 

J0B2 (KCNT)=KT E2 

CALL FEASt MlNSEQ*NEQ,KN,IES t KK» 

IF (KK.EQ.O)GO TO 140 
GO TO 130 

140 ZU8=ZZ 
KNZ=KN 

DO 7054 1=1, KNZ 
7054 KKEQU)=NEQU ) 

CO MJ 1=1, KN Z 
KPSE=KKEQ( I ) 

DO 746 J=l, KPSE 
KINSEQt I, J )=MINSEQ II , J) 

746 CONTINUE 

1.42 FORMAT ( IX, 13# 5X ,41 13 ) PROCX*- 

141 CONTINUE . I 

GO TO :,30 

135 CONTINUE PROO'iA- 

130 CONTINUE 

IFlKONT.gQ.OIGC TO 2000 
CALL MIN(Z,KONT ,KENT, MKONT) 

ZL£=MKONT 

IF IZLB.EQ.ZUB )G0 TO 2000 



KTf -* jcb;<ksnt» 

KTF:6 = JCB2 ( KENT ) 

DO 1 95 1=1 » N 
DO 695 J=1 ,h 
CO ( I , J ) = 7 D { J, J) 

695 CONTINUE 

CC1KTE3, KTE4 >=9999 

CALL A5GNM{MlNSEQ,ZI» KN»NEQ, N) 

IFlKN.SQ.iOOO JGO TO 90Cg 
CO (9> 1=1 , N 
DO 691 J = 1,N 
ID (I f J ) = CD ( I » J ) 

69 2 CONTINUE 

C/'LL FFAS< MIN5EQ,NFQ,KN,IES,KK> 

IF i KK.FQ.O )€0 TO 3C00 
CO TC TOO 

3000 PRINT 3003 * ZZ 

3003 FORMAT l / OX,*FINAL TPURS AR C AS FOLLCW5*» r JX,*ZZ=*I5 > 

PRIM 60' 9 
CO 2C01 1=1, KN 
KPZS=NEC( I ) 

PRINT 3002,1, IMINSEQCI, J), J=1,KPZ5) 

300? FORMAT ( IX. IB, 5X,M<I’> > 

30CI CONTINUE 

GO TO 6000 

2CC0 IFIKNZ.EQ.QJGO TO 7808 
PRIM 2001, ZUB 

2 0C1 FORMAT! //20X, * FOLLOWING TCURS ARP OPTIMAL*. 10X,* KIN DIS =* , I 5 ) 
PRINT 6079 
DO 2CS0 1=1, KNZ 
KPS£=KKEQ( I J 

PRINT 2060,1, (KlNSFQd, v}J,J=l»KPSE> 

2 06 0 FORMAT ( U, 13, IX ,41 <13 ) ) 

2050 CONTINUE 

GO TO 6000 
1808 PRINT 16 

GO TC 6000 
9008 PRINT 9007 

9007 FORMAT {//IX,* MISTAKE IN FINDING MIN. LIN. REQC. CHECK SUBR .«» 

GO TO 6000 
9000 STOP 

END 


PROG! 


PR 001 

PR001 
PROC'I 
PR 00] 

PROG;: 


PROG] 

PR001 

prog: 

prog: 

prog: 


PROCl 



G250 



FORTRAN SOURCE LIST 


ISN 


SOURCE STATS ME' NT 


| 

0 

$IBFTC 

F EAS 


PRC 

1 


SJBROJTINM F ASININSL 

Q,K£Qt KLN, KES » M K ) 

Pkf 

2 


C3MM0N/NE? M/hV 


PK) 

3 


01 MENS I ON fsll MSE y { 2C ,4 

) , K£Q( 20) , Kc 51 20 ) 

P r ^ 

& 


M< =0 


P R ( 

5 


DO 5 1 = 1, KLN 


P R f 

6 


KPS=KSQ{ I ) 


Pk! 

7 


DO 6 J = 1 , KPS 


P P ! 

10 


KF T =NI NSEQ ( I , J) 


PR? 

11 


DO 7 J J=1 » NV 


P P 1 

• 12 


IF (KTT.EQ. JJ )G3 TO 5 


PRi 

15 

7 

CONTINUE 


' PRf 

17 

6 

CONTINUE 


P R( 

21 


M<=MK+1 


PR! 

22 


KES (MK J=I 


P kl 

23 

5 

CONTINUE 


p rj 

25 


RETURN 


P R! 

' 26 


tN D 


PPJ 

S250 



I 8 MAP ASSEMBLY FLAS 


NO MESSAGES 

FOR ABOVE ASSEMBLY 





-G250 

ISN SDURCt STaTlK Mf 


FORTRAN SOURCE LIST 


0 SIBFTC SHORT 

1 SJ 8R0UTI NE SHORT ( IE Si EQ t MTP » I MS* KIM ) 

2 DIMENSION IES(20), NEG(20) »MTP(*-0) , KIM (50) »MINSE : Q(2D, 

3 CD MMDN/SLM/MINSEQj KK 

4 KT =0 

5 DD II 1=1, KK 

6 m = lE5( I ) 

7 KI N( I ) =NLG( INS) 

10 II CD NT I M U~ 

12 CALL NIN(KIN|KK,KTi MF J 

15 JTP=ILS{ KT ! 

14 IM S=iNc L( JTP ) 

15 DD 2 J=1 , IMS 

16 MTP( J) =MIN5EQ( JTP, J ) 

17 2 CONTINUE 

21 RETURN 


22 END 

1LG250 


IB HAP ASSEMBLY SHORT 


) 


PR 
P R 
P R 
P R 
P R 
P R 
PR 
P ft 
? R 
PR 
PR. 
PR. 
PR 
P R 
PR 
PR 
? k 


MO MESSAGES FDR ABOVE ASSEMBLY 



G2 50 FORTRAN SOURCE LIST 


ISN SOURCE ST AT..- ML NT 

0 $1 BFTC J PRBND PHi 

1 SUBROUTINE. UPRB ND { K ZU t:> , L EQ , LES f LL, KMN 3 PP> 

2 C3MM3N/R.UM/IP HR' 

3 DIMENSION Li ;:S r w(2 r ,4 ) » KDD( 4C » *V0 ) ,Lf Q J 20) , LrrS 1 2D » 1 P < 40* • •<' ) PR- 

4 DIMENSION I Ft AS ( 20 > ,1 NFEAS ( 20 ) PR S 

5 COMMON/RNC/LlKSGQ, KCD, N PE 1 

6 D3 1 I = 1 , N PM 

7 D3 I J=1,N PR 

10 1 3 ( I , J ) =0 PR 

11 1 CONTINUE P El 

14 CALL PM AIR { KMM , L •* Q ) P R 

15 103 1=0 PR : 

16 K= 0 PR 1 

17 03 2 1=1, KMM . PR: 

20 DO 3 J = 1 , LL PR 

21 IF { I.EQ®LLS( J ) ) GO TO PR- 

24 3 CONTINUE PR: 

26 K= K+ 1 PR 

27 1~ EAS( K) =1 PR 

30 G3 TO 2 PR 

31 4 L= L + l PR 

32 IN FEAS ( L 5 = I PH 

33 2 CONTINUE PR 

35 CALL S EARCH( L , < , I NF EAS, I FEAS , LEQ, N 1,N2,N3 , N4, 1 FFST , I EERT 5 PR! 

36 I? ( N 1» N3 ) = 1 PH 

37 IP ( N 1, N 2 ) = 2 • PR 

40 IPIN4, N2) = l PH 

41 I 3 (N4, -M3) =2 PH 

42 CALL MNSQ( IP,KDD,N, LINSeQ, LEQ ,KMN» 1 2UB ) PR 

43 CALL FEAS(LINS=Q» LErCt KMN, LES, LL) PR 1 

44 IF { LL* LQ*0 )G0 TO 1005 PR 

47 G3 T3 100 PR 

50 1035 K.Z UB=I ZUB PR 

51 RETURN PR 

52 END 3 R 

MEG25Q 1 8 MAP ASSEMBLY UPRBND 


NO MESSAGES F3R ABOVE ASSEMBLY 



FORTRAN SOURCE LIST 


:.32 50 


ISN 


S3 URGE STATE Ml ..NT 


0 

$IBFTC 

PMATRX 

P R f 

I 


SUBROUTINE P?’,ftTRX{ KMN » LEQ) 

p R f 

2 


01 MENS I ON L I NSE Q ( 2 0 ,40 ) , I P ( 40 , 40 ) , L&'Q < 20 ) » KDD ( 40 , 40 ) 

P K i 

3 


COMMON/RNC/L > NS ~Q, K CD , N 

n 0 
r t\ 

4 


COMMON/RUH/I P 

Pftf 

5 


03 1 1=1, KHM 

PRt 

6 


KP S=LE Q( I ) 

PKS 

7 


K? N=KPS“ 1 

jl^ 1"? f 

10 


03 2 J = 1 , K p N 

PR! 

11 


KT T=L1 NSEQ ( I , J) 

P ft! 

12 


JJ =J + 1 

P Rt 

13 


KT M=LI NSfr Q ( I s JU } 

P K 

14 


1? { KTT , KTM ) = 1 

PR] 

15 


03 3 K = 1 , N 

PRi 

16 


IF ( K.EQ.KTM) GO TO 3 

PR 

21 


I 3 { KTT j K } =2 

PR' 

22 

3 

CONTINUE 

PR 

24 


DO 4 K. K = 1 , N 

PR 

25 


IF (KK. EQ. KTT )G0 TCi 4 

PR 

30 


IP ( KK, KTM) =2 

PR 

31 

4 

CONTINUE 

PR 

33 

2 

CONTINUE 

PR 

35 

1 

CONTINUE 

PR 

37 


RETURN 

PR 

40 


END 

PR 


'1LG250 1 8 HAP ASSEMBLY PMATRX 

NO MESSAGE S FOR ABOVE ASSEMBLY 



G250 ' FORTRAN SOURCE LIST 

ISN S3 URGE ST«TT-k:NT 


0 $1 BFTC A NSQ 


1 


SJSR3U7IN, M:;SQ<PiMD» ,MANSLQ, N£Q, KN, IZ 3 

P R ! 

2 


DIMENSION P(-C, T),FD(4C, 40 ) , MANSE Q( 20, 40) ,M“Q(20J 

P K 1 

3 


INTEGER P 

p r : 



DI MENS ION JI MS” 9( 100) 

P K ! 

5 


N1 = 2*N 

P R 

6 


KN=0 

P R 1 

? 


DO 401 K=1 ,N1 

PR! 

10 


JI NSE Q ( K ) =0 

: PR! 

11 

401 

CONTINUE 

P R( 

13 


KJ M=0 

? R { 

1 4 


K< =0 

PR! 

15 


!Z = C 

P Ri 

16 


DO 402 1 = 1, M 

P R 

17 


DO 403 K = 1 , ii 1 

P 

20 


I F { I . E 3 • J I NS Q { K ) 5 5 C TO 402 


23 

403 

CONTI NUE: 

P R 

25 


KK=KK+1 

P R 

26 


JINSfcQ(KK) =1 

P-. 

27 


DO 404 J = 1 1 N 

PR 

30 


J1 = J 

PR 

31 


IF { PCI * J } ® : Q « 1 ) GO TC >05 

PR 

34 

404 

CONTINUE 

PR 

36 

405 

K< = KK+- 1 

PR 

37 


JINSEQ(KK) =J1 

PR 

40 


I Z = IZ+-MD ( I , J 1 ) ' 

PR 

41 

406 

J2=J1NSEQ( KK ) 

PR 

42 


DO 407 J = 1 , N 

PR 

43 


I F ( P { J 2 , J ) » t Q . 1 5 GO TO 410 

PR 

46 

407 

CONTINUE 

PR 

50 

410 

K< = KK+ 1 

PR 

51 


JINSEGI KK)=J 

PR 

52 


IZ = IZ+-MD(J2, J) 

PR 

53 


IF ( J.EQ.I >G0 TO 431 

PR 

56 


GO TO 406 

PR 

57 

431 

KN = KN+ 1 

PR 

60 


Kl M=KK~KJM 

PR 

61 


KJM=KK 

P F 

62 


NEQ ( KN ) =KI M 

PR 

63 


KZ=KK-KN 

PP 

64 


IF (KZ.GE.N)GO TO 590 

PE 

67 

402 

CONTINUE 

PP 

71 

5 93 

JSS=0 

PE 

72 


DO 1001 1=1, KN 

pf 

73 


NERC=MEQ(I ) 

PI- 

74 


DO 1003 J = l» NERO 

PE 

75 


K= J+JSS 

PI 

76 


MANSEQ ( I , J ) = J I NScQ ( K) 

■ P 

77 

1003 

CONTINUE 

Pt 

101 


JSS=K 

Pi 

102 

1001 

CONTINUE 

P 

104 


RETURN 

PI 

105 

HEG250 


END 

IBMAP ASSEMBLY MNSU 

Pi 


MO MESSAGES FDR ABOVE ASSEMBLY 



FORTRAN SOURCE LIST 


G250 

ISM S3 URGE STATE Mr. M 


0 

$1 8FTC 

S EARGH 


r’ K 1 

I 


SJBR3JTXN,: S ARCH( L »Ki INFcAS » 

I Ft AS »NEa* Nl, N2» N ,M- f IFFST .13) 

PrJ 

2 


C3MMOM/RNC/L INS EQ» KCD, N 


p C ( 

3 


DIMENSION LI MS r QC 20 >4- 3 , KDDC 40,40) , I NFL A S ( 20 ) , I F.„ AS { 2C ) , 0(2'.,) 

p f j 

4 


01 MENS ION IX PC 40 5 , K IXPC40 ) » IPXC40 ! »KI PXC 40 ) s IF ST ( 20 ) , I . RT ( 20 ) 

P R i 

5 


DI MENS ION 1 1 PXC 40 5 , KKI PXC 40) , 

11 C 20 ), 12( 20) »r. (20) , I-vC20) 

PK! 

6 


01 MENS I ON ND I SC 50 5, MI , Dl S ( 50 ) 

, M I NM UM { 5 0 ) , I F R T ( 2 0 3 , HI N I NC ( 2 0 ) 

; P K i 

7 


DI MENS I ON J 1 ( 20 5 , J 2 C2<. ) , J 0 C 20 ) , J4 ( 20) 

i*? \ 

10 


KZ =0 


p K 

11 


03 3 I =1, L 


P R- 

12 


IM = 1 NF cAS ( 1 3 


p K 

13 


K?S=N2UC I M } 


P K 

1 4 


Kd=0 


P R 

15 


03 4 J = 1 , K 


PR 

16 


IM = IFE ASC J 5 


9 K 

17 


KPN=N£Q( I N) 


P R 

20 


KRN=KPM-1 


p b 

21 


Kv/ = 0 


PK 

22 


03 5 J J = 1 , KRN 


PR 

23 


Ktf =KV+ 1 


PK 

24 


KTT=LI rCSEQ C I 1 1 J J ) 


PR 

25 


JOJJ+1 


PR 

26 


KFM=LINSEQ( IN, JK) 


PR 

27 


IXP (K.V ) =KTT 


PR 

30 


KI XP ( < V ) = KTM 


PK 

31 


K3=0 


PR 

32 


D3 6 11=2, KPS 


PR 

33 


KT S=LI MSEQ C I M* 1 1 ) 


PR 

34 


I < = 1 1 - 1 


PR 

35 


KFN=LI NS HU C IM, IK) 


PR 

36 


I3IS=KDDCKTT,KTS) 


PR 

37 


J3I S=KD0 ( KTN, KTM ) 



40 


K3 I S = KOO ( KTN i KT S ) 



41 


K3=KQ+- 1 

"■ ■' ' ; . ' ' , ■ ; 


42 


N3IS(KQ)=IDIS+J0IS-KDIS 


m 

43 


IPXCKQ )=KTS 


PR 



KIPX(KQ)=KTN 


Pf 

45 

5 

CONTINUE 


. PI 

47 


CALL MIN(NDIS»KQ»KCTC t MET) 


P i' 1 

50 


II PXCK v)=IPX(K3Q) 


Pf 

51 


K<IPX(KV) =KI PXC KQQ) 


PF 

52 


M3 I S = KDD ( KIT, KT M 3 


Pf 

53 


mindisckv)=m: t-mdi s 


Pr 

54 

5 

CONTINUE 


Pf 

56 


CALL MIN(M1NDIS,KV, KVV , HER ) 


P! 

57 


K/^ =KW+ 1 


P? 

60 


MI NMUM ( KW ) =MHR 


Pf 

61 


11 ( KW ) = I XP ( K. VV ) 


Pf 

62 


12 CKW) =KIXP( KW ) 


Pf 

63 


13 C KW) =1 IPXC KW) 


Pf 

64 


14 CKW) =KKI PXC KVV ) 


Pi 

65 


I=RT(KW) =IH 


Pf 

66 

4 

CONTINUE 


Pi 

70 


CALL MIN ( MINIMUM , KW , KWW.MEP ) 


P i 

71 


KZ = KZ+ 1 


P 


G25C 


ISN 

S3 URGE STAi: 

ML i T 

12 

Ml IN I NO { K.Z ) = M 

IP 

7 3 

JI ( KZ ) =1 1 { KW W ) 

7 4 

J2 (KZ) =12 (KW 

w ) 

75 

J3I KZ ) =1 3 ( KWW 5 

76 

J4 ( KZ ) = 1 4 { KW W J 

77 

IFSTt <Z) =1 FR 

7 I K W W 5 

IOO 

IERT(KZ)=IM 


101 3 

C3 MTI NU'v 


103 

CALL MINIMI N 

I .INC i KZ } KZZ, MEM) 

10 4 

MI ra = M,;f-1 


105 

Ml = Jl ( KZZ ) 


106 

N2 = J2( KZZ ) 


107 

N3=J3( KZZ ) 


110 

N4=J4{ KZZ) 


111 

IF F ST= I FST I K 

ID 

112 

IE £RT= IERT ( KZZ } 

113 

RETURN 


114 

r. N D 



FORTRAN SOURCE LIST 


IB MAP ASSEMBLY SEARCH 


NO MESSAGES FOR ABOVE ASSEMBLY 



S3 URGE ST A 


G2 50 
ISN 



FORTRAN SDURCS LIST 


0 

SIBFTC 

A 5GMM 

P P. 

1 


SUBROUTINE 4SGNMCMANS Q, IZ, KN*Nf?Q, N> 

P g 

2 


IMTCGi- R DD*P, UNMRKR ,UC,UR,UZ,D 

P R 

T 


CD MM3N/VSS/DD 

PP 

4 


C3MM3M/JLM/P 

P K 

5 


C3MMQM/ZLH/D 

p K 

6 


DI HENS ION DD(40»40> , PC 40,40) ,UNMRKR<40 ) , MRRKC 150) 

? R 

7 


DIMENSION DC < 0,40), IP0KC40) , I SOKC40), N DI S U 00 ) , I X P ( 1 0 01.1PM1 0 2 ) 

PR 

10 


DI MENS I ON MACS? QC 20 ) , NE Q( 20) 

P R 

11 


D3 713 1 = 1, N 

PR 

12 


03 710 J = 1 , N 

PR 

13 


D( I , J) =ODU, J) 

PR 

14 

710 

CD NT IN Ut 

P R 

17 

163 

CALL Z i£ RO(L) 

PR 

20 


DD 107 1=1. N 

P.R 

21 


DD 107 J = 1 , N 

P R 

22 


PC I , J) =0 

PR 

23 

107 

CONTINUE 

PR 

26 


CALL ASSIGN!:;, IPOK, ISCiK,KOK) 

P R 

27 

995 

I N DEX=0 

P R 

30 


DD 1 23 1=1, N 

PR 

31 


D3 123 J = 1 , N 

P r 

32 


IF ( P( I , J ) . ,.Q. 1 ) I MDC X=IND£X+1 

PR 

35 

123 

CONTINUE 

PP 

40 


IF (INDEX. F«*"J)GO TO 110 

PH 

43 


CALL M I NL I M C UNMRKR , LR* MRRKC,UC,N ) 

PR 

44 


UZ=UR+UC 

PR 

45 


I F ( I NDLX • E: Q. UZ ) GO TC 150 

PF 

50 


IF (INDEX. GT.UZ)GD TC 118 

PR 

53 


G3 TO 141 

PF 

54 

113 

PRINT 119 

PH 

55 

119 

FORMAT! //10X, -SUM TFI:;G WRONG IN FINDING MIN® LINGS Rf'QD.*) 

PH 

56 


GO TO 5000 

PH 

57 

141 

IF (UZ. LT . N ) GU TO 15C 

PH 

62 


GO TO 166 

PH 

63 

153 

CALL DRZERCUM, UNMRKR, MRRKC,UR,UC) 

PI 

64 


GO TO 160 

P! 

65 

110 

CA LL MMSQ( P, DD, : 1 ,MANSLQ, NEQ, KN, IZ ) 

Pf 

65 


GO TO 500 

PI 

67 

1 66 

DO 923 1=1, N 

Pf 

70 


KA K 1 = 0 

PF 

71 


DO 921 J = 1 , N 

PI 

72 


1= ( P( I , J).L ; 0.1)GD TC 920 

Pf 

75 

921 

CONTINUE 

P! 

77 


DO 922 K=1 » N 

Pi 

100 


13 I S = 0 D ( I , K ) 

Pi 

101 


DO 923 KK= 1, N 

P 

102 


1= { f(KK,K).EQ.l )G0 TO 924 

P 

105 

923 

CONTINUE 

Pi 

107 


GO TD 925 

P 

110 

924 

DO 950 J = 1 , N 

P 

111 


IF ( FlKK, J ) .5U«1 )G0 TO 950 

P 

114 


DO 990 JP S=1 , N 

P 

115 


I = { PC JPS, J ) . I Q. 1 )G0 TO 950 

P 

120 

993 

CONTINUE 

P 



i..G2 30 FORTRAN SOURCE LIST A5GNM 


ISN 


SOURCE ST A" NT 


122 


JO I S=DD ( KK j J ) 

r • ; 

123 


GO TO 960 


124 

9 50 

CONTINUE 

PH - 

126 

960 

MO I S=DD ( KK } K ) 

■ 

127 


KAKI=KAKI+1 

• 

133 


I< P (KAKi } = KK 

pf 

IS 1 


PXCKAKI )=J 

- • , 

132 


NDISCKAKI )=IDIS+JDI S-MDIS 

: , ! 

133 


GO TO 922 

H j 

134 

92 5 

KAKI = <aKI + 1 

P 

13 5 


NOlStKAKI ) =00 ( I i K ) 

PR 

136 


IX P CKAKI ) = 100 

PR 

137 


IP X C K4KI ) =0 

P / ; 

140 

922 

CONTI MUE 

P K 

142 


CALL M I N U'i D 1 S * K A K 1 * KIXl,MOA) 

v f ; 

143 


P( IjKIKI ) = 1 

PH 

144 


DO 930 I K= 1 * .0 

P R 

145 


IF UK. EQ.IJGG TO 900 

PR 

150 


PC I K j < I KI i =2 

p-i 

151 

930 

CONTI MUE 

P:-| 

153 


DO 940 J = 1 j N 

PR 

154 


IF{ J.SQ.KIKI )G0 TO $40 

PR 

157 


PC I, JS =2 

Pi- 

160 

940 

CONTI MUE 

P R 

162 


LSP=IPXCKIK1 ) 

?.< 

163 


Lrp=lXPCKIKI ) 

PR 

164 


I=(LTP.EQ. 100)30 TO 920 

■ PR 

167 


PC LTP» LSP ) = 1 

Pi- 

170 


DO 970 II = l,,j 

p r< 

171 


IFCLTP.EQ.il ) GO ID 970 

■ PR 

174 


PC I ] j LSP ) = 2 

P F 

175 

970 

CONTINUE 

PR 

177 


DO 933 J J = 1 * M 

PR 

200 


I FCLSP.6Q. JJ)GO TO 9,.0 

pJ 

203 


PC LIP, J J ) = 2 

PH 

2D4 

98 0 

CONTINUE 

PH 

206 

920 

CONTINUE 

PP 

210 


GO TO 995 

PH 

211 

5000 

KM =1000 

PH 

212 

500 

RETURN 

PH 

213 


END 

PH 


MLG250 


IBMAP ASSEMBLY ASGNM 


MO MESSAGES FOR ABOVE ASSEMBLY 



G2 5u 

ISM SOURCE STATEMENT 


FORTRAN SOURCft LIST 


0 

1 

2 

3 

4 

5 

6 
7 

11 

12 

13 

14 

15 
20 
21 
23 

25 

26 
27 

31 

32 

33 

34 

35 

40 

41 
43 

45 

46 

MeG250 


SI8FTC l PRO 

5JBR0JTI Ml Z RDU'H 

IN! TEGE R D,BD, R3 W , CO L,CA» RMIN, CM] N 

DIMENSION D ( 4 C ? s-C I * BO t 50 )» ROW C 40.) , COL 1 40 ) 

C3MM3N/ZLM/U 

03 1 I=1,N 

03 2 J = 1 » M 

2 ED { J ) = D U f J ) 

CALL M1NCB0,N*<T,RMXN) 

ROW (I ) =RMIN 
03 3 J = 1 j H 
CA=D< I j J ) 

IF(CA.EQ. 9999)30 TO 3 

D( I , J) =D( I , Jl-RMIM 

3 CONTINUE 
1 CONTINUE 

D3 4 1=1, N 
03 5 J = 1 , N 

5 83 { v ) = D { J , I ) 

CALL MINI BD, N , K T, CM IN ) 

C3 L { I ) =CMI N 
DO 6 J = 1 , N 
CA = D ( J , I ) 

1= (CA.EQ. 9.999) GU TO 6 
D( J , I ) =D ( J 1 1 ) -CM I N 

6 CONTINUE 

4 CONTINUE 
RETURN 
END 

IBMAP ASSEMBLY ZtRU 


P R ! 

P K 
p 

P R 
PR 

p | 

P R 

P ft 

PR 

PR 

PR 

PR 

P R' 

PR 

P P 

? P 

P F 

Pf 

PS' 

Pf 

PI- 

PR 

PR 

Pf 

Pf 

Pf 

Pf 

Pf 

Pf 


NO MESSAGES FOR ABOVE ASSEMBLY 


G2 5P 
ISN 


0 SIBFIC 

1 
2 


4 

5 

6 
7 

10 

11 

12 

13 

1* 

15 

16 
20 
21 

24 

25 
30 
32 

35 

36 

37 
4 2 

43 

44 
47 

50 

51 

52 
55 
60 
61 
63 

65 

66 

71 

72 

73 
7 % 

75 

76 
101 
102 

104 

105 
110 
111 
113 

115 

116 
117 
120 
122 


99D 


201 


6 

11 


501 


102 

2 


8 

1 


302 


SOURCE STATEMENT 
ASSIGN 

SUBROUTINE ASSI GN { N ,1 POK, I SDK ,KOK) 

COMMON/ JLM/P 
COMMON/ ZLM/D 
INTEGER L)»P 

DI MENS 1 ON 0( '--Of -0 ) , P { 0, 40 ) f L COUNT ( 50 ) , 1 PUK (40 ) , 1 SUK ( 

KO K = 0 

L< =0 

I N DEX= 1 

NJ M=0 

MJM=0 

DO 1 I = 1 , N 

DO 201 J = 1 , N 

LCOUNT ( J ) = 0 

CONTINUE 

DO 11 K= 1 » N 

IF ( 0 II » K ) • u Q» 0 ) GO TO 

GO TU 11 

IF (P(I jKJ.tQ.O) NUM= NUM+ 1 
CONTINUE 

IF (NUN. LE. INDEX) GO TO 3 
GO TO 1 
10 2 J = 1 » N 

IF (D(I » J).EW.O)GO TO 
LCOCNT ( J ) = 10*0 
GO TO 2 

IF ( P( I , J ) . EQ.O ) GO TO 501 
LCOUNT ( J)=10*N 
GO TO 2 

DO 102 1 MC = 1 »N 
IF(D(IMC*J)«Nt • 0 ) GO TU 102 
IF(P(IMC,J).NL.O)GO TO 102 
LCOUNT ( J)=LCUUNT( J ) +1 
CONTINUE 
CONTINUE 

CALL MIN (LCOUNT , IN, J 1, MINX) 

IF (MINX. EQ. 10*N ) GO TU 1 

P( I,J1)=1 

KO K=KOK+ 1 

IPOK ( KOK ) = I 

I S OK { KQK ) = J 1 

DO 5 1 1 = 1, IN 

IF (II.EQ.I )G0 T 0 5 

P( II, Jl)=2 

CONTINUE 

DO 8 JJ=1,N 

IF(JJ.EQ.J1)G0 TO 8 

P( I,JJ ) = 2 

CONTINUE 

CONTINUE 

DO 301 J=1»M 

DO 302 I = 1 » N 

LCOUNT ( I )=0 

CONTINUE 

DO 303 K=1,N 


FORTRAN SOURCE LIST t- 


\ . G2 50 


ISN 


SDURC;; S T A T , rl. N T 

123 


IFCDCK, J). Q . 0 ) GO TC .04 

126 


GO T3 303 

127 

304 

1= (P(K,J). LQ.U) MUM= FUM+ 1 

132 

3 03 

CD NT I NIUE 

i:h 


IF ( MUM. LE. INDEX )G0 TO 305 

137 


GD TO 301 

140 

305 

DD 305 1=1* N 

141 


I F { D ( I , J ) . Q . o ) GO TC 07 

144 


LCOUNK I ) = 10#N 

145 


GD Tu 306 

145 

? 07 

I F ( P ( I * J ) » t Q . 0 ) GO TC 08 

151 


LCOUNTC I )=10*N 

152 


GD TO 306 

153 

303 

OD 309 JMC=1 t IS 

154 


IF (DU, JMO.NE.OIGO TO 309 

157 


IF ( PC I * JMC ) .Nt . 0 ) GO TO 309 

162 


LCGUNTC I ) = LCD UN T { 1 ) +1 

163 

3 09 

CD NT I NIUE 

165 

2 05 

CONTINUE 

167 


CALL MI N ( LCDUNT * N * I I, UI MX) 

170 


IF (NIMX.EU* 10*N ) GD TO 301 

173 


PC I1,J) = 1 

174 


KD K=KOK+ 1 

175 


I P OK ( KQK ) = 1 1 

176 


1 50KCKDK ) =J 

177 


DD 310 I K= 1 * N 

200 


IF C IK.EQ.I 1)G0 TC 3 10 

203 


PC IK, J)=2 

204 

310 

CD NT I MUE 

206 


DD 311 J K= 1 , N 

207 


I=CJK.t;Q.J)GD TO 311 

212 


PC 11, JK)=2 

213 

311 

C3NTI SUE 

215 

301 

CONTINUE 

217 


LP=0 

220 


DD 900 1=1, N 

221 


DD 903 J = 1 , N 

222 


IFCPCI, J).LQ-1)LP=LP+1 

225 

900 

CD NT! SUE 

230 


IF (LP.EQ.LK)GO TO 9G1 

233 


L< =LP 

234 


GD TO 980 

235 

901 

ISDEX = lNDf:X+l 

236 

9 80 

J3J=0 

237 


DD 960 1 = 1, N 

240 


DD 950 J=1 »N 

241 


IFCD(I,J).rO.O)GO TC 970 

244 


GD TO 960 

245 

970 

IFCPCI, JI.EQ.O) JOJ= JOJ+1 

250 

963 

CD NT I SUE 

253 


IF C JOJ.eQ.O )G0 TO 999 

256 


GD TO 990 

257 

999 

RETURN 

260 


LSD 


MfcG250 


IB MAP ASSEMBLY ASSIGN 


MO; MESSAGES FDR ABOVE ASSEMBLY 


■1..G2 5D FORTRAN SOURCE LIST 


SN 


SOURCE STATE M r N T 


0 

tl BFTC 

MINLIN 

P V 

1 


SUBROUTINE MI ML I N ( UNMRKR, UR, MRRKC, UC, N ) 

PF' 

2 


C3MM0N/ZLM/D 

P K 

3 


COMMON/ JLM/P 

PR 

4 


INTEGER D, P, UNMRKR, LC, UR 

PR 
0 ) P Pf 

5 


DI MENS ION 0(40, 40 ) , P ( 0,40 ), UNMRKR (40) , MRKC( 100 ) *MRKR ( lUi ), MR ( li. 

6 


01 MENS I ON MC ( 100) , JP(IOO) , NINA ( 100 ) , MRRKC ( 50 ) 

Pff 

7 


KM = 0 

P Rl 

10 


KL =0 

P Rf 

11 


IR =0 

P R 

12 


i:=o 

PR: 


c*#* * 

ID MARK ROWS HAVING No ASSIGNED ZERO 

PR 

13 


DD 1 1=1, N 

PF 

14 


DD 2 J = 1 , N 

PR 

15 


IF(P(I,J)*FQ«1)GQ TC 1 

PF 

20 

2 

CDNTINUE 

PF 

22 


KL =KL+ 1 

PF 

23 


MR (KL) =1 

PF 

24 

1 

CONTINUE 

PF 

26 


DD 101 I = 1 j K L 

Pf 

27 


IR = IR+ 1 

Pf 

30 


MR KR ( I R ) =MR l I ) 

Pf 

31 

101 

CDNTINUE 

Pf 


C* »**#TD MARK COLUMNS HAVING ZEROS IN ABOVE MARKE D ROWS Pif 

33 

34 

35 

40 4 

42 

43 5 
44 ' 

47 

50 7 

51 

52 

55 

56 

57 

62 21 

64 

65 22 

66 
67 

70 6 

72 3 

74 

75 

76 

77 102 

101 U«U-U ilcruinnujMir jruM w : * 

102 1= ( KTT.EQ.OJGO TO 11 ? 

105 DO 103 1=1, KTT : 

106 n=IR+l > 

107 MRKR(IR)=JP( I) > 


DD 3 1=1, N 

DD 4 J=1,KL 

IF ( I.EQ.MRC J ) >G0 TO 5 

CONTINUE 

G3 TO 3 

DD 6 J = 1 , N 

IF ( D ( I , J ) • FQ • 0 ) GO TO 7 
GD TO 6 
KM=KM+1 
MS (KM) =J 

1= (KM. (IQ. 15 GO TO 6 

KN=KM- 1 

DD 21 11=1, KN 

IF (MC( KM) •ilQ.MC ( II ) )GO TO 22 

CDNTINUE 

GD TO 6 

KM S=KM 

MS ( KMS )=0 

KM = KM- 1 

CONTINUE' 

CDNTINUE 
DD 102 1=1, KM 

i:=ic+i 

MRKCUG )=MC( ■ ) 

CDNTINUE 

TD MARK ROWS HAVING ASSIGNED ZEROS CORRESP ID ABOVE MARKED COL. 

^ a l 1 rfk'.Ol t M Air 1/ U IO 1/ T T 1. 


Pf 

Pi 

PI 

P 

P 

P 

pf 

p 

p 

p| 

p- 

p 

R: 

p 

p 

p 

p 

p 

Pi 

P 

P 

P 



S3 URGE STATEMENT 


FORTRAN SOURCE LIST NINLIN 


•1,G250 

ISN 

110 

112 

112 

116 

117 

120 

121 

122 

123 

125 

126 
127 
122 

133 

134 

135 

136 
1 40 

142 

143 

146 

147 

150 

151 

152 

153 

155 

156 

157 
162 

163 

164 

165 

166 
170 

172 

173 

174 

175 

176 
201 

203 

204 

205 
207 
210 
211 
212 
215 
217 
220 
221 
222 
224 


103 CONTINUE 

C*****T3 MARK COLUMNS HAVING ZEROS IN ABOVE MARKED LINES 
200 C*LL STEP2(Nf JP » KTT ,fuNA, MODE ) 

I F ( MODE. EQ . 0 3 GD TO 11 


1 3 = I C 
I E = I D + l 

DO 104 J = 1 s MODE 

i:=ic+i 

MR KC ( 1 C ) = N 1 N A ( J ) 


104 CONTI NUt 


DD 105 J = 1 f I D 
DO 105 J J = 1 u » IC 

I F { MRKC { J ) .cU.MRKCt JJ 3 3 GO TO 150 
GD TO 106 

iso j r = j j - 1 

NI NA( JT) =0 
MR KC ( J J ) =0 

105 CONTINUE 
105 CONTINUE 

C^ LL STi.. P 1 ( N » N I iNA i MCD. »JP,KTT) 

1= (KTT.EQ.05G0 TO 11 
I T = I R 
I J= I R+ 1 

DD 107 I = lj KIT 
I R = I R+ 1 ' 

MRKR( IR)=JP( I ) 

107 CONTINUE 

DD 109 11=1, KIT 
DD 103 1=1, IT 

IF (MRKR5I ) .EQ. JPUI 33GO TO 160 
GD TO 108 
160 IN = ITf II 

MRKR ( I M ) =0 
J?( I I ) =0 

108 CONTINUE 

109 CONTINUE 
GD TO 200 

11 UR =0 

DD 301 1 = 1 , N 

DD 302 J = 1 , 1R 

IF ( NRKR{ J) «EQ.I 3 GO TO 301 

302 CONTINUE 
UR =UR+ 1 
UNMRKR C UR 3 =1 

301 CONTINUE 

u:=o 

DD 303 J =1 , N 

DD 304 K = 1 , 1 C 

IF (MRKC(K) .EQ,J J GO TO 305 

304 CONTINUE 
GD TO 303 

305 u:=uc+i 

MRRKC(UC)=J 

303 CDNTINUE 
RETURN 


FORTRAN SOURCE LIST M I ML I 


i'1 > G 2 5 0 


ISN 

S3 URCl. S T A T _ M t. N * 


225 

M.-G250 

END 

IB MAP ASSEMBLY M I N L I 


NO MESSAGES FDR ABOVE ASSlMBLY 



M_ G25G 


FORTRAN SOURCE LIS T 

ISM 


SOURCE STAT, t MT 

0 

$1 8FTC 

STEP1 

1 


SJ BROUT I ME STtl 3 1 ( M , KR , KM* JP, KTT) 

2 


COMMON / JLM/P 

3 


COMMON / ZLM/D 

4* 


INTEGER D,P 

5 


DIMENSION .■( C, • 0 ) T PC 0, A-0 } , KR{ 1 00 ) » J P { 100 5 

6 


KT =0 

7 


00 9 J = 1 » N 

10 


00 10 IK*1, KM 

11 

ID 

CO NT I MUE 

13 


GO TO 9 

1A- 

11 

DO 12 1 = 1 » N 

15 


IF(P(I,J).ey.l)GO TO 13 

20 

12 

CONTI MUE 

22 


GO TO 9 

23 

13 

Kf =KT+ 1 

2b 


(KT) =1 

25 

9 

CONTINUE 

27 


KT T =KT 

30 


RETURN 

31 


END 

MLG250 


I8MAP ASSEMBLY STEP1 

MO MESSAGES 

FOR ABOVE ASSEMBLY 



M...G2 50 


FORTRAN SOURCE LIST 

ISM 


SOURCE STAVoHiiNT 

0 

$I8FTC 

, S7EP2 

1 


$ J BRQJT'I Mi, S' 1 P2(N, JP,KTT, NINA, MODE) 

2 


COMMON / JLM/P 

3 


COMMON / ZLM/D 

/# 


DIMENSION D { -,0 , >0 ) , P { -0 , 40 ) » J P ( 1 00 ) , N I N A ( 100 5 

' 5 


INTEGER D , V 

6 


KO =0 

7 


DO 14 1=1, N 

10 


DO 15 KT=1 , KIT 

11 


IF ( I.EQ.JP(KT) ) GO TO 16 

14 

15 

CONTINUE 

IS 


GO TO 14 

17 

IS 

DO 18 J = 1 , N 

20 


IF(DU, J).EU.O)GQ TO 17 

23 


GO TO 18 

24 

17 

IF{ P{I, J) .EQ.DGO TO 18 

27 


K0=KD+1 

30 


NINAUD) = J 

31 


IF (KD.fcQ.l )G0 TO IE 

34 


M_ A=KD-1 

35 


DO 30 11=1, MLA 

36 


IF (NINMKD).EQ.NINA (II ) )G0 TO 31 

41 

33 

CONTINUE 

43 


GO TO 18 

44 

31 

ML ZP=KD 

45 


NI NA ( MLZP ) =0 

46 


KD=KD- 1 


18 

CONTINUE 


14 

CONTINUE 

53 


MODE=K.D 

54 


RETURN 

55 


END 

McG250 


IBMAP ASSEMBLY STEP2 

MO MESSAGES 

FOR ABOVE ASSEMBLY 



FORTRAN SOURCE LIST 


MLG2 30 

. ISM SOURCE StAT M NT 


0 

$1 BFT 

C OR ZERO 

1 


SJBR3UTI NL CRZ RO (N »Ui MRKR , MRKC , UR ,UC) 

2 


COMMON / Z LM/D 

3 


INTEGER D»UMMRKR|UC »UR 

4 


01 MENS I DM D{<>0,'+0), UNHRKRt 40) ,MRKC{50) 

5 . 


KM =0 

6 


DO 1 1 = 1, N 

7 


DO 2 K = 1* UR ■ 

10 


IF (UNMRKRC K) .LQ.I )GC TO 1 

13 

2 

CONTINUE 

15 


KM =KM+1 

16 


K~ = 0 

17 


DO 4 J = 1 , N 

20 


DO 5 L=L»UC 

21 


IF { MRKC ( L ) «EQ* J ) GO TO 4 

24 

5 

CONTINUE 

26 


KL =KL+ 1 

27 


1 3 ( KM, KL ) = 0( 1 1 J ) 

30 

4 

CONTINUE 

32 

1 

CONTINUE 

34 


DO 7 1=1, KM 

35 


DO 8 J = 1 , KL 

36 


JB(J)=I8(I,J) 

37 

8 

CONTINUE 

41 


CALL MINI JB,KL,KLN,MNL) 

42 


J 3 J ( I ) =MNL 

43 

7 

CONTINUE 

45 


CALL MIN(J8J»KM,KMN,MNM) 

46 


DO 9 I =1 , N 

47 


DO 10 K = 1 j UR 

50 


IF (UNMRKR(K) .EQ.I )GC TO 11 

53 

13 

CONTINUE 

55 


DO 12 J = 1 » N 

56 


DO 13 L = l, UC 

57 


IF ( MRKC { L ) . EQ. J ) GO TO 14 

62 

13 

CONTINUE 

64 


IFIDU, J).EQ. 9999)00 TO 14 

67 


D( I , J ) =D ( I » J ) ~MNM 

70 


GO TO 12 

71 

14 

D( I, J) =D ( I , J ) 

72. 

12 

CONTINUE 

74 


GO TO 9 

75 

11 

DO 15 0=1, N 

76 


DO 21 L= 1, UC 

77 


IF {MRKCI L) .EQ.J )G0 TO 16 

102 

21 

CONTINUE 

104 


D( I,J)=D(I , J) 

105 


GO TO 15 

106 

16 

IF ( D( I , J ) * EQb 9999 ) GO TO 15 

111 


D( I,J)=D(I, J)+MNM 

112 

15 

CONTINUE 

114 

9 

CONTINUE 

116 


RETURN 

117 


END 


MEG 2 50 ^ , ' IB MAP ASSEMBLY CRZIftO 

NO MESSAGES FOR ABOVE ASSEMBLY 


D 0 ) f J b J { 



K : G250 


FORTF 

.AN SOURCE LIST 

I Si-4 


SOURCE STAT.'i 47 


0 

SI 3FTC 

MIN 


1 


SJ BRjJTI NL M ! U < P, N , K, KMI N) 


2 


INTEGER P* KMl.i» SUM 


3 


DIMENSION PCS.) 


4 


K= 1 


5 


S J M=P ( 1) 


6 


IF { N. EQ. 1 ) GO TO 2 


11 


DO 1 I = 2 j N 


12 


I F { P ( I ) . GE • SUM ) GO TO ' 1 


15 


K= I 


16 


S J M=P ( I ) 


17 

1 

CONTINUE 


21 

2 

RM IN=SUM 


22 


RETURN 


23 


END 


M2G250 


IBMAP 

ASSEMBLY MIN 

M Q MESSAGES 

FOR ABOVE ASSEMBLY 


M. G2 50 


IBLDR 

— JOB OOOOOO 


*** 

OBJECT PROGRAM IS BEING ENTERED 

INTO STORAGE AT 11 HRS. 22 MTS 



*•**•****•*«»#*#•#*** «««*#••#*••»•#****«*#*«**•** *«•*** 
* * 

*' PROGRAM FDR ADJUSTING INITIAL ROUTFS FOR CAPACITY * 

* AND TIMS CONSTRAINTS AND FCR REFINING THE ADJUSTED* 

*■ ROUTES. PROGRAM ALSO CONSIDERS OMITTED CUSTOMERS . * 

* * 
■********#***■*«#•» ■#*"*#***•***#■***#■*•****•**#*:*«****«**■*■*■**■* 


C*****THI5 PROGRAM ADJUSTS THE INITIAL ROUTS TC BE WITHIN CAPACITY 
C AND TIME CONSTRAINT FOR LEAST INCREASE IN DISTANCE 

C AND THEN IMPROVES THE ROUTES BY REFINING HfEURI STI3S 

C*** THIS PROGRAM IS APPLICABLE TO BOTH TYPE OF VEHICLES 
C*****N( I ) DENOTES NO. OF NOD' S (EXCLUDING DEPOT IN SUBTOUR I 
C*****TOTAL NO. OF NODES ARE NC. OF CUSTOMERS PLUS NO. OF ARTIFICIAL 
C DEPOTS WHICH ARE EQUAL T C NO. OF VEHICLES 

C*****B(I, J) IS DISTANCE MATRIX SHCWING DISTANCES BETWEEN ALL NODES 
C*****D( . ) DENOTES DEMANDS OF ALL NODES. DEMANDS FOR ARTIFICIAL DEPOTS 
C ARE ZEROtTHAT MEANS D (1 ) , D ( 2 ) ,0 ( 3 ) ETC. ARE ZFRO. 

C*****KRT(I, J) DENOTES JTH CUSTOMER IN ITH ROUTE FOR EXAMPLE IN SUB* 
C TOUR 1—4“ 6-5—1 » KRT( 1 * a > =4, KRT (1 ,2 ) =6 FTC. 

C*****INITIAL TOURS OBTAINED IS, SAY, IN THE FORM- 1-4-6-5*3-10- U- 1 THI 
C SUBTOURS ARE X-4-6-5-1 AND 3-1 C-ll-3 PUT DATA IN THIS FORM 

C IS EQUAL TO 1 THEN ROUT" TRAVERSES FROM I TC J NODE 

DIMENSION NU 0) ,MMX (*0,40), KRTf 10,30) ,D<40) ,B(40, 40) 

DIMENSION JMX (40,40 ,OVCt 10) , UVC (10 ) ,KAPO Cl 0) , KE F ( 1 0 ) ,ODVC ( 10 » 
DIMENSION NBK J (4 0,40 ) , JCB(40) 

DIMENSION KTO VC(10),C(10) 

'DIMENSION IMEK(IO) 

DIMENSION TCR ( 10 ) , JRTD ( J O) , JRTd C) , JVT ( 10 ) 

INTEGER 8, C ,T ,D , OVC, UVC , COVC 

DIMENSION J ST (40 ) » JAV ( 10 ) »CT( 10) ,FX (10 ) , JS£T( 10) , JAL ( 10 ) , CST( 10 J 
INTEGER TD 

COMMON/OPM/B, T , JMX»KR T 
COMMON/ JLM/D, C 
COMMGN/ STS/NV ,NN 
. 3081 FORMAT ( 213 ) 

READ 3081, (JAV(I), 1=1,2) 

READ 3082, ICT(I), 1=1, 2> 

3082 FORMAT C 2F6.1 ) 

READ 3083, ( FX ( II, 1=1, 2) 

3083 FORMAT ( 2F8 . I ) 

320 READ 101, NC 
101 FORMAT ( 131 

IF(NC.EQ. 100)60 TO 7010 

NCC=NC+i 

DO 1580 1=1, NC 

JJ=I+i 

READ 1582, ( NBK.J < I » II ) ,1 1 = J J *NCC ) 

1580 CONTINUE 


1582 FORMAT ( 2014 ) 

902 FORMAT! 1H1, 5X ,*SET POLICY*, I3/10X,*SET NUMBER*, I 3/10X,*VEHICL€S TO 
6 BE USED*, 1013) 

903 FOR MAT (/I OX,* NO. OF CUSTOMERS =*,I3) 

■DO 103 1=1, me 

DO 103 J=l, NCC 
IF( J.GT.UGC TO 103 
IF ( J.EQ. I ) GO TO 1588 
NBKJ(I,J)=NBKJ(J,I ) 

GO TC 103 

1588 N8K J! I » J > =9999 


a. a. cl lOEr 



103 CONTINUE PRO 

RE AC X57], ! JST(I ),I=? ,NCC) PRO 

1571 FORMAT 12014 ) PRO 


READ 340,1 JD81J),J=1,NCC) 
740 FORMAT (2014) 

321 READ 7123, NPQLCY.MIPCY, CODE 
7123 FORMAT ( 213 * F6 .1 ) 

IFINPOLCY.EC. lOOGO TO ?20 
950 READ 910, NV 
910 FORMAT ( 13 ) 

IF(NV.EQ.500)G0 TO 321 
PE AC 3088, (JSETU), 1 = 1, NV) 
NV V=NV+X 
NN=NV+NC 
CO 351 1=1, NV 
DO 351 J=1 , NN 
IF ( J.GT.NVJGQ TO 356 
8 ( I, J ).=9999 
GO TO 351 
356 J J = J—NV+1 

8(1, J)=NBKJ!i , JJ ) 

351 CONTINUE 

CO 352 I=NVV, NN 
I 1=1— NV+1 
DO 352 J=1 ,NN • 

IF ( J.GT .NV )G0 TO 355 
8! I, J) = N8KJl I 1,1 ) 

GO TO 352 
355 JJ= J—NV+1 

B < I » J ) =N8K J { I I , J J ) 

352 CONTINUE 

CO 9100 1=1, NV 
0VC(I)=0 
9100 UVC( I )=0 
365 CONTINUE 
367 DO 375 I=1,NN 

IF ! I » LE. NV ) GO TO 956 


PRO 

PRO 

PRO 

PRO 

PRO 

PRO 

PRO 

PRO 

PRO 

PRO 

PRC 

PRO 

PRO 

PRC 

PRC 

PRC 

PRC 

PRC 

PRO 

PRO 

PRO 

PRtj 

PRtj 

PRO 

PRC 

PRC 

PR(! 

PRI 

PRC 

PRI 

PR] 

PRI 

PR' 

PR 


I I=I-NV+1 
C< I) = JDB{ II ) 

GO TO 375 
956 D ( I ) =0 

375 CONTINUE 
377 CONTINUE 

RVP=BVP/100. 

DVP=RVP+4.0 
CO 7121 1=1, NN 
XLC=FLOAT! D(I )3 
D( I) = IFIXl DVP*XLD) 

PRINT 7122, BVP, (0(1), 1=1, NN) 
FORMAT i 1H1 »1QX, *DE HANDS AT EACH 
6/5X,20I5/5X,20I5 ) 


1121 

7122 

380 

102 

397 

1,592 

393 


READ 102, (NCI ) , I =1 ,NV ) 

FORMAT 11014) 

FORMAT C//5X,* INITIAL TOURS BY MULTIPLE TRAVELLING SALESMEN 
FORMAT! //20X, *DEMAND5 ARE AS FGLL0WS*/5X,20I5/5X,2015 » 
FORMAT! /8X,*SUBTOUR*,5X,*NQ. OF NODE S*»SX,*PATH SEQUENCE*/ ) 
CO 105 1=1, NV 
NPL=N(I) 

READ 106, (KRT !I,J),J=1,NPL) 


PR I 
PR 
PR 
PR 
PR 
PR, 
PRi 
PRj 
: P 

P 
P 
P 

NODE INCREASED BY*,F6.1,*P£RC€NT*, P 

P 

P 

P 

PROB*/)PI? 

pi 

p ! 

pij 

pi 

pc 


305 CONTINUE 
106 FORMAT (2513} 

391 FORMAT ( iOX, I? » 10X , 13 * 10X . 2513 ) 

C********* CODF=0» FOR CONSTANT CAPACITY VEHICLE 
C C ODE = 1 « FOR DIFFERENT SIZED VEHICLES 

JP0LCY=0 

8868 JPCLCY = J POLCY + 1 

IF (CODE. HQ. C. )G0 TO 9227 
IF ( JP0LCY.GT.3JG0 TO 95 C 
PRINT 902, NP(1LCY,MIPCY, ( JSETC I ), 1=1, NV) 

PRINT 34'!, N V 

399 FORMAT l SX* *CONS IDE RING VSHICLE-FCLICY*, 13) 

GO TO 922.8 

9227 IF ( JPCLCY.GT.l )G0 TO 95C 

PRINT 902,NPOLCY,MIPCY, (JSETII ) ,1 = 1, NV) 

PRINT 342, NV 

9228 RFAC 100, T, (C ( I) , 1=1, NV) 

100 FORMAT (1515) 

READ 3085, ( JAL( I ),! = !., NV) 

3085 FORMAT (1013) 

3088 FORMAT (1013) 

8890 FORMAT (5 X,*CA PAC I TIES OF VEHICLES*, 1 C 15 ) 

I F ( JPOLCY-GT. 1 >G0 TO 90’ 

DO 200 1=1, NN 
DO 200 J=1 » NN 
MMX ( I, J 1=0 

200 CONTINUE 

DO 201 1=1, NV 
MPL=N( I ) 

JS = I 

DO 202 J = 1 , MR L 
KRTT=KRT( I, J) 

MMX ( JS , KRTT ) = 1 
JS=KRTT 
202 CONTINUE 

MMX ( JS, I ) =1 

201 CONTINUE 
901 IPCLCY=0 

PRINT 399, J POLCY 

PRINT 8890, (C(I), 1=1, NV) 

900 DO 110 1=1, NN 

do no j=i , m 

JMXU, J)=MMX( I, J ) 

110 CONTINUE 

PRINT 8892 

8892 FORMAT ( 5X, 33( * -*) ) 

IPOLCY=IPOLCY+1 

IF ( I FOLCY. GT. 3 )G0 TO 7000 

C*****POLICY FOR CONSIDERING OVERLOADED ROUTES IN RANDOM ORDER 
PRINT 851,1 POLCY 

851 FORMAT* 5X, *R0UTE-P0LICY#, 13 1 
MOT 1=0 

CALL GUCIYA(JMX,KRT,N) 

1161 KEJ=C 
JE J=0 

1001 CALL CAPCONCKRT »N»OVC,OVC) 

IF ( IPOLCY. EC. i ) GO TO 1002 
GO TO 1100 


PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PRC 

PR! 

PR! 

PR! 

PR1 

PR! 


PR 

PR 

PR 


PR 

PR 

PR 


PR 

PR 

PF 



IF (OVC(I ). r C. G) GO TO 119 
JE J = JE J + i 
KAPO ( JEJ ) = I 

CALL SEARCH IN ,QVC , UVC , KAPC, JEJ,KFP,KEJ) 

CALL CAPCONf KRT,N,OVC,UVC ) 

IFaPCLCY.ee. 1)60 TO ICO? 

IF!IPCLCY.EC.2)G0 TO 12C0 
IF!IP0LCY- C C.3)GQ TO 13C0 
119 COMINUE 
6000 CONTINUE 

CALL CISCAL(TD) 

811 FORMAT ! /‘I X » ®T GTAL DISTANCE =*» 15 ) 

IFCKEJ.EQ.OGC TO 800 
DO 960 1=1 , KE J 
KAMG=KEP! I ) 

CO 961 11 = 1, NV 

IF (UVC(II ) .FQ.CIGO TO 982 

IMEK(II )=D!KAMO)~UVC!II ) 

C-0 TC 961 

982 IMEKIII l=5CC0 
961 CONTINUE 

CALL MIN! IMEK ,NV*NVH,I MK) 

IF! IMK.EQ.5COOGO TO 979 
KRZP=KRT( NVH, l) 

JMX { NVH,KAMO) =1 
JNX(NVH,KRZP) =0 
JMX ( KAMO* KRZP )=1 
CALL GUCIYA! JMX,KRT,N) 

CALL CAPCON ! K RT *N*OVC»UVC) 

960 CONTINUE 

M0TI=M0TI+1 

IFCM0TI.FQ.21 )G0 TO 979 
GO TO 1161 
979 PRINT 973 

973 FORMAT! /10X»* OMM USD CUSTOMER REMAINS*) 

GO TC 900 
800 CONTINUE 

850 FORMAT! //10X,*THERE IS NO CUSTOMERS LEFT UNASSIGNED*) 
CALL REFINE(N ) 

PRINT 1350 

1350 FORMAT !10X,*F INAL ROUTES OBTAINED ARE*) 

PRINT 1360 

PRINT A 6/. 8 
NVN=NV-1 
DO 1351 1=1, NV 
NE P=N ! I ) 

CO 1352 J=1»MEP 
KR IM=KRT II , Jl 
KRTCI, J)=KRIM*»NVN 
1352 COMINUE 

1351 CONTINUE 

DO 3040 1=1, NV 
LI ZI = N ( I ) 

NST=0 
N I D=0 
JRTK=0 
JSS = 1 

DO 3041 J=1 , L IZI 
KLIZ=KRT! I,JI 


PROQ132 1 

PR001329 

PR 001 331. 

PP001332 

PR00J351 

PR001352 

PR001353 

PROG 115^ 

PROOI 355 

P R 0 01 3 5 6 

PR001357 

PR001358 

PR0G1359! 

PR00136G 

PR001361 

PR001362 

PR0C1363 

PR0C1364 

PR001365 

PR 001 366 

PR001367 

PR001368 

PR001369 

PR00137C 

PR001371 

PR001372 

PROOX373 

PR001379 

PR0C1381 
PROQ1382 
PR00138: 
PR0G1384 
PR00138S 
PRO0138< 
PR001381 
PRO 013 8? 
PRO 013 8* 
PR00139! 
PR001393 
PROOI • 9: 
PP00139< 
PR00139! 

PR0013SI 

PR00139 J 
PROOI 39 
PRO 01 7 9 
PR0014C 
PR00140! 
PROOI 4cf 

, : ' ' ■ f 



NST=N$T+JST(KLIZ ) 

NIC = NID+JDB(KHZ ) 

JRTK=JRTK+N8KJ(JS5,KLIZ ) 

JSS=KLIZ 
3041 CONTINUE 

JRTD(I )=JRTK+NBKJ( JSS, 1 > 

RTD=FL0AT« JRTDtl ) )/4.0 
UVC( I )=C( I }*»N ID 
IF ( C ( I )■• ? Q .40 CC ) GO TH 2059 
CST( I ) =CT ( 2)*RT0 
TCR(I)=FX(2)+CST (I) 

JRT(I ) = (;’ 5*JRTD( I ) / JA V ( 2 ) J + NST 
JVT ( I )=JRT ( I > + (2*JSTU ) t 
GO TC 30? 1 

3059 CST ( I ) =CT ( 1 )«■ RTD 

TCR ( I ) = FX( 1 J+CST ( I ) 

JRT! I ) = ( 4 5*JRTD( I}/JAV( J ) + NST 
JVT ( I ) = JRT ( I ) -M2*JST(l ) > 

3051 I 1=1 

PRINT 1361 » I, JAL (I },C(I !,UVC( I), JRT(I),JVTd I ,11, 

6 ( KRT ( I » J ) , J =? ,L I ZI ) , J I 
B 040 CONTINUE 

TC = I FIX (RTD ) 

C S T R = 0. 

TCRT=0. 

DO 3029 1 = 1, MV 
CSTR=CSTR+C$T ( I ) 

TCRT=TCRT+TCR (I) 

3029 CONTINUE 

PRINT 8 lit TD 

PRINT 3033, CSTR,T CRT 

3033 FORMAT(5X,*COST OF TRANSPORTATION**, F9.2/5X»*T0TAL C0ST=*,F9 
GO TO 900 

11C0 IF(IPCLCY.EG.2)G0 TO 12C0 
GO TO 1300 
1200 DO 1110 1*1, NV 
OOVCl I J=~OVC( I ) 
tllQ CONTINUE 

1360 FORMAT (5X,*R0UTE*,2X,*VEH- ALLOC .*,3 X, *UVC*, 3X,* RT*,3X, *VT*, 
6UTE PATH* } 

1361 FORMAT ( 6X, 12, 2X, I4,4X ,1 4,3X, 14, IX ,14 ,1 X,I4 ,5X,i 0 13 1 

4618 F0RMAT(12X,*NC. CAP-* ,7X, * ( MTS . )* » IX, *( MTS > *> 

342 FORMAT! 10X»*NC- OF VEHICLES =*,13} , / 

1140 CALL MINtOOVC ,NV, NVO, MOVC) 

NOVC=->MOVC 

IF (NOVC.fQ. C) GO TO 6000 

JFJ=JEJ+1 

KAPCl JEJ)=NVO 

5000 CALL SEARCH (M , OVC ,UVC »KAPO,JFJ,KEP,KEJ} 

CALL CAPCONCKRT.N, OVC, UVCl 
GO TC 1100 
13 CO CO 1210 1=1, NV 
KTQVC ( I 7=0 VC ( I ) 

IF(OVC( I)-EQ-0)KT0VC(I } =999999 

1210 CONTINUE 

CALL MlN(KTOVC,MV»NVF » JCVC1 
IF (JO VC. FQ -999 9 9 9) GO TO 6000 

IF ( JC VC.FQ- 99 59 )G0 TO 6COO 
JEJ=JEJ+J 


procia-c: 
ppnoi4C< 
PFO0I40; 
PR0C14 C> 
PROOI4C 
PR00340; 
PROCl^c* 
P ROC 14 1 1 
PR 00 14 1 i 
' PR0014J! 

p roc i*i : 
PROG 34 !> 
PR00141! 
PR0014 li 
PR 00 14J‘ 
PR 00 141t 
PR0C14?.: 
PR00142! 
PR00142 
PR00142; 
PR0C1 2 
PR00342 
PR00142 
PROG? 42 
PR00142 
PRO 014 
PR001.4; ' 
PR00I4; 
PR0C14 
PR0014: 
PR0014 
.2) PR0034 

PR00143 
PROC: 42j 
PR00143 
PR001 4^ 
PROOX^B 
PR0C144 
5X,*R0PR00144 

PR0C144 
PR0014^ 
PR0C144 
. PR0014^ 


PR0016 


PR0016 


K/sPO C JF.J)*NVE 
GO TO 5000 
7000 GO TO 8888 
7010 STOP 
ENC 
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STAT: Ml: NT 


FORTRAN SOURCc LI S'!' 


.G 250 


ISM 


S 3 URC; 


$1 8FTC SEARCH 

SJBROUTIME S ARCH{ N »QVC » UVC, KAPO, J EJ » KEP ? KX J ) 

INTEGER B » D» OVC , UVC ,C, T , BB , A , QB J FN 
C3MM0M / STS /MV j N I 
COMMON / JLM/Dj C 
C3MMDM/DPM/B, T, JMX» KRV 

DIMENSION KRT ( 10 > 30 ), ■!{ 10 ) ,0VC{ 10) ,UVC{ 10) 

01 MENS I ON KAPO ( 10) , KEP (10 S » J MX( 40 » 40 5 

DIMENSION B { 0, ; 0 } , D{ 0 ) » J EST ( *:-0 ) , MINT < A-0 ) , i SC ( M; ) , I CS ( 40 ) 

01 MEM SI 3M MINMJM(40),i VI ( 4G) , 1V2 (4 0) , HINDIS* 43 ) , I w It *0 5 , Iw2( 40) 
Df MENS I ON IT? „ M ( 10, 30 ) , IB { 10 ) , At 10,40 5 , BB ( 10 ) , OB JFN (AC ), NX ! - X ) 

dimension. ndis(40),'nx:;< 40) ,c ( ioj 

M3 M=KA PQ { J <- J ) 

NV S=N( MGM) 

D3 2 J=1 , NVS 
K_ =0 

KY = KRT ( MGM , J ) 

03 3 J J = 1 » UN 
D3 1 JP=1,JFJ 
M3S = KAP 3 { JP 5 

i=(duy).le.uv:(mgs))go to i 

IF ( JJ. EQ.MGS )G3 TD 3 
NV F = N( MGS) 

D3 101 JPS=1,NVF 
KZY=KRT(MGS»JPS) 

IF ( JJ.lU.KZY)GO TO 3 
101 C3NTINUE 
1 CONTINUE 

1= (KEJ.EQ.O)GO TO 1520 
D3 1510 *JAJC*1,KSJ 
MGZK=K6P( JAJC) 

IF { J J. 6Q.MGZK ) GO TO 3 
1510 CONTINUE 
1520 ID IS = 3 ( J J* KY ) 

K.3 =0 

D3 5 JK= 1 » NN 

1= ( JMK( JJ. JK).EQ.0)G0 TO 5 
JDI S = 3 ( KY* JK ) 

M3 I S=I 01 S+JD1 S 
K3IS=B( JJ, JK) 

KQ=KQ+ 1 

N3IS(KQ) =MDIS~<DIS 
JESTtKO) =JK 
5 C3NTINUE 

CALL MIN(NDIS,KQ,LSD,IMINA) 

K_=KL+1 

HINT(KL)=IMINA 

ISCIKL 

ICS(KL )=JCST(LS0) 

3 CONTINUE 

IF (KL.EQ.O)GO TO 999 
995 CALL M I N ( M 1 NT , K L , L S T , 1 M I N ) 

IF ( IMIN.EQ.9999 )G0 TO 999 
MI MMUM ( J ) - 1 MIN 
I V 1 ( J ) =1 SC ( LSI) 



M...G2S0 


FORTRAN SOURCE LIST S .ARCH 

ISM 


SOURCE ST ATE Mu NT 

114 


IV2(J)=ICS(LST) 

115 


IF { J.EQ. 1) GO TO 998 

120 


JT A = J- 1 

121 


03 1082 JQA = 1 , J T4 

122 


1= (IVKJI.LQ.miJQ'AI )G0 TO 997 

125 

1032 

CONTINUE. 

127 


GO TO 998 

1B0 

997 

MI NT( LST ) =9999 

131 


GO TO 996 

132 

998 

KSK E=0 

133 


IF { JMX(MGM,KY). EQ.O )Gu TO 9 

136 


LKY=MGM 

137 


K- KE = K EKE + B ( LKY , KY ) 

140 


GO TO 8 

141 

9 

DO 6 JJ=1, NVS 

142 


L< Y=KRT { MGM » -J J ) 

143 


IF { JMX ( LKY* KY) .,,Q.O)Gu TO 6 

146 


KEKE=K£KE+B(LKY»KY) 

147 


GO TO S 

150 

5 

CONTINUE 

152 

8 

1= (JMX{KYtMGM),EQ.O)GO TO 10 

155 


M< Y=MSM 

156 


K5KE=KEKE+B(KY, MKY ) 

157 


GO TO 11 

160 

10 

DO 12 J J = 1 , NVS 

161 


MKY=KRT( MGM, JJ ) 

162 


IF { JMKiKY, MKY). EQ.O )Gu TO 12 

165 


KEKE=KEKE + B( KY, MKY ) 

166 


GO TO 11 

167 

12 

CONTINUE 

171 

11 

K4 KU=KEKE-B( LKY, MKY ) 

172 


MI NDIS ( J ) = IMIM-KAKIJ 

173 


2W 1 ( J ) =LK Y 

174 


14 2 ( J ) =MKY 

175 


GO TO 2 

176 

999 

IMIN=9999 

177 


NHJ) =0 

200 


IV2 (J) =0 

201 


GO TO 998 

202 

2 

CONTINUE 

204 


KIINN=l 

205 


IF { JEJ.EQ. 1)G0 TO 20 

210 


JET=JE J* 1 

211 


00 400 1=1, JET 

212 


KI N*0 

213 


M3ST=KAP0( 1) 

214 


INVFS=N{ HOST) 

215 


DO 500 J = 1 ,NVS 

216 


KZW=KRT ( MGM* J ) 

217 


DO 600 J J = 1 f NVFS 

220 


I F (, I VI ( J ) • LQ. KRT ( MG ST , J J ) ) GO TO TOO 

223 

600 

CONTINUE 

225 


DO 750 JJ=1»-.;VFS 

226 


IF(lV2(J).[-Q.KRT(MGS7 s JJ))GO TO 700 

231 

750 

CONTINUE 



M„; G2 50 FORTRAN SOURCE LIST S 


ISN 


SOURCE STATEMENT 

233 


GO TO 450 

234 

703 

KI N=KI N + l 

235 


ITEENt KINN, J ) 0{ KZW ) 

236 

453 

I T E E N ( KINN* J S=3 

237 

503 

CONTINUE 

241 


I _ ( K I N « L" 0 • U ) GO TO 400 

244 


B { KINN)=UVC(rt3ST) 

245 


KI NN=K INN+ 1 

246 

400 

CONTINUE 

250 

20 

n:ons=kinm 

251 


NV AR = N VS 

252 


DO 2203 J= 1 , NVS 

253 


KRTT=KRT(MGM, J) 

254 


A( 1* J) =D ( KRTT ) 

255 

2203 

CONTINUE 

257 


B3 ( 1)=DVC( MGM) 

260 


1= (NCDNS.EQ. 1 ) GO TO 2. 00 

263 


DO 2204 I = 2*1 ICONS 

264 


Jl N = I-1 

265 


D3 2206 J = l» NVAR 

266 


B3 (I )=IB ( JIN) 

267 


A( I ,J) = ITtiEN( JIN, J ) 

270 

2206 

CONTINUE 

272 

2204 

CONTINUE 

274 

2500 

DO 2202 J=1,NV4R 

275 


08 J FN( J ) =MI NDI S ( J) 

276 

2202 

CONTINUE 

300 


DO 840 1=1, NVS 

301 

840 

NX (I )*0 

303 


C4LL ZfcRDltNVAR, NOONS , A, Bti ,0B JFN, NXX, I AB ) 

304 


DO 825 1=1, 1AB 

305 


NXT=NXX( I ) 

306 


NX t I'XT ) = 1 

307 

8 25 

CONTINUE 

311 


DO 6000 J = 1, NVS 

312 


KR I A=KRT ( MGM, J ) 

313 


IF { NX ( J ) . EQ* 1 J GO TO 6001 

316 


GO TO 6000 

317 

6001 

IF (HINDIS! J) .GT. 900 0 GO TO 6002 

322 


IV V 1 = 1 VI ( J ) 

323 


IVV2 = IV2( J > 

324 


J4XUVV1, IVV2)=0 

325 


JNX(IVV1,KRTA)=1 

326 


J *1 X ( KR T A , I V V 2 ) = 1 

327 


GO TO 6003 

330 

6002 

KE j=ke J + l 

331 


KS P (KEJ) =KRTA 

332 


IV Vl = 3 

333 


I VV2=0 

334 

■ 6003 

UW1 = I Wl( J) 

335 


IrfW2=lW2( J) 

336 


1F(JMX(IWW1, KRT A ) « . Q. 1 ) GO TO 6050 

341 


1= ( JMX ( MGM»KRTA ).EQ.1)GQ TO 6069 

344 


DO 606 0 J J = 1 , NV S 

345 


KRAC=KRT( MGM, JJ ) 


RCH 



G250 



FORTRAN SOURCE LIST SEARCH 

ISM 


SOURCE STATEMENT 


3 %6 


1= ( JMX ( KRAC » KRTA ) .2 C. 1) GO 

TO 6070 

351 

6060 

CONTINUE 



5069 

KR AC=MGM 


35 % 

5070 

IW W1=<RAC 


355 

5 05 C 

JMX(IWW1»KRTA) =0 


3 56 


JMX ( KRTA » 1 WW2 ) =0 


35 7 


JMX(IWW1»IWW2)=1 


360 

6000 

CONTINUE 


E 62 


CALL GUDIYAI JMX,KRT»N) 


363 


RETURN 


3 64 


END 


.G250 



IBMAP ASSEMBLY St ARCH 

NO MESSAGES 

FOR ABOVE ASSEMBLY 




FORTRAN SOURCE LIST 


M . G25D 

ISN S3 URGE STATEMENT 


0 

1 

2 

3 

4 

5 

6 
7 

10 

11 

12 

13 

14 

15 

16 
20 

23 

24 

25 

26 

30 

31 

MEG250 


SI 8FTC 3APC0N 

SUBROUTINE 3 ARSON ( KRT, N, OVC ,U VC) 

C3MM0N/STS/NV,NN 

C3MM0N/JLM/D,C 

INTEGER. D,nvc, UVC,C,SUM 

DI MENS1DN KRT( 10,30), D{40> ,SUM{ 10) ,QVC ( 10 ) , UVC ( 10) , ill 10 ) 

D3 1 1=1, NV 

OVC (I)=0 

UVC(I)=0 

SJM(I) =0 

NV S=N( I ) 

D3 2 J = 1 , NVS 
KRTT =KRT ( I , J ) 

S J MCI ) =SUM( 1 ) +D { KRTT ) 

2 CONTINUE 

IF ( SUM ( I ).U .3(1) ) GC TO 3 
□V C ( I ) =SUM ( I )~C(I) 

GO TO 1 

3 UVC(I)=CU ) * SUM ( I ) 

1 CONTINUE 

RETURN 

END 

IBMAP ASSEMBLY CAPCDM 


NO MESSAGES F3R ABOVE ASSEMBLY 



Uj fNJ -*4 CT' VJI N) 1/" Ui 4> iv* f\> h~* o 


FORTRAN SOURCE. LIST 


M i . G 2 5 0 

I SM 


MLG2 50 


S3URCL S TAIL Me NT 


SIBFTC MIN 

SJ8R3 JT I NK MIM{ P,N, K f RMIN) 

INTEGER P t RMIN, SUM 
DIMENSION P(100) 

K= I 

SJM=P( 15 

IF (N.EQ.15G0 TO 2 
D3 1 I = 2 , N 

IF ( P( I ) • GE*SUM ) GO TO 1 
K= I 

SJM=P( 1 5 

1 CONTINUE 

2 RM I N = S UM 
RETURN 
FN D 

IBMAP ASSEMBLY MIN 


NO MESSAGES F3R ABOVE ASSEMBLY 



M..G2 50 


FORTRAN SOURCE LIST 

ISN 


SOURCc STATLflcMT 

0 

$1 BFTC 

. OUDIYA 

I 


SUBROUTINE GUOI YA( J MX i KRT , N) 

2 


COMMQN/STS/NVfNN 

3 


DIMENSION JMXI40, 40 ),KRT( 10,30), N( 10)' 

4 


DO 1 1=1, NV 

5 


K= 0 

6 


DO 2 J =1, NN 

7 


1= ( JMX ( I , J ) . Q. I )G0 TO 3 

12 

2 

CONTINUE 

14 

3 

U2 = J 

15 

5 

Jl = J2 

16 


K= K + 1 

17 


KRT(I,K)=J1 

20 


DO 130 J J= I, UN 

21 


IF { JMX{ Jl, JJ) . E Q . 1 ) GO TO 4 

24 

130 

CONTINUE 

26 

4 

J2 = J J 

27 


IF ( J2. EQ.I )GQ TO 6 

32 


GO TO 5 

33 

5 

m i )=k 

34 

1 

CONTINUE 

36 


RETURN 

37 


END 

MEG25Q 


IBMAP ASSEMBLY GUDIYA 

NO MESSAGES 

FOR ABOVE ASSEMBLY 



4. G25C 


FORTRAN SOURCE LIST 

ISN 


SOURCE STATEMENT 

0 

$1 BFT 

C l ERQl 

1 


SJBROJTINF L RO 1 { MM ,MM , A , B» C , JOPT j XAB } 

* 2 


INTEGER BMQST 

3 


INTEGER Y, A, B,C ,V, SGRI 

4 


INTEGER TEMP, ZMIN, C JSUM, ZOPT » SUMCR ,APOS , YASUM, CF SUH 

5 


DI HENS I ON J ST 4 0 ) » Y { 10} jJSUL(40! ,NS <40 ) » l-1~( 40) , V(--Q ) , MF { 40 ) 

6 


DIMENSION JP(a.D ) , JR ( 10 , 40 5 , NR < 10 ) , SORT ( 40 ) , JQPT( t u) 

7 


DI MENS I ON IN..Q( 10) , A( 10, 40) , B ( 1 0) , C C 40) , I NVC ( 40 ) 

10 


MS MAX= 100 

11 


MM MX = 3 

12 


INEQtl )=1 

13 


1= (MM. L:Q. 1 ) GO 7 0 5010 

16 


DO 5000 1=2, MM 

17 


IMEQ(I )=2 

20 

5000 

CONTINUE 

22 

5010 

DO 43 1=1,40 

23 


JP(I)=0 

24- 


JS < I ) = 0 

25 


JSULU ) =0 

26 


SORT (I )=0.0 

27 


DO 42 J = l, 10 

30 

42 

JR < J, I ) =0 

32 

43 

IN VC ( I ) =0 


C*4* 

READ PROBLEM, WRITE INITIAL TABLEAU 

3% 


CALL' G ERO ( NN , MM , MMM X, I NEQ , A, B ,C , INVC, MSMA X ) 

35 


DO 49 1=1, MM 

36 


NR (I) =0 

37 


DO 45 J=1 , NN 

40 


IF { A ( I » J ) ) 45, 44, 44 

41 

44 

NR (I) = NR ( I ) +1 

42 


IABC = NR ( I ) 

43 


1= ( A( I , J ) ) 4400,4400,4401 

44 

4400 

SORT ( I ABC ) = 7777777 

45 


GO TO 4402 

46 

4401 

SORTIIABC) =C ( J ) / A ( I, J) 

47 

4402 

JR ( 1 1 1 ABC ) = J 

50 

45 

CONTINUE 

52 


IF ( NR ( I )~ 1 ) 49, 49 , *6 

53 

46 

I A B=NR ( I ) - 1 

54 


DO 48 K= 1 , IAB 

55 


K? 1=K+ 1 

56 


I A BC=NR { I ) 

57 


DO 48 L = KP 1 , 1 ABC 

60 


IF (SORT ( K )-"SORT( L ) ) 40, 48, 47 

61 

47 

TEMP=SORT { K) 

62 


SORT ( < ) = SORT t L) 

63 


SORTt L ) =TEMP 

64 


ITEMP= JR < I ,K ) 

65 


JR C I , < ) = JR ( I ,L ) 

66 


JR ( I , L) = IT EMP 

■' 67 

43 

CONTINUE 

72 

49 

CONTINUE 

74 


ZmtN=77777?7 

.75 

i itii , 

MS = 1 ■ 

76 

300' 

BNOST = B ( MM ) 



FORTRAN SOUR CO LIST Z.:RU1 


i y L.G2 50 

ISM SOURCE STAT.-.-KMi' 

7T IMAX=1 

100 D3 30 A 1 = 1 , Mu 

101 IF ( BMliST- B { J ) ) ! 32, . CK, SO*- 

102 302 BM OST = b ( I ) 

103 1M AX=I 

10% 30% C3 NT1 ''i Ur. 

106 3 05 IA BC='4k( I ) 

107 03 325 I = 1,I,.BC 

110 JSC I)* JR (I MAX, I ) 

111 D3 319 K = 1 * M i'1 

112 Y{ K ) =- B ( K. ) 

113 D3 313 L=1 , I 

11% 1 A 8 = JS ( L ) 

115 313 Y(K)=Y(K)+A(K,IAB) 

117 1= (Y( K ) ) 32 , 319,3 19 

120 319 CDNTIMUF 

122 CJ5UM=0. 

123 D3 323 J=1,I 

12% I A B=JS { J ) 

125 CJ SUM=C JSUM+C ( I AB ) 

126 323 CONTINUE 

ISO JRME=I 

131 G3 TO 72 

132 325 CONTINUE 

13% JRME=NR( IMAX) 

135 G3 TO 100 

136 52 MS =MS+ 1 

C*** I TERATION COUNTER, M,\IT if msmax exceeded 

c 

137 IF (MS -MSMAX) 5%, 5%, 53 


1%Q 

53 

CONTINUE 

1%1 


IF ( JRME-39) 59, 59, j 1 

1%2 

51 

CDNTINUE 

1%3 

59 

K= 1 

1%% 


DD 35% 4=1, NN 

1%5 


IF (INVCCJ)) 35%,35%,:>52 

146 

352 

J 3 (K)=J 

1%7 


K= K + 1 

150 

35% 

CD NT I Si UC 

152 


IF CNN-39) 356,356,355 

153 

355 

CD NT I NIUE 

154 

356 

GD TO 993 

155 

5% 

DD 57 I Y = 1 , NN 

156 


1= IY 

157 


IF ( JS ( I ) ) 55, 58, 55 

160 

55 

IF (NN-1) 57 , 56, 57 

161 

56 

1=1 + 1 ' " ■ ■ ■ ■ , . : 

162 


GD TO' 58 ' ■ ■ ' ■ ■■ : 

163 

57 

CONTINUE : ' '■ :■ 

165 

58 

JRME=T~L 

166 


CJ SUM=0 ' " 

167 


DD 6% 1=1,, JRME ' ' , , 

170 


IF USTD) 6 %, 6 %, 62 

iT't. 

'62 

I A 8=45 Cl ) 

172 


CJ SUM = CJ SUM+CC I AB ) 



FORTRAN SOURCE LIST l HOI 


rt,.G2 50 

ISM SOURCE STATEMENT 


173 

64 

CONTINUE; 

175 


DO 70 1=1, MM 

176 


Y( I >=-Bt I ) 

17 7 


D3 68 J=l, JRM_ 

200 


IF (JS(J)I 6b, 68, 67 

201 

67 

14 B=JS { J ) 

202 


Y( I )=Y (I ) +A( I, I AB) 

203 

63 

CONTINUE 

205 

70 

CONTINUE 

207 


DO 71 I =1 , MM 

210 


1= ( Y( I ) ) 100, 71,71 

211 

71 

CONTINUE 

213 


IF (ZMIN-CJSUM) 73, 7 , 72 

214 

72 

ZM I N= CJSUM 

215 

73 

CONTINUE 

216 

83 

DO 76 J= 1 » JR ME 

217 

76 

JO PTC J )=JS( J ) 

221 


14 8= JR ME + 1 

222 


DO 79 J=I AB, NN 

223 

79 

JO PTC J ) = 0 

225 

74 

IF CJSUL(JRME)) 75, 92, 92 

226 

75 

DO 80 1=1, JRMfc 

227 


J= JRME+l-I 

230 


IF ( JSULt J ) ) 80, 85 , 85 

231 

80 

CONTINUE 


o ** 

P 

PROBLEM COMPLETE, WRITE OUTPUT 

233 

L 

K= 1 

234 


DO 84 J=l, NN 

235 


IF (INVCC J 3 ) 84, 84 , 3 2 

236 

82 

JP ( K ) = J 

237 


K= K + 1 

240 

84 

CONTINUE 

242 

88 

CONTINUE 

243 


DO 400 K=1,NN 

244 

400 

JP C K > = 0 

246 


DO 405 J = 1 , NN 

247 


1= ( JOPTC J) ) 401 , 405,402 

250 

401 

14 B=~J DPT ( J ) 

251 


J 3 ( 143 ) = JOPT ( J ) 

252 


GO TD 405 

253 

402 

14 B = JO PT ( J ) 

254 


J? { IAB)=JOPTC J) 

255 

405 

CONTINUE 

257 


K= 1 

260 


DO 412 J=1*NN 

261 


IF ( JPC J) ) 408,406,407 

262 

406 

IF ( INVCC J ) ) • jT 2 , 1 2 , - , 0 9 

263 

407 

IF (INVCC J ) > 409,4-0 9, 12 

264 

409 

JO PTC < ) = J 

.'2:65, 


K= K. + 1 

266 

.412 , 

CONTINUE 

270 

,aV , . 

IF (K-l-39) 413,413,414 

2?1 

,,413 

HD=K“1 

272 


GO TO 415 



G250 


FORTRAN SOURCE LIST Z-.RG1 

ISN 


SOURCE STATEMENT 

273 

414 

CONTINUE 

27 k 


I A B=K- 1 

2?5 

415 

Z3 PT = D 

276 


i \ B=K~ 1 

277 


DO 418 J = 1 , 1 A B 

.,00 


U BC = I ABS( JQPTU) ) 

301 


IF { I N VC ( 1 ABC 3 ) 417,417,416 

302 

416 

Z3PT = Z OPT - C ( I ABC ) 

303 


GD TO 418 

304 

417 

Z3PT = ZOPT + C ( I ABC! 

305 

418 

CONTINUE. 

307 


1= (MNMX) 420,420,419 

310 

419 

ZOPT=-ZQPT 

311 

420 

CONTI NUE 

312 


GO TO 993 

313 

83 

JSUH J )=-l 

314 


JS (J)=-JS( J) 

315 


1 3 C = J + 1 

316 


DO 90 K=IBC,JRMS 

317 


JS ( K ) = 0 

320 

90 

JS UL( K 3 = 0 

322 


GO TO 52 

323 

92 

JS ( JRMc)=-JS IJRMF) 

324 


JSULt JRME )=“1 

325 


GO TO 52 

326 

100 

DO 102 J = 1 , N N 

327 

102 

NS { J)=l 

331 


DO 106 J= 1 , JRME 

332 


IF ( JS { J ) 3 103, 106, 104 

333 

103 

JN EG-- JS ( J 3 

334 


NS ( JNEG) =0 

335 


GO TO 106 

336 

104 

I A B C= J S ( J 3 

337 


NS { IABC3 =0 

340 

106 

CONTINUE 

342 


DO 120 1=1, NN 

343 


1= (NS(I 3 ) 120, 120, 108 

344 

106 

IF (C( D-ZMI.N+: JSUH ) 120,110,110 

345 

113 

NS (13=0 

346 

123 

CONTINUE 

350 

121 

DO 122 1=1, NN 

351 

122 

ME ( I > = 1 

353 

' . ,7 

DO 140 1=1, MM 

354 


1= ( Y ( I ) 3 132,140, 140 

355 

132 

DO 135 J = 1 , N N 

356 

.* ’ ; 

IF (A( 1, J ) 3 135,135 ,134 

357 

134 

ME ( J ) =0 

360 

135 

CONTINUE 

362 

140 

CONTINUE 

364 


DO 150 J = 1 « NN 

365 


IF' (NS ( J 3 3 150, 1 50 , 14 

366'' 

144 

IF (ME ( J) 3 150,150, 14-; 

367 

146 

NS ( J 3 = 0 

370 : 

150 

CONTINUE 

372 

|J1| 

DO 160 J=1 , NN 



M . G2 50 


FORTRAN SOURCE LIST Z'RCJl 

ISM 


SOURCE STATEMENT 

•; 75 


1= (MSI J ) ) 160, 160, 167 

374 

160 

CONTINUE 

i 76 


GO TO 74 

377 

167 

M6 RKF = 0 

400 


DO 182 1=1, MM 

401 


1= ( Y ( I ) ) 169, 133 , 183 

402 

169 

SJ MCR = 0 

403 


AP OS = 0 

404 


1 A 6C= MR ( 1) 

405 


DO 160 N = 1 » I ABC 

406 


IIB = JFU I,N) 

407 


IF ( MS ( I AB ) ) 130,180,171 

410 

171 

HB=JR(I ,N) 

411 


SJMCR=SUMCR+CC I AB) 

412 


A? OS = A POS + A ( I , I AB ) 

413 


IF ( SUMCR-ZMIM+CJSUR) 176,177,177 

414 

176 

IF (APOS+Y (1)3 150, 131,183 

415 

177 

1= ( APOS+Y ( I ) ) 74, 74,183 

41 6 

180 

CONTINUE 

420 


GO TO 74 

421 

181 

1= (N-NR(I)) 173,182,182 

422 

173 

II BC=NR( I) 

423 


DO 174 K=IM » I ABC 

42 4 


IIB=JR(I ,K) 

425 


IF ( MS ( I AB ) ) 174, 174,183 

426 

174 

CONTIMUE 

430 

182 

Ml RKF= 1 

431 

133 

CONTIMUE 

433 


1= ( MARK F) 190, 190, 240 

434 

190 

DO 210 J = 1 , NN 

435 


IF (MSI J) ) 210,210, 192 

436 

192 

YISUM=Q 

437 


DO 198 1 = 1, MM 

440 


IF ( Y ( I ) +A ( I , J ) ) 193,193,198 

441 

192 

Yl SUM= YASUM+ Y ( I )+A( I, J) 

442 

198 

CONTIMUE 

444 


V( J ) =YASUM 

445 

210 

CONTIMUE 

447 


DO 230 I Y= 1, NN 

450 


J = I Y 

451 


IF (MSI J ) ) 230,230,218 

452 

218 

IF ( NN- J ) 235,235,219 

453 

219 

L= J + l 

454 


DO 225 K=L,NN 

455 


IF C MS ( K ) ) 226,226,220 

456 

220 

IF (V(J)-V(K)) 230,226,226 

457 

225 

CONTIMUE 

461 


DO 229 I =L, NN 

462 


IF (NS ( I ) ) 229, 229,22: 

463 

225 

IF ( V ( J ) - V £ 1 )) 227,227,229 

46 : 4 
■ 46 5' 

22? 

IF (CU)-C(I)) 229,229,223 

223 

J= I 

466 

229 

CONTINUE 

470 


GO TO 235 

471 

: 230. 

CONTINUE 








M...G250 


FORTRAN SOURCE LIST Z -R'jl 

ISN 


SOURCE STATE Mi- NT 

i-,Tl 

2S5 

JS( JRMc+l)=J 

4 74 


GO TO 52 

475 

240 

03 24l J»1,N!N- 

476 

241 

Mr ( J ) = 0 

500 


DO 256 1=1 k MM 

501 


1= ( V ( 1 ) ) 243,256,256 

502 

243 

A? QS=G 

50 3 


DO 248 J = 1 , NN 

504 


IF { NS ( J ) ) 248,248,245 

505 

245 

1= ( A ( I , J ) ) 243 , 248 ,2s 6 

506 

245 

A? OS = APOS + A( I , J ) 

507 

248 

CONTINUE 

511 


IF (APOS + Y ( I ) ) 256, 249, 256 

512 

249 

DO 255 K=1,NN 

513 


IF t NS ( K ) ) 255,255, 250 

514 

25D 

1= { A { I , K 5 ) 255,255 ,231 

515 

251 

MF { K ) = 1 

516 

255 

CONTINUE 

520 

255 

CONTINUE 

522 


CFSUM=0 

523 


DO 262 1=1, NN 

524 


1= (MFU )} 262,262,260 

525 

260 

CF SUM=CFSUM+C{ I ) 

526 

262 

continue 

530 


1= (CJSUM+CFSUM-ZMIN) 263,74, 74 

531 

263 

IF (MS-1 ) 264,264,265 

532 

264 

J= 1 

533 


GO TO 267 

534 

2 65 

J = JRME + 1 

535 

267 

DO 263 1 = 1, NN 

536 


IF (MFU)) 268,268,266 

537 

266 

JS ( J)=I 

540 


j = j + 1 

541 

268 

CONTINUE 

543 


GO TO 52 

544 

9 93 

RETURN 

545 


END 

M6G250 


IBMAP ASSEMBLY ZER01 

MO MESSAGES 

FOR ABOVE ASSEMBLY 



;%„G2 5D 


FORTRAN SOURCE LIST 

I SN 


SOURCE STATEMENT 

0 

$1 BF 1 

X 3ERQ 

1 


SJBROJTINF GcRO ( NN, MM, MNMX , I NEQ , A, B, C , I NVC jMSMAK ) 

2 


DIMENSION INIQC 10), AC 10,4Q), BUG), CC40) , INVCC40) 

3 


INTEGER A, B,C,SUMA 

4 


MM M = MM 

5 


1= ( MM MX ) 2,4,2 

6 

2 

03 3 J = 1 , H N 

7 

5 

C( J)=-C{ J) 

11 

4 

D3 25 J=l, MM 

12 


IF { 1 N E Q { J ) - 1 ! 25,25,5 

13 

5 

IF CIMLQ(J) — 2) 25,8,6 

14 

5 

1= ( I N c Q ( J ) ~ ) 25,15,15 

15 

• 8 

DO 10 K= 1 , NN 

16 

ID 

A( J ,K) = ~ A ( J, K ) 

20 


B( J ) =- B ( J ) 

21 


GO TO 25 

22 

15 

MMM=MM+1 

23 


DO 20 K= 1, NN 

24 

23 

AC MMM, K) =“ At J, K ) 

26 


B( MMM) =-B ( J ) 

27 

25 

CONTINUE 

31 


MM =MMM 

32 


DO 30 J = 1 , NN 

33 


IF ( C ( J ) ) 27,30,30 

34 

27 

C( J ) = - C ( J ) 

35 


IN VCC J ) = 1 

36 


DO 28 1=1, MM 

37 

28 

AC I , J) =~A ( I, J) 

41 

33 

CONTINUE 

43 


DO 35 1=1, MM 

44 


SJMA=0 

45 


DO 34 J = 1 , NN 

46 


IF ( IN VC ( J ) ) 34,34,32 

47 

32 

SJMA=SUMA+AC I , J) 

50 

34 

CONTINUE 

.52 


BC I ) =B ( I ) +SUMA 

53 

35 

CONTINUE 

55 


RETURN 

56 


END 

M^G250 


IBMAP ASSEMBLY GERO 

NO MESSAGE! 

i FOR ABOVE ASSEMBLY 



FORTRAN SOURCE LIST 


M..G2 50 


ISN S3 URGE ST AT tMLMT 


0 $IBFTC 3ISCAL 

1 S J BRDUTI NE DISC AL { TD) 

2 C3MMOSl/STS/f4V*NN 

3’ C3MMON/OPM/B»T» JMX» KRT 

4 INTEGER TD,B,T 

5 DIMENSION JM/. ( 40, 40 ) , o C 40, 40 ) , KRT 1 10 , 30 ) 

6 T3 =0 

7 D3 I 1=1, NN 

10 D3 2 J = 1 , NN 

11 IF ( JMX ( I * J ) • LQ» 1 ) GO TO 3 

1A 2 CONTI NUE 

16 3 J1=J 

17 TD=TD+B( I, J1 ) 

20 1 CONTINUE 

22 RETURN 

23 END 


Me, G2 50 


IB MAP ASSEMBLY OISCAL 


NO MESSAGES FOR ABOVE ASSEMBLY 



XL G230 

ISN 

0 

1 

2 

3 

4 

5 

6 
7 

10 

11 

12 

13 

14 

15 

16 

17 

20 

21 

22 

23 

26 

27 

30 

31 

34 

35 
40 

43 

44 

45 

46 

47 
50 
53 
56 
60 
61 
62 
65 
67 

70 

71 

72 

73 

74 

75 

76 

77 
102 

105 

106 
107 
110 
111 


FORTRAN SOURCE LIST 

SOURCE STATEMENT 

$1 BFTC REFINE 

SUBROUTINE R6.FIN£(N) 

C*****p*oGRAM FOR 10 7 EH I CL S AND 30 CUSTOMERS 

DI MENS ION N ( 1 0 ) , D ( 4 0 } , B { 40 , 40 ) , KRT ( 10 , 30 ) * MMX (40,*-0) , J MX ( , - 0 5 
D[ MENS I ON ND I S { 40 ) , LI { 40 ) , L2 < 40 ) »L (40 ) » L3 < 40 ) T OVC ( 10 ! , UVC ( 1 ) 

DI MENS 1 0 H MC U S ( 3 0 ) , MN U S ( 3 0 ) , M RO UT ( 30 , 1 0 ) , H W 1 ( ? 0 ) , N wZ C 0 ) * H K 0 1 S C 
DIMENSION L4 ( 40 ) 

01 MENS I DM NS 1 ( 30 ) » N S2( 40 ) , NSS 1( 30) ,NSS2( 20 ) , N£01 S ( SG ) 

01 MENS ION C( 10) 

DIMENSION N20T(40) 

INTEGER C 

INTEGER TD,D,b,OVC f LVC jT 
C3 MM DM / 3 P M / B ? T » J MX , KRT 
COMMDN/STS/NVjNN 
C3MM0N/JLM/D,C 
CALL OISCAL(TD) 

ITD=TO 

03 10 1=1} NN 
03 10 J = 1 » NN 
MMX ( I » J ) = J MX ( I t J ) 

10 C3NTIMUE 

03 103 1*1, NV 
KS P = N( 1 ) 

KS PP=KSP 

IF (KSP.EQ* 1)G0 TO ICO 
D3 101 J = 1 j KSP 
IF C N 1 1 ) . EQ. 1 ) GO TO 10G 
IF ( KSPP«LT « J ) GO TO 100 

Kf =KRT ( I , J ) 

K3 = 0 

03 102 J J = 1 * NN 
03 103 KK= 1, KSPP 
KTT=KRT ( I »KK ) 

IF € JJ.EQ.I ) GO TO 102 
1= ( JJ. fcQ.KTDGD TO 102 

103 C3NT iNUE 

I D I S=B { J J, KT ) 

03 104 JK= 1 f NN 

1= ( JMX( JJ, JK).EQ.1)G0 TO 105 

104 CONTINUE 

105 J3 I S=B { KT, JK) 

K3 I S=B ( J J, JK J 
K3 =KQ+ 1 . 

MDI St <Q) -IDI S+JDI S-KDIS 

LI ( KQ } = J J 

L2 C KQ) = JK 

L( KQ ) = KT 

03 106 JJK = 1*!MN 

IF ( JMX( JJK,K7).£Q.l 160 TO 107 

I c ( JMX(KT»JJK). :.0.1)GO TO 108 

G3 TO 106 

107 L3 ( KQ) = J JK. 

■ G 3 ' TO 106 

108 L4 ( KQ) = JJK 

106 C3JNTI NUE 



M;. G250 

ISN 

113 

115 

116 
11 ? 
120 
121 
122 

123 

124 

125 

126 
127 
130 
13 1 
132 
133 
13 % 
13 5 
136 
141 
14% 
1%6 
150 


151 

152 
153 
15% 
155 
160 
161 
162 
163 
166 
167 

170 

171 


173 

174 
177 
200 
201 
202 

203 

204 

205 
210 

213 

214 


SOURCE ST AT . M, NT 


FORTRAN SOURCE LIST R. FIT 


102 CONTINUE 

CUL MIN(NDIS|<Q»KQQ|IMIN) 

N 3 S=L ( KQQ ) 

N 3 1 = L 1 { KQQ ) 
r-P2 = L2(KQQ} 

N 3 3 = L3 { KQQ ) 

N? 4=L4 ( K QQ ) 

JNX( NP3, NFS) =0 
JMX { NPS, NP4 ) =0 
JMX ( MP 3 , NPA ) = 1 
JMX ( NP 1 » NP5 ) =1 
JMX( NPS» NP2)=1 
JNX ( NP 1 s NP2 ) =0 
CUL 3 U D I Y A ( JMsCjKRT » N ) 

DO 109 I 1 = 1, NV 
N<ST=N( 1 1 ) 

UD 110 I K = 1 , NKS T 
KMT=KRT ( II, IK) 

I- (Lit KQQ ) . l Q . K MT ) G 0 TO 111 
1= {L2(KQQ)*EQ,KMT)GC TO 111 

110 CONTINUE 
109 CONTINUE 

1 11 INR T = I I 

C ***** N > S IS THE CUSTOMER Oa ROUTE 1 WHICH IS TO BE SHIFTED TU SOME OTl- 
C ***** N ;> 1 AND HP 2 ARC NEW CONNECTIONS 

C *****nrT IS THE ROUTE TO WHICH CUSTOMER L ( KQQ ) IS SHIFTED 
KT Z=N( II ) 

CALL DAPCONI KRT , N»0 VC, UVC ) 

7539 FORMAT! //5X» *SUBTQU F NO. * ,5X , *OVC { . ) *» 5X, *UVC ( . ) *5X , *PATH S £0.*/) 
7339 FORMAT (10X, 12, 8X, 15 ,6X, 13, 8X, 3013) 

IFfOV: ( II ) .EU.OJGO TO 4889 
K3 =0 

DO 112 I R = 1, KTZ 
KST = KRT ( 1 1 , 1 R ) 

IF ( D( KCT) eGT.OVC { 1 1 ) ) 60 TO 113 
GO TO 112 
113 K3=KB+ 1 

MS US ( KB ) =KCT 

112 CONTINUE 

C FOR ACCOMODATING NPS -IN IT) MCUSU ) IS ACTUAL CUSTOMER ON ROUTE 

C*****k3 IS THE NO. OF CUSTOMERS PERMISSIBLE FOR SHIFTING FROM ROUTE I 
C** NOW CHECK WHICH CUSTOMERS OUT OF ABOVE MCUSU) CAN BE SHIFTED TO 
C OTHER ROUTES AND FOR EACH OF THESE CUSTOMERS, NOTE THE POSSIBLE 

KB B=0 

1= { KB. LQ.03G0 TO 5879 
K3R=0 

DO 114 IU=1» KB 

m:ust=mcus< iu) 

KB B=K8B + 1 
KRi =0 

DO 115 IT=1,NV 
IF(IT.EQ.II)GO TO 115 
1= (UVC(IT) . G . . 0 { MCUST ) ) GO TO 116 
GO TO' LI 5 
116 K3L=KSL+1 



0 


FORTRAN SQURC; 

LIST R FIN 

ISM 


S3 URGE STATS ML NT 


215 


IF {K8L.GE.23GQ TO 1 13-9 


220 


K3R=KBR+1 


221 


NZOKKBR ) = KBL 


222 


MN US { KBR 3 = MCUS T 


223 

1139 

MROLT(KBR s KBL)=IT 


224 

115 

CONTINUE 


226 

1 14 

CONTI N Ut 



0***4 

KER IS THE NO. OF C1ST0MERS ON ROUTS II 

WHI CH CAN Be 


C 

OTHER ROUTES AND DEL TI ON OF WHICH CAN ACCOMODATE NFS 

230 


IF (KBR.EQ.O)GD TO 5 879 r 



o** 

ABOVE STATEMENT MEANS NO CUSTOMER CAN BE 

DELETED FROM 


c 

SO CONSIDER NEXT CUSTOMER ON ROUTE I. 


233 


DO 117 I Z = 1 j KBR 


2 34 


MS C = MN US { IZ) 


235 


03 118 I WK= 1 f N!N 


236 


IF (JMX(IWKfMSC) . Q. 1) GO TO 119 


241 


GO TO 118 


242 

1 19 

NW1(IZ)=IWK 


243 

118 

CONTINUE 


245 


DO 111S I WKQ = ly NN 


246 


IF ( JMX(MSCslWKQ) .EG.13G0 TO 120 


251 


30 TO 1118 


252 

123 

NW2(IZ)=IWKQ 


253 

1 118 

CONTINUE 


255 

117 

CONTINUE 


257 


DO 121 IZ = 1» KBR 


260 


MSC = MNUS( IZ) 


261 


N 3 LM=NZOT ( IZ ) 


262 


KC ONT = 0 


263 


DO 122 I ZK= 1 » NN 


264 


DO 123 1UK*=1| NPLM 


265 


KROUT=MROUT(IZ, IUK) 


266 


NJW=N( K ROUT 3 


267 


DO 124 JUK»1, NOW 


270 


KARR=KRT(KROUT, JUKI 


271 


IF ( IZK.EQ. KARR) GO TO 125 


27 4 


GO TO 124 


275 

125 

M< D I S= B ( KARR » MSG 3 


276 


DO 126 I WKK = 1 1 NN 


277 


IF { JMK{ KARR, IWKK) .EQ.13G0 TO 127 


302 

125 

CONTINUE 


304 

127 

J<DIS=B(MSC, IWKK) 


305 


K< DI S = B « KARR j I WKK 3 


306 


k:ont=kcqnt+i - •••' ; 


307 


NKDI S ( KCOINT 3 =MKDI S+ JKOI S-KKDI S 


310 


NS 1 ( KCONT ) -KARR 


311 


NS 2 { KCONT ) =1 WKK 


312 

124 

CONTINUE 


314 

123 

CONTINUE 


316 

122 

CONTINUE 


320 


CALL M I N { N K D Z S , K C 0 N T , K C £ M T , J I N 3 


321 


NS S 1 ( I Z 3 =M S 1 { KC ::N1 ) 


322 


NS S2 ( I Z ) =NS2 ( KC : NT 3 


323 


NATC=NW1 (IZ) 


324 


MAT0=NW2 ( IZ). 



SHIFT 

FROM 

R OUT 



FORTRAN 


I SN 


SOURCE STATEMENT 

325 


MOP IS = B{ NATO, MSC > 

326 


NOPIS=B(MSC|MATO) 

327 


KO P IS= B ( NATO j MATO ) 

330 


MI 01 S = MOP IS + NQP I S-KCP1 S 

3 s 1 


NI 0 IS( 1 Z)=JX N-MZDIS 

33 2 

121 

CONTINUE; 

334 


ULL MINtNZDIS, KBR, KBRR, J INE ) 

3 3 5 


M7 =MMUS ( KBRR ) 

3 3 6 


Mf 1 = N/U(KBRR) 

337 


NY 2=NW 2 ( KBRR ) 

340 


NY 3»NSS1( KBRR ) 

341 


NY4=NS S2 I KBRR ) 

3 42 


J1X(NY1*NY2) = 1 

343 


JMX ( NY 1 j NY ) = C 

344 


JMX ( NY t NY2 ) = '? 

345 


J«tX(NY2 »NY ) = 1 

346 


JMX(NY ,NY4 ) = 1 

347 


JMXINY3 »NY4) =0 

350 


CALL GUOIYA( JMX»KRT »N) 

351 


CALL C A PCOM ( KRT » N » 0 VC » UVC ) 

352 

4839 

CONTINUE 

353 


CALL OISCAUTD) 

354 


NT D=TD 

355 


IF ( NTD. LT . ITD )G0 TO 160 

360 

5879 

DO 140 I£S=1 »NN 

361 


DO 140 JES = 1 » NN 

362 


JMX( IES { JIiS)*MMX( ISSi JES) 

3 63 

143 

CONTINUE 

366 


GO TO 150 

367 

160 

DO 200 I FZ=1 »NN 

370 


DO 200 1 £ W = 1 j f 1 N 

371 


MNXdFZ, IEW)=JMX( IFZ,1EW) 

372 

200 

CONTINUE 

375 


I T D=NTD 

376 

150 

CALL SUDlYAt JMX,KRT»N) 

377 


CALL CAPCON{KRT,N,OVC»UVCJ 

400 


KSPP=N(I ) 

401 

101 

CONTINUE 

403 

100. 

CONTINUE 

405 


CALL GUDI YA ( JMX* KRT »N } 

406 


CALL C APCON ( KRT » N j 0 VC j UVC ) 

407 


RETURN 

410 


END 


M.,G250 


IBMAP ASSEMBLY REFINE 


NO MESSAGES FDR ABOVE ASSEMBLY 
MEG250 


IBLDR 


JOB OOOOOO 




*** OBJECT PROGRAM IS BEING ENTERED INTO STORAGE AT II HRS. 2? M 


ftSSlSBiSlW^ 









»#**#♦*»*.•#•*** #***#*#*#**********•**»***■*■*•»******* ****** 
* * 

* SHORTEST PATH ALGORITHM * 

* * 
#******•*«■**♦* ******************************** ********* 


DIMENSION DO C,3C ) 

INTEGER 0 
* P«A0 ICO , N 
CC FORMAT ( IB ) 

N1 N=N~1 

DO 101 I = 1 * N I N 

JJ=I+i 

REAL 102*ID(I,II ),II=JJ,NI 
101 CONTI NLF 
:: 02 FORMAT (2014) 

CG 2 GO 1=1, N 
DO ICO J-l , N 
IF (J.GT.I)GC TO 200 
IFl J.EQ.IJGG TO 201 
CCI, JJ=D( J»I) 

GO TO 200 
201 D ( I , J )=0 
200 CONTINUE 

DO 5CI0 I TR=1 *8 
PRINT 400,(1, 1=1, N) 

400 F0RMAT(iHl,////25X,*DI STANCE NAT RIX*//11X, 3014/ > 

'DO 1.03 1=1, N 

PRINT 104, I*(C(I,J),J=1,N) 

103 CONTINUE 

104 FORMAT ( 5X, 12, 5X,30I4> 

PRINT 500 

5CC FORMAT (//10X,*N0TE~ D(I ,J) MEANS SHORTEST CIRgCT LINK BETWEEN 

61 AND J NODE* ,/10X,* DC I, J) =999 MEANS NO DIRECT LINK BETWEEN 

, 61 AND J NODE* > 

5010 CONTINUE 

CALL FLOYD ( D, N! 

CO 6010 JTR=i ,6 
‘ PRINT 300, (I, 1=1, Nf 

300 FORMAT ( 1H1 , // //25X,*SH0RTEST DISTANCE MATRIX*//! IX, 30 14/ I 
DO 105 1=1, N 

PRINT 106,1, (0(1, J),J='UN) 

1.05 CONTINUE 

106 FORMAT ( 5X, 12, 5X,3CI4) 

6010 CONTINUE 
PRINT 700 

700 FORMAT ( iHl, // *FINISH* ) 

'■ • , STOP 

■ ' gNC : 

f IBFTC FLOYD 

SUBROUTINE FLCYD(D,N) 

INTEGER D, S 
DIMENSION DC*IC,30) 

DIMENSION I PC 30,30 ) 

" DO' 1.5 1=1, N 
DO IS J®1 , N 
1 : . IPCI,J)=0 

. 15' ' CONTINUE 



0 0 . X ~ f N 

!FU{I,K). J *C.S99)G'J TO ’ 2 
CO 11 J=V » N 

IF (DtK* J).EC. S 99 ) GO TO , 1 
S=C(I,K)+D(K, J) 

IF (S.LT.DC I»J ) )GG TO FO 
GO TO 1 1 
5G DII,J)=S 
IP (I * J )~K 

11 continue 

12 CONTINUE 
RETURN 

EMC 



u u 


* •* 4 4?*#*****#* **#■* •♦♦♦♦■•••If**#***-*#*#***-**#'*#**#********# 

* * 

* CLARKS' AND WRIGHTS METHOD * 

c * 

#*•* ********** **#**#**« #?##***#*##.»****##«***#***■****** 


frC******« THIS PROGRAM IS FOR 34 CUSTOMERS AND l DEPOT 

D I F INSIQN DO 3 > , B i 35 , " 3 ) ,C( 15 ) , KCP (35 } » JMX (35,35), IB (35 ),JB(35) 
CINfNSICN DROUTU5) ,Kfi ( • 5) , KI P (35 ) , JGO ( 35 ) 

DIMENSION KRT (13*35)»N( 5) 

INTEGER 8, D» OROUT » C 
COFFCN/SIT &/NT,NN« JMY * ° 

***** N C IS NU* " CF CUSTOMERS AND NT IS TOTAL NO . CF VEHICLES AVAIL6BL 
WHICH FIGHT Bfc LESS OR FCR5 THAN ACTUAL VEHICLES NlrSf Df-D 
RE AC I, NC* NT 
? FORMAT ( 213 ) 

IF (NC.EQ. iOO) GO TO 9000 
NN=NC+I 

C****«**FGR DEPOT DEMAND IS ZERO, D(i)=0 
CO 3 1=1, NC 

JJ=I+i 

R'iAC 4, (B( I, I I), II=JJ,NN) 

3 CONTINUE 

4 FORMAT ( 2014 ) 

CO 1020 1=1, NN 
CO 1C2C J=1,NN 

IF ( J.GT. I ) GC TO 1020 
I F ( J.EQ . I ) GC TO 1021 
B(IfJ)=B{J,I) 

GO TO 1020 
1021 e(I, J>=999 
1020 CONTINUE 
IPO 1=0 

8049 IP0L=IP0L+1 

IF ( IPCL.GT . 4) GO TO 9000 
‘ . PRINT 250 

25 0 FORMAT ( IH1 » f\ CX,*ROUT£S OBTAINED BY CLARKE AND WRIGHT METHOD* ) 
READ 2, ( D( J ) , J=1 ,NN) 

1 2 FORMAT (2014) 

JPQLCY=0 

8 JPOLCY= J POLCY +1 
PRINT 5050 

5050 FORMAT (5X,33(* * *) ) 

IF t JPOLCY. Gt, 3) GO TO 8049 
1 060 READ 1001, (C( I), 1=1, NT) 

1001 FORMAT ( 1016 ) 

PRINT !5 ! :<0t JPCLC Y 

1550 FORMAT ( h OX,* CONSIDER INC Vf.H ICL> P0LICY»,I2 » 

PRINT 2005,(0 (I ) ,1=1, NT) 

2 0 C; FOR MAT (, f OX,*CAPACITI’"S«»I 01 5 / ) 




“inn Tfin ) aqqoq 

DO : C I«:,KN 
CO 103 JJK-% NN 
IF (KCP(JJK) . • C. I )GD TO 03 
102 CONTINUE 
GO TO 1 04 
: *03 18(1 ) =9999 

GO TO 10 

jO As 18 ( I ) = Bllt I ) 

1C CONTINUE 

CALL MIN( IB,NN»NK,IBB) 

IF (IB8.EC.9S99) GO TO 4CG0 

C*-***-» NK IS THE NEAREST CUSTOMER TO DEPOT AMONG UNAS SIGNED CUSTOMERS 
KR0UT = KR0UT+I 
IF (KR0liT.G7.NT )GC TO 40..Q 
JMX ( 1 »NK 1 =1 
DROUT ( KROUT )=D ( NK ) 

500 K=K+1 


KOPCK)=NK 

C**«* JB(I) IS SAVINGS OBTAINED BY CONNECTING NK TO 
JB Cl >=*’9999 
«JI E8=B ( X, NK ) 

DO 11 1=2, NN 
DO 12 KK=1,K 

IFtl.EQ.KOP(KK) )GO TO 13 

12 CONTINUE 

JB ( I 1 = JlBB+Bl I,1)-8CNK,I) 

GO TC 11 

13 OB (I )=*»9999 
Xi CONTINUE 

DO 107 1=1, NN 
107 KB(I)=~J5C I) 

MKP = C 

16 CALL PIN(KB,NN,NKK,K88) 

IFCKBB.6Q.999S)G0 TO 23 



IFCKEB.EQ. S99)G0 TO 15 

MKP=MKP+I 

KIPC KKP)=NKK 

KB(NKK>=999 

GO TO 16 

CONTINUE 

CO 17 1=1, PKP 

NKT=KIP(I> 

JOCtI ) = JBl NKT ) 

CONTINUE 
CO 20 1*1, m? 

NKT=Kl PC 1 1 

DR OUT C KR OUT ) = DROUT C KR OU?)+D{ NKT ) 
IFCCRCUTt KROUT ) .GT.C ( KROUT) ) GC TO 21 
GO T 0 22 

DROUT ( KROUT ) = DROUT f KROUT KDC NKT ) 
CONTINUE 








I 




f 's ) =<JpOQ 


CO ’01 JJK= , NN 

IF{KCF( JJK).- C.DGO TO ' 03 

2 CONTINUE 
gc t c :c ! 

3 T P ( J ) =9999 
GO TC ;io 

4 IECI ) = B t i » I) 

C CONTINUE 

CALL PIN{ IByMNfNKflBB) 

IF(lB6.gC.999S)G0 TO 4000 

*** N* IS the NEAREST CUSTOMER TC DEPOT APONG UNASSIGNED CUSTOMERS 
KROUT=KR OUT+1 
IF (KROUT.GT.NT JGG TO 4050 
JMX«ltNK)=X 
DROUT ( KR OUT ) = D (NK ) 

C K = K+ 1 

KCP ( K) =NK 

** JBCI) IS SAVINGS OBTAINED BY CONNECTING NK TG I 
JB ( 2 )=~9999 
jiee=Bti,NK» 

DO II 1=2, KN 
DO 12 KK=1 , K 
IFCI-EC.KOPCKKHGO TO 1* 

2 CONTINUE 

JB ( I J= JI8B+B( J,U- B(NK,I) 

GO TC II 

3 Jeill=**9999 
I CONTINUE 

00 IC7 1=1, NN 
? KBCn=«JB( I) 

NKP = C 

6 CALL PlNt K8,MN,NKK, KBB ) 

IFCKBB*EQ»9999)GG TO 23 
IF (KE8. EQ. 999 ) GO TC 15 
MKP = PKP + .. 

K I P { PK P ) =NKK 
KB (NKK } =999 
GO TO 16 

5 CONTINUE 
CO 17 1=1, NKP 
NKT = KIP { I ) 

: JCCII >=J6(NKT ) 

7 CONTINUE 
CQ 20 1= , PKP 
NKT=Kl P ( I ) 

DR CUT I KR OUT ) = DROUT C KR OU )+D( NKT ) 

IFtCRCUTl KROUTJ.GT.C(KROUT) JGC TO 21 
'GO TO 

:i DROUT ( KR3UT ) = DRDUT (KRGUT I« DINKT) 

!0 CONTINUE 

' feft ' ' ft ■ '%% ' 

- , wW it 1 # m 

ft'' JlKlWf NKTi=l 











DO -10 I = i » KN 
N P S = N ( I ) 

PR t 20* I» I I» l KRT { I , J),j=l,NpS ),H 

FQ- K AT € I 0X,-*SUBT0UR*tI3,5X»30I3 j 
IC CONTINUE 

»rint 2eo, IOIS 

'0 F3Ry.AU/ OX,* TOTAL DI STANCE =*,I4) 

PRIM ?20tKN 

c FORMAT { 10 X»*REQD. VEHICLES =*, 14 ) 

GC TC 8 
)50 PRINT 290 

JC F0RWATt//20X,*RECD. VEHICLES? ROUTES ) ARE FORE THAN AVAILABLE VFH*) 

GO TC 3 
ICC STOP 
END 

I FTC BA £80 

SUBROUTINE BABBG (KRT, N t I CIS* KM 
COPMON/SlTA/NT«NNt JMX V B 
INTEGER B 

DIMENSION 8? 35,35),JM>'«'5 f 35),NU5>,KRTU3,35) 

I-C 
ID I S=0 
D0^2 4-1, NN 

IFIJMXU,JJ.EC.1)G0 TO 3 

1 = 1 + 1 

ID IS-ICIS+B (1 ,J2) 

5 J1=J2 

K=K + 1 

KRT ( I , K) =J1 
DO 130 JJ=l,NN 
XF{ JWXt Ji, JJJ .EC.DGO TC 4 
C CONTINUE 

4 J 2 — J J 

I0IS=IDIS+B(Jl,J2) 

IFtJ2.EO.DGO TO 6 

« *T 1*1 S 

i y zt 

2 CC MINUS 
KN = I 

RETURN 

END 

JFTC WIN 

SUEPGUTI 'F WlMP,N,K,RWiN} 

INTEGER P» RHI NfSUM 

Dlftf Jf'SION PtSCi 

k = :; 

■■ ' su'P»PTii / ' ' ' . 

-"V/:' IFlN.iQ.lJGC TO '2 




/• 31705 ! 






